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Résumé 

De nombreux problèmes combinatoires peuvent être formulés comme des Problèmes de Satis¬ 
faction de Contraintes Valuées (VCSPs) où les contraintes sont définies à l'aide de fonctions coût 
qui reflètent différents degrés de satisfaction. Résoudre un VCSP revient à trouver une attribution 
de valeurs aux variables ayant un coût global optimal. Pour plusieurs problèmes, cette tâche est 
difficile du point de vue de la complexité algorithmique puisqu'il s'agit d'un problème NP-difficile. 

Cette thèse présente trois contributions qui consistent en des techniques visant à améliorer la 
résolution des Problèmes de Satisfaction de Contraintes Valuées : 

- La première contribution consiste en une technique de résolution qui s'appuie sur la notion 
de substituabilité directionnelle. Cette dernière notion a prouvé son efficacité dans la réso¬ 
lution des Problèmes de Satisfaction de Contraintes CSPs classiques [91]. Pour une classe 
de VCSPs binaires, où les domaines de valeurs sont réduits à des chaînes de valeurs dites 
directionnellement substituables, un algorithme d’identification de complexité 0(nd 3 ) 1 et 
un algorithme de résolution de complexité 0(nd ) ont été proposés, ce qui fait de cette 
classe une classe polynomiale [92], 

En s'appuyant sur cette classe, nous avons également proposé une méthode de décomposi¬ 
tion de VCSP, dont les valeurs ne sont pas forcement directionnellement substituables, en 
un ensemble de sous-problèmes polynomiaux. 

- La deuxième technique de résolution s'appuie sur la classe des fonctions modulaires [45]. 
Pour une classe de VCSPs binaires où toutes les fonctions de valuation sont exclusivement 
modulaires, un algorithme d’identification de complexité 0(ed 2 ) 2 et un algorithme de ré¬ 
solution de complexité 0(ed ) ont été proposés [93, 58]. 

Nous avons également proposé une méthode de décomposition de VCSP s'appuyant sur 
cette classe polynomiale. 

- La troisième technique de résolution s'appuie sur les fonctions sous-modulaires [14]. Une 
super-classe des VCSPs sous-modulaires est obtenue en calculant, en temps polynomial, des 
permutations des domaines des variables d’un VCSP quelconque en vue de le transformer en 

1. Notons que : 

- n est le nombre de variable d’un VCSP ; 

- d est le nombre de valeur du domaine le plus large dans un VCSP. 

2. e est le nombre de contrainte dans un VCSP. 
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un VCSP sous-modulaire. Pour les VCSPs de cette super-classe, que nous appelons VCSPs 
sous-modulaires par permutation, une solution algorithmique de complexité 0(n 2 d 4 ) a été 
proposée [56, 57, 59]. 

Nous avons également proposé une méthode de décomposition de VCSPs binaires quel¬ 
conques s’appuyant sur la classe de VCSPs sous-modulaires par permutations. 
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Les problèmes de satisfaction de contraintes ou CSPs (Constraints Satisfaction 
Problems) constituent un cadre simple et assez général pour modéliser divers pro¬ 
blèmes combinatoires. Un CSP se définit par la donnée d’un ensemble de variables, 
chaque variable pouvant prendre une valeur parmi un ensemble de valeurs possibles 
appelé domaine. Un CSP est aussi défini par un ensemble de contraintes exprimant 
des restrictions sur les combinaisons de valeurs que peuvent simultanément prendre 
les variables. 

Résoudre un CSP revient à affecter aux variables, des valeurs de leurs domaines 
respectifs de telle sorte que toutes les contraintes soient satisfaites. Cette tâche est 
difficile du point de vue de la complexité algorithmique puisqu’il s’agit de résoudre 
un problème NP-complet [82], 

En général, un CSP peut avoir plusieurs solutions, comme il peut ne pas en 
avoir du tout. Pour les CSPs dits classiques, on cherche une solution qui satisfait 
obligatoirement toutes les contraintes (solution réalisable). Si plusieurs solutions se 
présentent, alors elles sont toutes supposées avoir la même “qualité”. De ce fait, la 
résolution de tels problèmes se résume souvent à trouver une seule solution réalisable. 

Cependant, dans plusieurs situations réelles, telles que h ordonnancement dans le 
secteur industriel [16], les solutions ne peuvent pas être considérées comme équiva¬ 
lentes au sens que certaines d’entre elles sont jugées meilleures que d’autres selon 
un critère bien déterminé. Il s’agit désormais de trouver la solution qui satisfait les 
contraintes et qui, en plus, optimise un critère mesurant la qualité des solutions. 
Les problèmes ainsi définis sont alors des problèmes de satisfaction de contraintes 
flexibles ou Soft CSPs (SCSP). Dans les SCSPs, les contraintes utilisées sont dites 
flexibles par opposition aux contraintes dures que l’ont trouve dans les CSPs clas¬ 
siques. 

Les CSPs flexibles fournissent un cadre plus général que celui des CSPs clas¬ 
siques permettant ainsi de résoudre plus de problèmes combinatoires. En effet, dans 
plusieurs cas réels, on a besoin d’exprimer divers degrés de satisfaction pour refléter 
la spécificité du problème. 

Le cadre des Problèmes de Satisfaction de Contraintes Valuées ou VCSPs (Va- 
lued Constraint Satisfaction Problcm) [106] permet de bien prendre en compte les 
subtilités d’une telle situation. Les VCSPs sont des CSPs flexibles qui se distinguent 
par le fait que les valuations (coûts) utilisées sont supposées faire partie d’un en¬ 
semble totalement ordonné. 

Un VCSP est un ensemble de variables prenant des valeurs dans des ensembles 
discrets appelés domaines. Une contrainte valuée est définie à l’aide d’une fonction 
de valuation qui associe un degré de satisfaction à chaque combinaison de valeurs. 
Résoudre un VCSP revient à trouver une affectation de valeurs aux variables avec un 
coût global qui soit optimal. Trouver une telle affectation est une tâche NP-difficile 
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[17]. 

Dans le cadre de cette thèse, nous nous intéressons exclusivement aux VCSPs dits 
binaires. Dans les VCSPs binaires, une contrainte peut porter sur une (contrainte 
unaire) ou deux (contrainte binaire) variables à la fois, mais pas plus. Pour ce qui est 
des méthodes de résolution, nous nous intéressons aux méthodes dites exhaustives, 
c’est-à-dire, des méthodes qui permettent de trouver une des solutions optimales, 
si une telle solution existe. Ces méthodes sont opposées aux méthodes dites appro¬ 
chées ou heuristiques, ou encore par recherche locale. Ces dernières ne garantissent 
pas l’obtention de la solution optimale, mais peuvent s’en approcher en des temps 
relativement courts [118]. 

Les VCSPs à contraintes binaires dures (VCSP-CbC) sont des CSPs classiques 
aux quels on ajoute un ensemble de contraintes unaires flexibles. Motivés par le fait 
que le problème dual 3 d’un VCSP quelconque est un VCSP-CbC [23], nous nous 
sommes attelés à la tâche de proposer des solutions algorithmiques optimisées pour 
la résolution de ces derniers. Notre point de départ était les travaux exposés dans 
[42] où les auteurs ont proposé la notion de substituabilité de voisinage et l’ont utilisé 
comme un moyen de simplification des CSPs binaires. Plusieurs variantes de cette 
notion ont été proposées depuis. Parmi ces variantes celle de substituabilité direc¬ 
tionnelle [91]. Cette dernière est une forme faible de la substituabilité de voisinage 
qui, au départ, visait à améliorer la résolution des CSPs binaires. 

Comme première contribution, nous proposons, dans le Chapitre 4, une méthode 
de décomposition de VCSPs en passant par leurs duaux, c’est-à-dire, les VCSPs 
binaires avec Contraintes binaires Dures (VCSPs-CbC). Plus précisément, nous pro¬ 
posons une extension du travail présenté dans [91] en adaptant la notion de sub¬ 
stituabilité directionnelle aux VCSPs-CbC. La méthode de décomposition proposée 
ramène la résolution d’un VCSPs-CbC à la résolution de plusieurs VCSPs-CbC po¬ 
lynomiaux. 

La résolution d’un VCSP par le biais de son dual présente, toutefois, un incon¬ 
vénient c’est que le nombre de variables, dans le problème dual d’un VCSP binaire, 
est de l’ordre du carré du nombre de variables du VCSP primai. Il est donc intéres¬ 
sant d’envisager de trouver des solutions algorithmiques qui résolvent directement 
les VCSPs sans passer par leurs duaux. 

La complexité algorithmique de la recherche d’une solution optimale d’un VCSP 
a été étudiée dans plusieurs travaux de recherche et plusieurs classes polynomiales 
de VCSPs ont été identifiées et résolues [24, 70, 17]. Par exemple, les VCSPs limités 

3. Le problème dual d’un VCSP, est le problème dans lequel chaque contrainte du problème primai est 
représentée par une variable dont les valeurs possibles correspondent aux combinaison de valeurs admises 
par la contrainte du problème primai (par analogie au dual d’un CSP [48]). 
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à des fonctions de valuations binaires dites sous-modulaires 4 forment une classe de 
problèmes polynomiaux. Les VCSPs binaires dans cette classe peuvent être résolus 
en 0(n 3 d 3 ) étapes où n est le nombre de variables et d est la taille du domaine le plus 
large [14]. Ceci a été possible en exprimant le VCSP binaire comme un problème de 
recherche de coupe minimale dans un graphe orienté et pondéré. 

Notre deuxième contribution, qui est exposée dans le Chapitre 5, consiste en une 
stratégie de décomposition de problème visant les VCSPs binaires qui ne se limitent 
pas à des fonctions de valuations sous-modulaires. Le processus de décomposition 
produit un ensemble de sous-problèmes dont les fonctions de valuations sont toutes 
modulaires 5 et dont la résolution est polynomiale. 

Le choix des fonctions de valuation modulaires dans notre schéma de décompo¬ 
sition est justifié par la simplicité des traitements qu’elles impliquent. Plus précisé¬ 
ment, nous proposons un algorithme qui résout de tels VCSPs en 0(ed) oit e est 
le nombre de contraintes et d est la taille du domaine le plus large. Dans le cas où 
le problème à résoudre ne se limite pas à des fonctions de valuations modulaires, 
nous décrivons comment il est possible de le décomposer en un ensemble de sous- 
problèmes impliquant, exclusivement, des fonctions de valuations modulaires. De 
tels sous-problèmes peuvent être résolus en temps polynomial. 

Enfin, nous décrivons, dans le Chapitre 6, une super-classe de VCSPs binaires 
sous-modulaires. Cette super-classe regroupe les VCSPs binaires dit sous-modulaires 
par permutation. Cette super-classe, est obtenue en appliquant, de manière indé¬ 
pendante, des permutations sur chaque domaine de valeurs. Nous montrons que les 
VCSPs binaires appartenant à cette classe peuvent être identifiés en 0(n 2 d 4 ) étapes, 
et résolus en 0{n 3 d 3 ) étapes où n est le nombre de variables et d est la taille du 
domaine le plus large [14]. 

Nous proposons également un algorithme de décomposition de VCSP dont le 
but est d’améliorer la résolution des VCSPs binaires qui n’appartiennent pas à la 
super-classe proposée. Le processus de décomposition, produit un ensemble de sous- 
problèmes appartenant à la classe polynomiale des VCSPs binaires par permutation. 

Les trois contributions décrites ci-dessus seront détaillées dans les chapitres suivants. 


4. En mathématique une fonction binaire / : R 2 —> R est sous-modulaire si et seulement si pour toutes 
u, v, u', v' G R tel que u < v et v! < v', nous avons f(u, u) + f(v,v') < f(u,v') + f{v,u') 

5. Une fonction binaire / : R 2 —> R est modulaire si et seulement si pour toutes u,v,i i.v' G R, nous 
avons f(u, u') + f(v,v') = f{u,v') + f(v,u') 
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Résumé 

Les Problèmes de Satisfaction de Contraintes Valuées ou VCSPs (Valued Constraints Satisfac¬ 
tion Problems) sont des problèmes où l'on cherche la meilleure combinaison d'objets satisfaisant 
un certain nombre de contraintes valuées. 

Dans ce chapitre nous présentons le Problème de Satisfaction de Contraintes Valuées, une 
classification de ses divers variantes ainsi que les algorithmes classiques qui permettent de le 
résoudre. 


1.1 Introduction 

Nous rencontrons régulièrement des contraintes dans nos vies quotidiennes : un 
espace limité de disque dure dans nos PCs, un nombre de places limité dans le métro, 
un nombre de fréquences limité par support de transmission, une somme d’argent 
limitée dans notre compte bancaire, un nombre d’épouses limité par la loi. Et nous 
nous engageons régulièrement dans la résolution des problèmes de satisfaction et 
optimisation de contraintes : Comment avoir un vaste réseau de télécommunication 
sans avoir d’interférences. Comment vivre en liberté tout en respectant les pensés 
des autres ? Comment réussir un projet sans trop s’endéter ?... En dépit de leur 
simplicités, ces situations donnent lieu à des problèmes combinatoires. La plupart 
du temps, nous n’exigeons pas que des ordinateurs sophistiqués comprennent s’il 
est possible d’écouter la musique à haut volume au milieu de la nuit ou s’il faut 
s’arrêter ail feu rouge. Par contre, quand le nombre de contraintes qui doivent être 
satisfaites et les variables impliquées, commencent à grandir, il devient nécessaire 
de faire recours à une modélisation combinatoire efficace permettant d’optimiser les 
ressources, de réduire les coûts, d’améliorer la qualité et de maximiser les profits. 

Les problèmes de satisfaction de contraintes (CSPs) constituent un cadre simple 
et assez général pour modéliser divers problèmes combinatoires. Pour les CSPs clas¬ 
siques, toutes les solutions sont supposées avoir la même “qualité”. Ce qui fait que la 
résolution d’un CSP se résume souvent à trouver la première solution. En général, un 
CSP n’admet pas une solution unique, et se pose donc les questions suivantes : Ces 
solutions sont-t-clles équivalentes ? Offrent-t-elles réellement la même satisfaction ? 

Désormais, les CSPs classiques ne permettent ni d’exprimer ni de traiter des pré¬ 
férences ou des incertitudes ce qui a créé le besoin de diverses extensions du cadre 
CSP. L’ensemble de ces extensions ont été définies dans [106] comme des Problèmes 
de Satisfaction de Contraintes Flexibles (Soft CSPs). 

Le cadre des CSPs valués (VCSPs) est un peu moins général que celui des Soft 
CSPs dans le sens que les valuations (coûts) utilisées dans les VCSPs représentent 
des niveaux de satisfactions des contraintes et doivent être totalement ordonnées, 
alors que les valuations utilisées dans les Soft CSPs représentent des préférences et 
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peuvent n’être que partiellement ordonnées. 

Les VCSPs font l’objet de recherches intenses à la fois en intelligence artificielle 
et en recherche opérationnelle. La raison est que ces problèmes ont une forte capa¬ 
cité de modélisation vu leur aspect général. En effet, les VCSPs englobent plusieurs 
problèmes plus spécifiques tels que : 

- les CSPs Pondérés (WCSPs), 

- les CSPs Probabilistes (PCSPs), 

- les CSPs Flous (FCSPs), 

- les CSPs Valuées à Contraintes Dures (VCSPs-CC) : 

- les MAX-CSPs, 

- les CSPs classiques. 

Ces extensions des CSPs classiques permettent de modéliser des situations réelles de 
manière plus efficace. Dans le présent chapitre, nous les définissons et nous précisons 
leurs différences. 

Ce chapitre est organisé comme suit : Dans la Section 1.2, nous décrivons les 
Problèmes de Satisfaction de Contraintes Valuées, nous justifions le besoin d’une 
modélisation de plus en plus générale en partant des CSPs classiques vers les VCSPs. 
Dans la Section 1.3 nous présentons une définition précise du VCSP et de ses va¬ 
riantes, et nous étudions la complexité de leur résolution. 


1.2 Problèmes de Satisfaction de Contraintes Valuées (VCSPs) 

1.2.1 Du CSP au VCSP 

Pour justifier le besoin de passer d’une formulation qui considère que toutes 
les contraintes doivent être satisfaites (CSP) vers une formulation qui considère les 
contraintes comme des fonctions de valuations qui représentent des niveaux de sa¬ 
tisfactions (VCSP), nous nous appuyons sur un exemple connu. 


Exemple 1.2.1. (Un problème d’emploi du temps) Imaginons la difficulté de réser¬ 
ver des séances d’enseignement pendant un semestre d’enseignement universitaire. 
Nous devons allouer un créneau horaire, une salle et un professeur pour chaque 
séance d’enseignement. 

En général, les Problèmes de Satisfaction de Contraintes (CSPs) incluent deux 
composantes principales : des variables avec leurs domaines associés et des contraintes. 
Définissons chacune de ces composantes. D’abord, chaque problème de contraintes 
doit inclure des variables. 
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Définition 1. Les variables sont des objets ou des articles qui peuvent prendre une 
variété de valeurs. L ’ensemble de valeurs possibles pour une variable donnée est 
appelé domaine de valeur. 

La deuxième composante d’un problème de satisfaction de contraintes est l’ensemble 
des contraintes. 

Définition 2. Les contraintes sont des règles qui imposent une limitation aux va¬ 
leurs qu ’on peut, simultanément, affecter à une variable, ou à un sous-ensemble de 
variables. 


Exemple 1.2.2. Pour notre problème d’emploi du temps, nous devons aussi satis¬ 
faire plusieurs contraintes. Deux séances d’enseignement ne peuvent pas être tenues 
dans une même salle en même temps. Aucun professeur ne peut assurer deux séances 
d’enseignement différentes en même temps. Une séance d’enseignement ne peut être 
prévue au milieu de la nuit. On doit prévoir à toutes les séances d’enseignement des 
salles de taille convenable ou des amphithéâtres. Les cours de mathématiques doivent 
avoir lieu pendant la matinée. Dans la semaine, les séances de cours doivent avoir 
lieu avant les séances de travaux dirigés correspondantes, etc. 

Dans un CSP, les contraintes sont soit satisfaites soit non satisfaites. Cependant, 
dans plusieurs situations réelles, certaines contraintes peuvent être satisfaites uni¬ 
quement dans la mesure du possible. 


Exemple 1.2.3. Reprenons notre problème d’emploi du temps considéré précédem¬ 
ment. Plus précisément intéressons nous aux deux contraintes qui stipulent que les 
cours de mathématiques doivent être effectués le matin et que, dans la semaine, 
les cours doivent précéder les travaux dirigés associés. Pour ces deux contraintes, 
on peut se trouver dans une situation où, compte tenue des ressources de l’établis¬ 
sement, satisfaire ces contraintes entraine l’élimination de toutes les solutions du 
problème. Alors, si on veut rester dans le cadre des CSPs classiques, on est obligé à 
faire un choix entre ignorer ces contraintes pour espérer avoir une solution qui ne 
les satisfait pas forcement, ou les garder au risque de perdre toutes les solutions. 

Comme pour un CSP, un VCSP peut être représenté par un graphe de contraintes 
qui reflète la structure du problème. 

Définition 3. Un graphe de contraintes d’un VCSP binaire associe, à chaque va¬ 
riable, un sommet et relit chaque paire de sommets reliés par une contrainte. 

Des extensions des CSPs classiques permettent de mieux prendre en compte des 
situations impliquant des contraintes problématiques telles que celles évoquées dans 
l’exemple précédent. Notons que ces extensions nécessitent l’introduction d’une fonc¬ 
tion coût afin de comparer les solutions entre elles. 
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Exemple 1.2.4. Reprenons notre problème d’emploi du temps. 

- Pour la contrainte qui stipule que les cours de mathématique doivent être ef¬ 
fectués le matin, nous pouvons associé un coûts à chaque horaire possible des 
cours de mathématiques. Ce coût doit être croissant en passant de la première 
séance le matin à la dernière de Vaprès midi. Ceci peut se faire à l’aide de 
contraintes unaires flexibles, qui sont des contraintes qui permettent d ’associer 
des coûts à de simples valeurs de variable. Les Problèmes de Satisfaction de 
Contraintes Valuées binaires à Contraintes biaires Dures (VCSPs-CC) 1 uti¬ 
lisent de telles contraintes. 

- Pour la contrainte qui stipule que, dans la semaine, les cours doivent précéder 
les travaux dirigés correspondant, il est possible d’associer un coût à chaque 
travaux dirigés qui précède le cours correspondant. Ceci peut se faire à l ’aide 
de contraintes binaires flexibles. Ces dernières permettent d’associer des coûts 
à des combinaisons de deux valeurs comme il sera précisé dans la section sui¬ 
vante. 


1.2.2 Définition d’un VCSP 

Un Problème de Satisfaction de Contraintes Valuées (VCSP) est alors une exten¬ 
sion du Problème de Satisfaction de Contraintes classique (CSP) en associant des 
coûts aux différentes combinaisons de valeurs que peuvent prendre les variables. Ces 
coûts sont appelés valuations [106]. 


Définition 4. Une valuation est un objet mathématique (habituellement un nombre) 
associé à chaque combinaison de valeurs. Ces valuations expriment l’impact du ni¬ 
veau de satisfaction d’une contrainte sur la “qualité” de la solution. 

Pour définir un VCSP, et en plus des composantes définissant un CSP, on intro¬ 
duit une structure de valuation S qui associe à chaque combinaison de valeurs que 
peuvent simultanément prendre des sous-ensembles de variables, une valuation. 

Définition 5. Une structure de valuation est un triplet S = (E, ©, ^) où : 

- E est l’ensemble de valuations contenant un élément minimal _L et un élément 
maximal T ; 

- P est un ordre total sur E ; 

- © est un opérateur d’agrégation qui est commutatif, associatif et monotone. 

Nous présentons ci-dessous trois remarques importantes à propos de la structure 
de valuation S : 


1. Ces problèmes seront définis dans la sous-section suivante. 
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- Une structure de valuation S est dite propre si pour chaque paire de valuations 
a, (3 G E, tel que a ^ (3, il existe une différence maximale entre fie ta dénotée 
par /3 © a. 

- Un opérateur d’agrégation © est strictement monotone si pour toute a, (3, 7 
dans E tel que a -< (3 et 7 ^ T, nous avons a © 7 -< (3 © 7 . 

- Une structure de valuation S est strictement monotone si elle possède un opé¬ 
rateur d’agrégation strictement monotone. 

Le Problème de Satisfaction de Contraintes Valuées (VCSP), étudié dans cette 
thèse, est alors défini comme suit : 

Définition 6. Un Problème de Satisfaction de contraintes Valuées (VCSP) est dé¬ 
fini par un quadruplet (X, V , C, S) où : 

- X est un ensemble fini de variables ; 

- V est un ensemble fini de domaines, D x e V étant le domaine de x G X ; 

- S = (U,©, X) est une structure de valuation propre et strictement monotone. 

- C est un ensemble de contraintes valuées. Chaque contraintes est un couple 
(a, fi) où 

- a, la portée de la contrainte, est une liste ordonnée d’éléments de X et 

- c t> : U xea D x —y E est une fonction de valuation. 

Pour des raisons de simplicité, nous supposons que les domaines sont des sous 
ensembles d’entiers. Une variable x ne peut être affectée qu’à une valeur de son 
domaine. L ’arité d’une contrainte valuée est la taille de sa portée. L’arité du pro¬ 
blème est l’arité maximale parmi toutes ses contraintes. Dans ce travail, nous nous 
limitons aux VCSPs binaires. Ce sont des VCSPs dont les contraintes sont unaires 
ou binaires. Dans un VCSP binaire, deux variables sont dites voisines si elles sont 
reliées par une contrainte binaire. 

Nous supposons aussi que les fonctions de valuation sont évaluées en temps 
constant. 

Dans le cas des VCSPs binaires, l’ensemble des contraintes peut être exprimé 
comme l’union d’un ensemble de contraintes unaires dénoté par C 1 et un ensemble 
de contraintes binaires dénoté par C 2 . Nous avons aussi besoin du sous-ensemble des 
contraintes, dont les portées sont ordonnées, qui est défini comme suit : 

C = C 1 U{((x,y),fi)eC 2 \x^y} 

où —>■ désigne un quelconque ordre total sur les variables. Le problème (X, V, C , S) 
est alors appelé le VCSP ordonné selon —> associé à (À", V, C , S). Le fait de ne consi¬ 
dérer que le sous-ensemble de contraintes C empêche, comme nous le verons ci-après, 
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les valuations binaires d’être comptées deux fois. De plus l’utilisation de C au lieu de 
C est une restriction qui va intervenir dans, pratiquement, toutes les contributions 
proposées dans cette thèse. 

A présent, nous précisons comment les valuations données par les diverses fonc¬ 
tions de valuations sont combinées et comparées entre elles en vue d’identifier les 
solutions optimales d’un VCSP. 

Une distanciation est l’ensemble des affectations d’une valeur à chaque variable. 
La valuation d’une distanciation t d’un sous-ensemble de variables V Ç X est donnée 
par 


= ® (p c (tia c ) ( 1 . 1 ) 

cGC,ct(c)ÇV 

où t l U dénote la projection de t sur les variables de U Ç V et C est définie par 
rapport à un ordre total quelconque sur X. Une solution optimale pour un VCSP 
V est un n-tuple t* tel que n = |X| et é-p(t*) est minimal selon l’ordre A. 

1.2.3 Exemple 
Les mariages stables 

Nous définissons dans ce paragraphe le problème test connu sous le nom du problème 
des mariages stables [117]. 

Chikh Salah est père de quatre hiles à l’age de mariage. Pour chacune d’elle, il 
a reçu deux demandes en mariage. Chikh Salah souhaite proposer à ses quatre hiles 
des mariages “stables” en tenant compte des préférences de ses hiles, mais aussi de 
l’entente entre les époux pour garantir une belle ambiance dans la famille. Pour ce 
faire, il consulte une agence matrimoniale aux méthodes modernes. 

Pour tenir compte au mieux des désirs de Chikh Salah, l’agence a demandé 
à chacune des hiles de donner un degré de préférence pour chaque homme qui a 
demandé sa main. Désignons chaque hile par P) : i : 1,2,3,4 et chaque éventuel 
époux par Hj , j : 1,2,..., 8. Les degrés de satisfaction des éventuels époux sont 
donnés dans le Tableau 1.1. 

En coopération avec le père et en étudiant les dossiers des éventuels époux, 
l’agence a réussit à exprimer les ententes entre les éventuels époux par des valuations 
qui sont données dans le Tableau 1.2. 

Ainsi, l’agence a modélisé le désir de Chikh Salah comme un Problème de Sa¬ 
tisfaction de Contraintes Valuées binaire. Selon cette modélisation, les variables 
représentent les éventuels époux de chaque hile (époux de P)). Pour chaque va¬ 
riable, l’agence a associé le domaine Di composé des hommes qui ont demandé 
la main de la hile P). L’ensemble des contraintes unaires C 1 n’est autre que les 
degrés de satisfaction des éventuels époux. Tandis que l’ensemble des contraintes 
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éventuel époux 

satisfaction 


1 

h 2 

2 

h 3 

2 

h a 

1 

h 5 

2 

h 6 

1 

h 7 

2 

h 8 

1 


Tableau 1.1- Les degrés de satisfaction des éventuels époux 


Entente 

Hs 

H 4 

h 7 

h 8 


1 

1 

1 

2 

h 2 

3 

2 

3 

1 

h 5 

3 

2 

3 

2 

h 6 

1 

1 

1 

2 


Tableau 1.2 - Les valuations de l’entente entre les éventuels époux 


Epoux de F2 



Figure 1.1 - Une solution qui maximise les désirs du Chikh Salah. (Les valeurs qui 
constituent la solution sont indiquées par des ellipses à fond noir) 
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binaires C 2 est l’éventuelle entente entre Hj et Hji avec 1 < j, j' < 8 ,j 7^ j'. 
Ainsi la solution la plus stable est celle qui est associée à une valuation globale 
maximale. L’une de ces solutions (Voir Figure 1.1) est celle formée par les couples 
(Fi, H 2 ), (F 2 , H 3 ), (F 3 , H 5 ), (F 4 , H 7 ) et dont la valuation est de 20. 

1.3 Classification des VCSPs 

Dans cette section, nous présentons une classification des VCSPs et leurs va¬ 
riantes. Cette classification est résumée dans l’arborescence de la Figure 1.2. 


VCSP 



FCSP 

Sf =([0 1],MAX,<) 


WCSP 

Sw(k) fc},©,r<) 


PCSP 

Sp =([0 1],x,<) 



FCSP-CC MAX-CSP 

Sf=({0,1},MAX,<) Sm =(IN,+,<) 


WCSP-CC PCSP-CC 

Sw =({0,“},+,<) Sp=({0,1},x,<) 



Figure 1.2 - Arborescence des VCSPs 


1.3.1 Problèmes de Satisfaction de Contraintes Floues (FCSP) 

Les Problèmes de Satisfaction de Contraintes Floues (Fuzzy Constraint Satisfaction 
Problems) FCSPs [101] prolongent la notion des CSPs classiques en permettant des 
contraintes floues, c’est-à-dire des contraintes définies à l’aide de relations floues. 
Les contraintes floues utilisent des valuations comprises entre 0 et 1. Une valuation 
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Ûij 

1 

2 

1 

0 

0 

2 

1 (si i = j [ 2 ]) 

0.75 


0 (sinon) 



Tableau 1.3 - La fonction de valuation 

de 1 représente la meilleure valuation (c’est-à-dire que la combinaison est permise 
avec certitude). Une valuation nulle indique que la combinaison n’est pas permise. 
L’ensemble des valuations E est donc l’intervalle [0 1]. Dans cet ensemble de valua¬ 
tion, 0 est l’élément neutre et 1 l’élément absorbant. L’ordre dans cet ensemble de 
valuations est l’ordre total classique sur les réels (<). L’opérateur d’agrégation n’est 
autre que l’opérateur max qui renvoi le maximum de ces arguments et qui, en plus 
d’être monotone, commutatif, associatif, est aussi idempotent 2 . 

La structure de valuations permettant de définir les FCSPs est donc la suivante : 

Définition 7. La structure de valuation floue est le triplet 
Sf = ([0 1 ], max, <) 

La solution d’un CSP fkm est alors définie comme le n-uplet t* qui a une valua¬ 
tion maximale. 


Exemple 1.3.1. Soit V un VCSP composé de trois variables (entreprises) e 1: e 2 
et e 3 . (Voir Figure 1.3 gauche). Les domaines de chaque variable e; sont formés 
de deux valeurs (quantité à produire) Di = { 1 , 2 }. Ce VCSP est composé de trois 
contraintes Ci = ((ei, e 2 ), </>i 2 ), c 2 = ((e 2 , e 3 ), </> 2 3 ) et c 3 = ((ei, e 3 ), <j> 13 ) dont la 
fonction de valuations faj est définie par le Tableau 1.3. 

La solution de ce problème en utilisant la structure de valuation Sp est repré¬ 
sentée dans la Figure 1.3 droite. 

1.3.2 CSP Probabilistes (PCSPs) 

Dans les Problèmes de Satisfaction de Contraintes Probabilistes PCSPs [105], chaque 
contrainte c a une certaine probabilité p(c), indépendante des probabilités des autres 
contraintes. En réalité, la probabilité n’est pas propre à la contrainte, mais à la si¬ 
tuation qui correspond à la contrainte. Dire que c a la probabilité p signifie que la 
situation correspondante à c a une probabilité p d’apparition dans le problème réel. 
Evidement, p(c) est toujours entre 0 et 1. L’ensemble de valuations E est donc l’in¬ 
tervalle [0 1], où 0 est l’élément neutre et 1 est l’élément absorbant. L’ordre dans cet 
ensemble de valuations est l’ordre classique sur les réels (<) et l’opérateur d’agréga¬ 
tion n’est autre que la multiplication (x). 


2. En mathématique, la propriété d’idempotence signifie qu’une opération, appliquée une ou plusieurs 
fois, a le même effet (l’opération est réappliquée sans modification de résultat). L’opérateur max vérifie la 
propriété d’idempotence. On dit qu’il est idempotent. 
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Figure 1.3 - VCSP et sa solution par la structure de valuation Sf- (Les valeurs qui 
constituent la solution sont indiquées par des cercles à fond vert) 


Toutes ces données se résument dans la structure de valuation probabiliste suivante : 
Définition 8. La structure de valuation probabiliste est le triplet 

S P = ([0 1], x,<) 

Exemple 1.3.2. Reprenons le même VCSP défini dans l’exemple précédent (Voir 
Figure 1.3 gauche). La solution de ce problème en utilisant la structure de valuation 
Sp est représentée dans la Figure 1.4. On peut constater que ce n’est pas la même 
solution qu’avec la structure floue. 



Figure 1.4 - Maximisation de profits des entreprises du VCSP de l’exemple précédent 
par la structure de valuation Sp. (Les valeurs qui constituent la solution sont indiquées 
par des cercles à fond vert) 
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1.3.3 CSP Pondérées (WCSPs) 

Les problèmes de satisfaction de contraintes pondérées, (WCSPs) [75], associent 
un coût à certaines combinaisons de valeurs. Ce coût est toujours compris dans 
l’ensemble de valuations ordonnées {1,..., k}, où k est une constante prédéterminée. 
L’opérateur d’agrégation © est défini par ctffi/3 = min{/c, a + b} et son inverse partiel 
© est défini pour a ■< (3 par /3Qa = ksi/3 = ket/3 — a sinon. Notons que Sw(k ) 
n’est pas strictement monotone pour k -< T. L’ordre dans cet ensemble de valuations 
est dénoté W La structure de valuations associée au CSP pondérées est donc définie 
comme suit : 

Définition 9. Une structure de valuation pondérée est un triplet 

Sw{k) = ({1, • • •, k}, ffi, <). 


Ces WCSPs présentent un modèle pour les problèmes d’optimisation où le but est 
de réduire au minimum le coût total de la solution recherchée. Ce coût peut être un 
temps, une distance, un nombre de ressources... 

Ainsi, pour les WCSPs, le coût d’une solution est obtenu en sommant les pondé¬ 
rations associées aux combinaisons de valeurs utilisées. Et le but est de trouver un 
n-uplets t* ayant un coût minimal. 

On peut voir les FCSPs et les WCSPs comme deux approches différentes per¬ 
mettant de donner une signification à la notion d’optimisation. Les deux modèles 
correspondent à deux définitions de protection sociale dans la théorie utilitaire [83] : 

- L’égalitarisme 3 , qui optimise l’utilité individuelle minimale pour les FCSPs. 

- L’utilitarisme 4 , qui optimise la somme des utilités individuelles pour les WCSPs. 

Les FCSPs sont basés sur l’approche qui vise à maximiser le niveau de cohérence 
en équilibrant les niveaux de toutes les contraintes, tandis que les WCSPs ont une 
approche qui vise à obtenir le coût minimal à l’échelle globale, bien que quelques 
contraintes puissent être négligées même si elles présentent un coût important. 

1.3.4 CSPs classiques 

Un CSP se définit par la donnée d’un ensemble de variables, chaque variable pouvant 
prendre une valeur parmi un ensemble de valeurs possibles appelé domaine. Puis, par 
un ensemble de contraintes exprimant des restrictions sur les combinaisons de valeurs 
permises. Dans un Problème de Satisfaction de Contraintes, les contraintes doivent 
être obligatoirement satisfaites. Donc, toutes les contraintes sont dures. Ainsi, l’en¬ 
semble des valuations E est la paire {0, +oo}, où 0 est l’élément neutre et +oo est 

3. En matière politique, économique ou sociale, l’égalitarisme peut être défini comme l’égalité des ci¬ 
toyens. L’égalitarisme a pour objectif politique suprême l’égalité matérielle de tous. 

4. L’utilitarisme prescrit d’agir (ou ne pas agir) tout en optimisant le bien-être global de l’ensemble des 
êtres sensibles. 
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Homme 

ne désire pas 

Ft\ 

F 3 et F 4 

h 2 

Pi 

h 3 

f 2 

h 4 

Fi et F 3 


Tableau 1.4 - La liste de conflits des hommes 


Femme 

ne désire pas 

Fi 

H 2 , H 3 et H 4 

f 2 

h 3 

Fs 

H\. H 2 et H \ 

Fs 

Hi 


Tableau 1.5 - La liste de conflits des femmes 

l’élément absorbant. L’ordre dans cet ensemble de valuations est l’ordre donné par 
0 < Too. L’opérateur d’agrégation n’est autre que l’addition (+). Ainsi on a : 

Définition 10. La structure de valuations pour un CSP est le triplet 
Scsp = ({0, +oo}, +, <). 

La valuation d’une affectation t d’un ensemble de variables V Ç. X est donnée 
par (1.1). Une solution optimale pour un CSP V est un n-tuple t tel que n = |X| et 
= 0 . 


Les mariages réalisables 

Il s’agit d’une version simplifiée des mariages optimaux réalisables. Cette fois, notre 
agence matrimoniale moderne souhaite proposer à ses clients des mariages “réali¬ 
sables”. Elle demande pour cela à chacun de ses membres de déterminer une liste 
des membres du sexe opposé qui sont non désirés, établissant de la sorte une liste 
de conflits. Il s’agit alors, en tenant compte de ces listes de conflits de former des 
couples possibles. Par possible, on entend que personne ne doit épouser une personne 
de sa liste de conflit. 

Supposons que l’agence matrimoniale a toujours 4 hommes : Hi, H 2 , H 3 et H 4 , 
et 4 femmes : F \, F 2 , F 3 et F 4 . Leurs conflits exprimés sont donnés dans les Tableaux 
1.4 et 1.5. 

Une modélisation possible de ce problème consiste à représenter chaque client, 
homme ou femme, par une variable : Xi,... ,x%, tandis que les domaines des variables 
sont : pour les époux D h = {1,... ,4} et pour les épouses Df = {5,... ,8 }, et la 
contrainte sur chaque paire de variables est que Xf ^ Xh si / G Df ne désire pas 
h G D h ou si h ne désire pas /. 
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Epouse de H2 


Epouse de Hi 



Figure 1.5 - Une solution qui satisfait toutes les contraintes de l’exemple de mariages 
réalisables. (Les valeurs qui constituent la solution sont indiquées par des ellipses à fond 
gris) 


En Tunisie, chaque homme ne peut pas avoir plus d’une épouse, de même pour les 
femmes. D’où on ne peut pas avoir une même épouse (époux) pour deux hommes 
(femmes) différents (différentes) donc, on doit avoir Xh ^ Xh> et Xf ^ Xf 
pour tout /, /' G Df et / ^ f et pour tout h, h' G Df t , et h ^ h'. 

Néanmoins, sachant que l’on peut retrouver le mari d’une femme en cherchant 
l’homme dont elle est l’épouse, 011 peut ne garder comme variables que les épouses. 
La modélisation de ce problème sera de représenter chaque épouse de Hf par une 
variable, d’où les quatre variables xi,... ,x±. Les domaines de ces variables sont les 
épouses possibles Dh = {Fi,..., F 4 }. 

Une solution réalisable (voir Figure 1.5) consisterait, par exemple, à former les 
couples suivants (//i,iq), (H 2 ,F 2 ), ( H 3 , F 3 ), (H 4 , F A ). 

1.3.5 CSP Valuées à Contraintes binaire Dures (VCSPs-CbC) 

Rappelons qu’un VCSP à Contraintes binaires Dures [29] est un CSP classique 
augmenté d’un ensemble de contraintes unaires flexibles. 

Donc, en plus des composantes qui définissent le problème de satisfaction de 
contrainte sous-jacent, un VCSP-CbC fait intervenir une structure de valuation qui 
permet de définir des fonctions de coût unaire permettant de mesurer la qualité de 
la meilleure solution. 

Dans un VCSP-CbC, les contraintes binaires sont toutes des contraintes dures. 
Ainsi, l’ensemble des valuations pour ces contraintes est la paire {0, + 00 }, Par contre 
les contraintes unaires sont flexibles. Un coût est associé à chaque valeur. Ainsi, l’en- 
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semble des valuations pour ces contraintes unaires est M + . L’opérateur d’agrégation 
n’est autre que l’addition (+). Ainsi la structure de valuation d’un VCSP-CbC est : 

Définition 11. La structure de valuations pour un VCSP-CbC est définie par : 
S C c = (Ei U E 2 , +, <) tel que : 

- Ei = E est l’ensemble des valuations flexibles des contraintes unaires ; 

- E ‘2 — {0,+oo} est l’ensemble des valuations des contraintes binaires dures; 

Notre intérêt pour les VCSPs-CbC provient du fait que le problème dual d’un 
VCSP (même non binaire) est un VCSP-CbC. 

Les mariages optimaux réalisables 

Reprenons l’exemple des mariages stables et définissons les mariages optimaux réali¬ 
sables. Supposons que dans l’agence matrimoniale, il y a 4 hommes (que l’on appel¬ 
lera H i, H 2 , H 3 et H fi) et 4 femmes (que l’on appellera Fi, F 2 , F 3 et Ffi). Les clients 
de notre agence sont les hommes. Après avoir demander à chacun des hommes de 
déterminer une liste des membres du sexe opposé non désirés, établissant de la sorte 
une liste de conflits. Notre agence matrimoniale souhaite dans ce cas de proposer à 
ses clients (les hommes) des mariages “optimaux réalisables”. Elle doit alors proposer 
pour chaque homme une épouse. Elle détermine pour cela pour chacune des femmes 
un score selon sa beauté, son intelligence, son statut social ... aux yeux de chaque 
client. 

Il s’agit alors, à partir des listes de conflits et en tenant compte des scores de 
chaque femme, de former quatre couples optimaux réalisables. Par couples optimaux 
réalisables, on entend qu’aucun homme ne doit être soumis à la tentation de divorcer 
pour trouver mieux ailleurs : si Antar est marié à Laila et Abla n’a pas de marie, 
et si Antar préfère Abla à Laila alors ce mariage ne sera pas optimal réalisable car 
Antar sera tenté de divorcer pour trouver sa Abla. Ou encore, si Antar est marié à 
Laila et Kais est marié à Abla et si Antar préfère Abla à Laila et Kais préfère Laila 
à Abla, alors ces mariages ne seront pas optimaux réalisables car Antar sera tenté 
de divorcer pour trouver sa Abla et Kais sera tenté de divorcer pour trouver sa Laila. 

Les conflits et les degrés de préférences sont les suivants : 


Hommes 

ne désire pas 

Hi 

F 3 et Fi 

h 2 

Fi et F 3 

H-i 

Fi et F 2 

Hi 

Fi et F 3 


Tableau 1.6 - Les conflits des hommes 
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Score de Satisfaction pour 

Hi 

H 2 

h 3 

H a 

Fi 

2 

0 

0 

0 

f 2 

1 

1 

0 

2 

f 3 

0 

0 

2 

0 

Fa 

0 

2 

1 

0 


Tableau 1.7 - Les degrés de satisfaction des femmes 


Epouse de H2 


Epouse de Hi 



Figure 1.6 - Une solution qui satisfait toutes les contraintes et qui maximise la satisfaction 
des clients de l’exemple de mariages optimaux réalisables. (Les valeurs qui constituent la 
solution sont indiquées par des ellipses à fond gris) 
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Au vu de ces préférences, le mariages de H 2 avec F 2 n’est pas optimal réalisable 
car H 2 préférera F 4 à F 2 à cause des scores donc, H 2 va tenter de divorcer pour 
épouser F 4 . En revanche, une solution optimale réalisable consisterait (voir Figure 
1.6), à former les couples (H 2 ,F 4 ), (. H 3 ,F 3 ) et (H 4 ,F 2 ). 

La satisfaction de cette solution est maximale et elle a comme coût : 8. 

1.4 Conclusion 

Dans ce chapitre nous avons montré le besoin d’enrichir le modèle standard des pro¬ 
blèmes de satisfaction de contrainte en vue de tenir compte des subtilités que l’on 
peut rencontrer lors de la modélisation de situations réelles. Nous avons, ensuite, 
présenté une classification des différentes variantes du Problème de Satisfaction de 
Contraintes Valuées, un formalisme offrant la possibilité de modélisation des pro¬ 
blèmes combinatoires impliquant la satisfaction et l’optimisation de contraintes. 
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Résumé 

Dans ce chapitre nous présentons les principales techniques de résolution des problèmes de 
satisfaction de contraintes valuées. L'idée fondamentale sur laquelle s’appuient ces techniques 
est la réduction de l’espace de recherche à explorer. Ces techniques tirent profit du coût de la 
meilleure solution trouvée pour qu’un algorithme du type “recherche en profondeur d'abord” explore 
un espace de recherche de taille réduite. Ces techniques permettent de guider la recherche vers la 
solution optimale, accélérant ainsi la résolution des problèmes. 


2.1 Introduction 

En théorie, pour résoudre des VCSPs dont les domaines sont finis, il est possible, 
afin de trouver une solution minimisant la fonction objective, d’énumérer toutes 
les possibilités représentant l’espace de recherche à explorer. Pour ce faire, on fait 
généralement recours à un algorithme de recherche arborescente s’appuyant sur le 
mécanisme du retour arrière chronologique. Cependant, vu le grand nombre de com¬ 
binaisons possibles, on fait souvent recours à des techniques, dites de propagation de 
contraintes (ou filtrage). Ces techniques permettent de réduire l’espace de recherche 
en éliminant de chaque domaine les valeurs qui ne peuvent pas participer à une so¬ 
lution réalisable, et encore moins, à une solution optimale. Pour ce faire, on intègre 
un algorithme de filtrage à celui assurant la recherche arborescente. 

Bien que les algorithmes de résolution des VCSPs diffèrent de ceux dédiés aux 
CSPs classiques, la plupart en sont des adaptations. Cependant, alors que dans le 
cadre CSP classique on se limite, en général, à trouver une seule solution, pour ce 
qui est des VCSPs, la tâche est beaucoup plus difficile car il s’agit de trouver la so¬ 
lution de meilleure valuation. Donc au lieu d’utiliser des algorithmes de type retour 
arrière simple, qui seraient trop coûteux en temps de calcul, il est plus intelligent 
de profiter du coût de la meilleure solution trouvée (appelée borne supérieure) pour 
épargner l’exploration de branches inutiles. Des branches dont le coût dépasse cette 
borne supérieure. Un tel algorithme est connu sous le nom de l’algorithme de branch 
and bound. L’algorithme de filtrage intégré peut profiter de la borne supérieure per¬ 
mettant ainsi de guider, plus rapidement et plus efficacement, la recherche vers la 
solution optimale. 

Ce chapitre est organisé comme suit : dans la Section 2.2, nous décrivons les tech¬ 
niques de filtrages utilisées pour résoudre les problèmes de satisfaction de contraintes 
classiques (CSPs). Ensuite, dans la Section 2.3, nous présentons les principales amé¬ 
liorations aux algorithmes de filtrage et d’exploration de l’arbre de recherche afin 
de les rendre plus efficaces lors de la résolution des VCSPs. Nous discutons les 
heuristiques de choix de variables et de valeurs dans la Section 2.4. Une solution 
algorithmique existante pour la résolution des VCSPs sera présentée dans la Section 
2.5. 
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2.2 Propagation de contraintes 

Afin de pouvoir résoudre un CSP dans un temps raisonnable, on fait souvent recours 
à une combinaison d’algorithmes. Pour trouver la solution d’un CSP, il faut trouver la 
solution permettant de satisfaire toutes les contraintes. L’exploration d’une partie 
de l’espace de recherche peut être évitée seulement si on peut prouver qu’aucune 
solution n’appartient à cette partie. Les parties infructueuses sont identifiées par 
des algorithmes spécifiques, dit de propagation de contraintes, plus connus sous le 
nom générique d’algorithmes de filtrages. Ces algorithmes s’appuient sur la notion 
de consistance locale. 

La propagation de contraintes joue un rôle important dans la résolution efficace 
des problèmes de satisfaction de contraintes [82]. La propagation de contraintes est 
définie par le fait de réduire le domaine d’une variable en éliminant l’ensemble des 
valeurs infructueuses : celles qui sont incohérentes avec les contraintes du problème. 
Dans le cas où la propagation permet d’établir la cohérence globale du problème 
alors elle permet de dériver une solution à ce problème. O 11 peut distinguer plusieurs 
niveaux de consistances locales, plus ou moins fortes : 

- Consistance de nœud (1-consistance). 

- Consistance d’arcs (2-consistance). 

- Consistance de chemin (3-consistance) 

- /s-consistance. 

2.2.1 Consistance de nœud 

Ce niveau consiste à ne considérer qu’une seule variable à la fois. On retire alors 
toutes les valeurs pour lesquelles, au moins, une contrainte est impossible à satisfaire. 
Ce filtrage est très rapide, mais peu efficace car les suppressions des valeurs ne 
peuvent avoir lieu que s’il y a des contraintes unaires dans le problème. Ce qui n’est 
pas toujours le cas. 

2.2.2 Consistance d’arcs 

La consistance d’arcs est définie comme suit : 

- On dit qu’une contrainte binaire c xy est arc-consistante (AC) si et seulement si 
pour toute valeur appartenant au domaine de la variable x, il existe au moins 
une valeur dans le domaine de la variable y qui satisfait c xy . 

- On dit qu’un problème est AC si toutes ses contraintes sont AC. 

Il existe plusieurs algorithmes permettant d’établir l’arc-consistance d’un CSP. 
Dans ce qui suit, nous présentons un algorithme d’arc-consistance AC3 qui est très 
utilisé et qui est simple à implémenter [82], 
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L’algorithme AC3 

Le principe général de AC3 est le suivant : 

Soit la contrainte binaire ({x,y) ,(j> xy ), u parcourant toutes les valeurs possibles 
de la variable x et v parcourant toutes les valeurs possibles de la variable y. Si, pour 
un u donné, (p(u,v) = +oo pour toutes les valeurs v de la variable y , alors il faut 
enlever u du domaine de la variable x. Si une valeur u a été enlevée du domaine 
de x, alors il faut répéter ces vérifications pour toutes les contraintes impliquant la 
variable x. Cet algorithme nécessite donc de gérer une pile contenant les variables 
dont le domaine est modifié. On peut montrer que si d est la taille maximale des 
domaines des variables et e le nombre de contraintes, alors : 

- Une contrainte binaire est examinée au plus d fois (car une contrainte est exa¬ 
minée à cause de la suppression d’une valeur dans le domaine de sa première 
variable, et qu’on ne peut pas supprimer plus de d valeurs du domaine). Il y a 
donc au plus ed examens de contraintes. 

- Chaque examen de contraintes nécessite d 2 vérifications car il y a d 2 couples 
dans le produit cartésien D x x D y . L’algorithme a donc une complexité en 
0 (ed 3 ). 


Algorithme 1: AC3(U) 

// avec P = (X, D, C, S) 

PlLEG- 4> 

pour chaque ({x,y) ,(/>) G C faire 
Pile g- Pile u{(x,y)} 

tant que Pile / 0 faire 

SÉLECTIONNER^, y) 

SUPPRIMER^, y) (de la Pile) 

Réviser^, e- faux 
pour chaque u G D x faire 

si il n’existe pas de support ( u , v) alors 
supprimer(u) 

RÉVISERD x VRAI 
si RÉviserZU alors 

Pile g- Pile U{(c, x)\((z, x) , (f>) G C} A ((z, y) , <f>) G C 

RENVOYER V 
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2.2.3 /c-consistance 

La /c-consistance consiste à considérer k variables, et à tester tous les fc-uplets 
de valeurs possibles afin de vérifier s’ils satisfont aux contraintes. Plus k est grand, 
plus le filtrage est efficace. Cependant, du fait du grand nombre de combinaisons à 
tester le fitrage devient vite trop lourd pour k dépassant 4. La 3—consistance donne 
de bons résultats, cependant du fait de la complexité de son implémentation, elle 
n’est que très rarement présente dans des solveurs de contraintes. La 2—consistance 
correspond en fait à la consistance d’arc. C’est le niveau de consistance le plus utilisé 
en pratique. 

En intégrant une technique de propagation de contraintes à un algorithme de 
recherche arborescente, on obtient l’algorithme de résolution connu sous le nom de 
Maintient dArc Consistance (MAC) [103]. 

2.3 Résolution d’un VCSP 

Afin de pouvoir résoudre un VCSP dans un temps raisonnable, nous faisons souvent 
recours à une combinaison d’algorithmes. Pour trouver la solution optimale, il faut 
trouver la solution permettant de satisfaire au mieux les contraintes c’est-à-dire celle 
la moins coûteuse. Une partie de l’espace de recherche peut être élaguée seulement 
si on peut prouver que la solution optimale n’appartient pas à cette partie. Les 
parties infructueuses sont identifiées par des algorithmes spécifiques connus sous le 
nom générique d’algorithmes de filtrages. Ces algorithmes s’appuient sur la notion 
de consistance locale flexible. 

2.3.1 Consistance locale flexible (Soft local consistency) 

La consistance locale flexible, que nous étudions ci-dessous, joue un rôle impor¬ 
tant dans 1a, résolution efficace des VCSPs [73, 75, 76]. Calculer une borne inférieure 
du coût de 1a, solution optimale permet d’éviter l’exploration de branches inutiles et 
d’accélérer la, convergence vers la solution optimale. 

Le calcul de cette borne ainsi que l’élimination des branches inutiles sont assurés 
par des algorithmes de consistance locale flexibles. 

WAC3 

L’algorithme WAC3 [75] est un algorithme qui applique l’arc-consistance flexible 
sur les WCSPs, et ceci en s’appuiant sur le principe de fonctionnement de AC3 
[82], L’algorithme utilise deux sous-programmes : la fonction SupprimerVar (voir 
Algorithme 2) qui supprime les valeurs inconsistantes dans D x et la procédure 
TrouverSupportsAC3 (voir Algorithme 3) qui projette la, fonction de valuation 
( t> xy sur la fonction de valuation unaire <f> x . 

Tout d’abord, commençons par définir la notion de support simple. 
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Définition 12. Une valeur b G D y est un support simple pour a G D x si 4> xy (a, b) = 
_L. 


Algorithme 2: SupprimerVar(x) 

CHANGER G- FAUX 

pour chaque a G D x A 4> x ( a ) = T faire 
Dx G- D x \ {a} 

CHANGER G- VRAI 

RENVOYER CHANGER 


Algorithme 3: TrouverSupportsAC3(x, y) 

pour chaque a G D x faire 

î oc i min 6eDy {4> xy {a,b)} 

2 4> x {a) G- </> x (a) © a 

pour chaque b G faire 

3 |_ 4>xy(a,b ) G- ( i>xy(a,b)Qa 


Pour chaque valeur a de D x , WAC3 détermine une valeur v E D y tel que (j) xy (a , u) 
est minimale. Cette valuation minimale sera stockée dans une variable de travail a 
(voir Algorithme 3, ligne 1). Si la valeur v est complètement compatible avec a 
(c’est-à-dire, (p xy (a,v ) = _L) les lignes 2 et 3 n’auront aucun effet. Sinon, la valeur 
v devient un support simple pour a en envoyant le coût a de la fonction binaire 
(j) xy (a,v) vers la fonction unaire (p x (a) (lignes 2 et 3). 


Algorithme 4: WAC3 (X,D,C) 


Pile g- {1,2,... ,rt} 
tant que PlLE / 0 faire 
y g- pop(Pile) 
pour chaque cj} xy G 4> faire 

TrouverSupportsAC3(x, y) 
si SupprimerVar(x) alors 
Pile g- Pile U {x} 


La pile utilisée dans WAC3 contient les variables dont les domaines ont été ré¬ 
duits, et par conséquent, les variables adjacentes peuvent avoir de nouvelles valeurs 
qui n’ont pas de support simple. La Pile est initialisée avec toutes les variables 
(Algorithme 4, ligne 1 ), car chaque variable doit avoir un support simple initial 
pour chacune de ses valeurs et pour chacune des contraintes où elle est impliquée. 
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La boucle principale est répétée tant que la Pile n’est pas vide. La variable y est ex¬ 
traite de la tête de la pile (ligne 2) et pour toute variable x de la contrainte ((x, y) ,<j>), 
de nouveaux supports sont trouvés dans D x , si nécessaire (ligne 3). Comme trouver 
de nouveaux supports dans D x peut augmenter les coûts de 0 X , la consistance de 
nœud dans D x doit être vérifiée et les valeurs incohérentes doivent être supprimées 
(ligne 4). Si D x est modifié, la variable x est ajoutée à la pile, car les variables 
connectées à x doivent avoir leurs supports révisés. Si pendant le processus certains 
domaines deviennent vides, l’algorithme peut être interrompu avec la certitude que 
le problème ne peut pas être résolu avec un coût inférieur à T. 


Théorème 1. (Larosa et al. [75]) La complexité temporelle de WAC3 est de l’ordre 
de 0(ed 3 ), où e et d sont respectivement le nombre de contraintes et la plus large 
taille des domaines. 

2.3.2 Étoile arc-consistance (AC*) 

L’algorithme arc-consistance étoile (AC*) [75] est un algorithme de filtrage dédié 
aux WCSPs. 

Notons par 0 0 la borne inférieure du problème définie par : 

0 0 = V min 4> x (v) 

xCX VeD * 

Cette équation suppose qu’une contrainte unaire pèse sur chaque variable du pro¬ 
blème. Ceci ne constitue pas une limitation, car on peut toujours associer une 
contrainte unaire à chaque variable sans changer le problème. 

Le principe de fonctionnement de AC* est le suivant : 

(i) 00 doit être mis à jour après les projections de contraintes binaires sur les 
contraintes unaires 1 * . 

(ii) chaque fois que 0 0 est mis à jour, tous les domaines doivent être revérifiés afin 
d’identifier d’éventuelles valeurs inconsistantes. 

Définition 13. Une variable x est étoile nœud consistante si pour toute valeur 
a G D x , 00 ©0x(o) -< T, et il existe une valeur a G D x tel que 0 æ (a) = _L. Un VCSP 
est nœud consistant (NC*J si toutes ses variables sont étoile nœud consistantes. 

1. Une projection d’une contrainte binaire c xy sur deux contraintes unaires c x et c y doit vérifier que : 

- Toute solution optimale de c xy est une solution optimale de c x et de c y 

- Toute solution optimale de c x peut être étendue à une solution optimale de c xy 

- Toute solution optimale de c y peut être étendue à une solution optimale de c xy . 
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Avant l’exécution de AC*, le problème doit être NC* 2 . Un problème AC* est 
défini comme suit : 

Définition 14. Une variable x est étoile arc-consistante si toute valeur a G D x a un 
support simple 3 dans ép xy . Un VCSP est arc-consistant (AC*) si toutes ses variables 
sont étoile nœud consistantes et étoile arc-consistantes. 

La démarche de AC* est similaire à celle de WAC3. Nous discutons seulement 
les principales différences. La Fonction SUPPRIMERVar diffère par le fait que (p 0 
est utilisée en vue de supprimer des valeurs. La Fonction TrouverSupportsAC*(x, y) 
projette (p xy sur (j) x (lignes 2-7) et par la suite (ft x est projeté sur 0 0 (lignes 8-10) 
pour restaurer le support pour la consistance de nœud de x. Elle renvoie VRAI si 0 0 
est augmentée. Dans la boucle principale, lorsque SupprimerSupportsAC*(x, y) 
renvoie VRAI, chaque variable doit être vérifiée pour l’inconsistance de valeurs (lignes 
1). La raison est que les augmentations de (p 0 peuvent causer des inconsistances de 
nœud dans n’importe quel domaine. 

Théorème 2. (Larosa et al. [15]) La complexité temporelle de AC* est 0{n 2 d 2 + 
ed 3 ), où n,e et d sont le nombre de variables, contraintes et la plus large taille des 
domaines. 


Algorithme 5: TrouverSupportsAC*(x, y) 


1 

2 

3 

4 

5 


TROUVER G- FAUX 

pour chaque a € D x faire 

a G- min bGDy {(j) X y(a, &)} 

v g- Nil 

si <j) x (a) = lAa!^l alors 
TROUVER V- VRAI 


6 <j> x {a) G- <j> x (a) © a 

pour chaque b G D y faire 

7 |_ 4>xy(a,b ) V- 4> xy (a, b) Q a 


s ol i min aeDx {<^(a)} 


pour chaque a G D x faire 

îo |_ 4> x (a) G- 4> x (a) © a 
RENVOYER TROUVER 


AC* a été conçu pour les WCSPs binaires. Rappelons qu’un WCSP utilise la 
structure de valuations Sw(k) = ({ 1 ,..., k}, ©, <), où k est une constante prédéter¬ 
minée. L’opérateur d’agrégation © est défini de telle sorte que affi/3 = min{fc,a + f>} 

2. En général, pour une affectation complète, la combinaison minimale des coûts relatifs aux contraintes 
unaires représente une borne inférieure. Cette idées intégrée à la définition de consistance de nœuds produit 
la consistance de nœuds étoile dénotée par NC*. 

3. Rappelons que b G D y est un support simple pour a G D x si <j> xy (a, b) = _L. (Voir Définition 12) 
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Algorithme 6: TrouverSupportsCompletsAC*(x, y) 


1 

2 

3 

4 

5 


TROUVER G- FAUX 
pour chaque a G D x faire 

P (a) G- mm b&Dy {(j) X y(a, b) © <j) x {a)} 
v g- Nil 

si <f) x (a) =1A P(a) >- _L alors 

TROUVER -U- VRAI 


pour chaque b G D y faire 

6 |_ E(a) G- max aeDx {(j) X y(a, b) 0 P (a)} 


7 

8 


pour chaque b G D y faire 

Mb) <- ^(6) 0Ê(i) 

pour chaque a G faire 


( b) 


M( a ’ b ) © E{b) 


pour chaque a G D x faire 

9 (j) x (a) G- M a ) © P ( a ) 

pour chaque b G faire 

io |_ 4> xy ( a ’b) t- 4> xy ( a , b) © P(a) 

n ce G- min aGDx {çi a; (a)} 

12 00 •<— 00 © O! 

pour chaque a G D x faire 

13 |_ 0 æ (a) G- 0a;(a) © a 

RENVOYER TROUVER 


et son inverse partiel © est défini pour a © /3 par /3©a = /csi/3 = /cet/3 — a sinon. 
Notons que Sw{k ) n’est pas strictement monotone pour k -< T. Ceci est nécessaire 
pour permettre le filtrage de valeurs ayant une valuation finie. 

Les algorithmes de filtrage sont, en général, incrémentaux. En d’autres termes, 
si une consistance est établie dans un nœud de l’arbre de recherche alors l’établisse¬ 
ment de la consistance locale dans un nœud fils peut se faire en ne considérant que 
les modifications entre celui-ci et le nœud père. Cette propriété est très utile car le 
filtrage est effectué à chaque nœud de l’arbre de recherche. 

En pratique, AC* est intégré à un algorithme du type retour sur trace chrono¬ 
logique et utilisé pour résoudre les VCSPs, qui sont plus généraux que les WCSPs. 
Toutefois, pendent cette intégration le paramètre k doit varier avec le coût de la 
meilleure solution déjà trouvée. Ce point est essentiel pour que ces algorithmes 
puissent effectuer des filtrages qui tirent profit du coût de la meilleure solution déjà 
trouvée. 
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Algorithme 7: TrouverSupportsExtentielAC*(:e, y) 


1 

2 

3 

4 


TROUVER G- FAUX 

a G- min a 6 23 a! {(^ x (a) ©0 x1/e ^ | min 6eDi/ {^ 2/ (6) © <j) xy (a,b)}} 

si a >- _L alors 

pour chaque (f) xy G cj) \ y — >• x faire 

L 


TROUVER G- TROUVER V TROUVERSUPPORTSCOMPLETSAC*(x, y) 


RENVOYER TROUVER 


Algorithme 8: SupprimerVar*(x) 

CHANGER G- FAUX 

pour chaque a G D x \ 4> x (a) © (f >0 = T faire 
D x G- D x \ {a} 

CHANGER G- VRAI 

RENVOYER CHANGER 


2.3.3 Améliorations de AC* 

En se référant à [76], l’algorithme de filtrage de WCSP binaire DAC* (pour Direct 
AC) est une première amélioration de AC*. Pour décrire le fonctionnement de cet 
algorithme, on a besoin d’introduire quelques nouvelles notions. 

Définition 15. Soit ({x,y) ,4> x , y ) une contrainte binaire flexible. Une valeur b G D y 
est un support complet pour a G D x si <j) xy (a , b) © <j> y (b) = _L. 

Définition 16. Une variable x est étoile complètement arc-consistante si toute 
valeur a G D x a un support complet dans le domaine de chaque variable voisine, 
et il existe a G D x tel que <f> x (a) = _L. Un VCSP est étoile complètement arc- 
consistant si toutes ses variables sont étoile nœud consistantes et étoile complètement 
arc-consistantes. 

DAC* consiste à combiner les coûts binaires et les coûts unaires dans le calcul 
de la valuation minimale à projeter. Pour ce faire, DAC* fait appel à un nouvel 
algorithme : Algorithme 6 au lieu de l’Algorithme 5. Cette nouvelle opération de 
base consiste à établir des supports complets pour les valeurs de la variable x pour 
chaque contrainte ({x,y) ,(j> xy ). Ceci peut se faire par l’extension du coût unaire de 
4>y vers (p xy et la projection de (f) xy sur (p x . Cependant, l’extension de tous les coûts 
unaires peut détruire les supports de y sur (f) xy . 

En se référant à [76], FDAC* est une amélioration de AC* et DAC*. Cet algo¬ 
rithme suppose un nombre total sur les variables noté —X 

Définition 17. Une variable x est arc-consistante directionnelle 
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Algorithme 9: AC *(X,D,C) 


PileQ -g- X 

tant que PileQ A 0 faire 
y <r- pop(PileQ) 
pour chaque 4> xy G <j) faire 

si TrouverSupportsAC*(x, y) alors 
pileR g- pileR U {x} 

pour chaque x G X A SuPPRlMERVÀR*(x) faire 
1 pileQ g- pileQ U {x} 


- si toute valeur a G D x a un support complet dans le domaine de chaque variable 
voisine y tel que x —> y. 

- il existe une valeur a G D x tel que 4> x (a) = A. 

Elle est arc-consistante directionnelle complet (FDAC) 

- si, en plus, chaque valeur a G D x a un support simple dans chaque contrainte 
((x, y) , 4 >) tel que y — > x. 

Un VCSP est arc-consistant directionnel complet (FDAC*,) 

- si toutes ses variables sont nœud consistantes et arc-consistantes direction¬ 
nelles complets. 

En effet, l’intégration simultanée de AC* et DAC* permet d’établir FDAC* 
[76], les supports complets pour la variable x sur (f> xy sont obtenus en préservant les 
supports pour toutes les valeurs de y sur (j > xy . Ceci est obtenu par l’extension du 
coût minimal de (f) y requis pour la projection ultérieure sur (j) x . 

L’Algorithme 5 permet de trouver des supports complets pour toutes les valeurs 
de x sur <j) xy tout en prenant soin des supports pour les valeurs de D y . Cet algorithme 
retourne VRAI lorsque le coût de la valeur a G D x a été augmenté de _L. L’Algorithme 
10 a été conçu pour être utilisé tout seul pour établir DAC* ou conjointement avec 
l’Algorithme 9 pour établir FDAC*. Par conséquent, chaque fois qu’une valeur est 
supprimées, DAC* insère cette variable dans pileQ pour informer AC* de la sup¬ 
pression. DAC* utilise PILER : une file de priorité qui contient les variables telles 
que le coût unaire a été augmenté par rapport à _L. Dans ce cas, plusieurs valeurs 
dans des variables d’ordre antérieur peuvent avoir perdu le support complet et de 
nouveaux supports doivent être trouvées. La boucle principale itère tant que pileR 
n’est pas vide. A chaque itération, la dernière variable y est récupérée à partir de 
pileR. Les valeurs incohérentes (en raison du coût unaire et l’incrémentation de 
la borne inférieure) sont éliminées en utilisant SUPPRIMERVAR* (Algorithme 8) et 
les variables supprimées seront insérées dans pileQ. Puis des nouveaux supports 
complets seront recherchés pour chaque variable inférieure à y. 

Enfin, toutes les variables seront traitées pour établir NC* qui peut être perdu 
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au cours du processus, en raison de l’incrémentation de la borne inférieure. Les va¬ 
riables supprimées seront insérées dans pileQ. FDAC* établit simultanément AC* 
et DAC* : rétablissement de AC* videra pileQ mais il pourra ajouter des variables 
à pileR, et l’établissement de DAC* videra pileR mais il pourra ajouter des va¬ 
riables à pileQ. FDAC* terminera lorsque pileR et pileQ seront simultanément 
vides. 

La complexité temporelle de DAC* est de O(ed 2 ) et la complexité spatiale est de 
O (ed). Quand à FDAC* la complexité temporelle est de O (end 3 ) et la complexité 
spatiale est de O (ed). 


Algorithme 10: DAC *(X,D,C) 
pileR g- X 

tant que PILER / 0 faire 
y g- pop(pileR) 

pour chaque (j) xy G (j> \ x —> y faire 

si TrouverSupportsCompletsAC*(x, y) alors 
pileR g- pileR U {x} 
pileS g- pileS U {x} 

pour chaque x G A A SuPPRlMERVÀR*(x) faire 
î pileQ g- pileQ U {x} 


Algorithme 11: EAC *(X,D,C) 
pileR g- X 

pileP g- {z\x G S,z > x, <j> xz G 4>} U S 
S ^0 

tant que PILEP A 0 faire 
x g- pop(pileP) 

si TrouverSupportsExtentielAC*(x) alors 
pileR g- pileR U {x} 

pour chaque cf> xy G (j)\x —> y faire 
î pileP g- pileP U {y} 


Les supports complets peuvent être établis en deux directions si leur établissement 
produit une augmentation dans la borne inférieure. Il s’agit d’une propriété de co¬ 
hérence locale naturelle, appelée arc-consistance existentielle (EAC* [76]). 


Définition 18. Une variable x est arc-consistante existentielle s’il y a au moins 
une valeur a G D x tel que (p x (a) — _L et elle a un support complet dans le domaine 
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Algorithme 12: FDAC *(X,D,C) 

pileR g- A 
pileQ g- X 

tant que PILeQ V 0VpileR V 0 faire 
AC *{X,D,C) 

DAC *(X,D,C) 


de chaque variable voisine. Un VCSP est arc-consistant existentiel (eac*) si toutes 
ses variables sont nœud consistantes et arc-consistantes existentielles. Un VCSP est 
edac* s’il est fdac* et EAC*. 

Si une variable x n’est pas eac*, alors pour toute valeur a tel que <j> x (a) = _L, 
il y a une variable y tel que pour toute valeur b G D y , (p xy (a, b) © (j> y (b) > _L. Ceci 
nécessite de rétablir les supports complets de x et par conséquent de briser la pro¬ 
priété de NC*. Il sera donc possible d’augmenter la borne inférieure <p 0 . 

En se référant à [76], il est possible d’intégrer eac* avec fdac* afin d’exploiter 
les bénéfices de chacun. Ce qui nous donne la consistance d’arc directionnelle exis¬ 
tentielle (edac*). 

EDAC* exige que chaque valeur a un support complet dans une direction et un 
support simple dans l’autre direction (pour satisfaire FDAC*). En plus, au moins une 
valeur par variable doit avoir un support complet dans les deux directions (pour 
satisfaire eac*). Cette valeur spéciale est appelée la valeur support complet. 


Algorithme 13: EDAC* (A, D, C) 


pileR g- X 
pileQ g- X 
pileS a- X 

tant que pileS ^ 0VpileQ ^ 0VpileR / 0 faire 
EAC *(A,D,C) 

DAC *{X,D,C) 
tant que pileQ V 0 faire 
y G- pop(pileQ) 
pour chaque <f) xy G (f>\y —> x faire 

si TrouverSupportsAC*(x, y) alors 
pileR g- pileR U {x} 
pileS g- pileS U {x} 

pour chaque x G A A SupprimerVar*(x) faire 
pileQ g- pileQ U {x} 


En pratique, les algorithmes AC*, DAC*, FDAC* et EDAC* sont intégrés à un al- 
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gorithme du type retour sur trace chronologique et utilisés pour résoudre les VCSPs. 

2.4 Heuristiques de sélection de variables et de valeurs 

Les heuristiques de sélection des variables et des valeurs jouent un rôle fondamental 
dans la résolution des VCSPs. En effet, une bonne sélection permet de définir une 
stratégies de recherche efficace. 

Dans [60] les auteurs ont étudié différents ordres sur les variables, nous présentons 
ci-dessous les deux types les plus importants : 

1. En amant, avant de débuter l’exploration de l’espace de recherche, des ordres, 
dits statiques, sont calculés et maintenus fixes pendant la résolution. De tels 
ordres ont l’avantage d’éviter des calculs supplémentaires pendant la résolu¬ 
tion comme ils n’interfèrent pas avec l’algorithme de filtrage. Ces algorithmes 
exploitent des informations collectées lors de l’étude du problème original, on 
ne prend pas compte des modifications faites au cours de la recherche par l’af¬ 
fectation de variables. La limite de tels ordres et la perte de l’exactitude de 
l’information à la suite des hltages qui ont lieu pendant la résolution. 

2 . A l’inverse des ordres statiques, des ordres dynamiques sont calculés à chaque 
nœud de l’arbre de recherche en utilisant les informations du sous-problème 
courant. De tels ordre ont l’avantage d’être plus précis. La limite de tels ordres 
et le besoin de plus de temps de calcul. 

Dans [60] les auteurs ont montré que la sélection, statique, de variable la plus 
robuste est celle qui est basée sur l’ordre inverse du degré des variables c’est-à-dire 
le nombre de sommets adjacents à la variable dans un graphe de contraintes 4 . 

Pour les heuristiques dynamiques d’ordre des variables, celle qui semble la plus 
efficace est celle qui s’appuie sur la somme des coûts unaires. Dans [27] les auteurs 
ont obtenu de bonnes performances pour résoudre des WCSPs aléatoires en amé¬ 
liorant une sélection dynamique d’ordre. Ils ont appelé cette heuristique, “2-sided 
Jeroslow-like”, les variables sont ordonnées selon un ratio entre la taille de leur do¬ 
maine et la somme des moyennes des coûts unaires et binaires des contraintes qui 
sont liées à chaque variable. 

Pour l’heuristique de sélection de valeur dans un problème de minimisation de 
coût, la valeur instanciée en premier sera souvent celle qui a le coût unaire le plus 
faible ou la combinaison des coûts la plus faible. 

4. Un graphe de contraintes d’un CSP binaire associe, à chaque variable, un sommet et relit chaque 
couple de sommets { x , y} tel que C xy € C. 
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2.5 La solution algorithmique 

Les opérations élémentaires effectuées au cours du filtrage par arc-consistance di¬ 
rectionnelle existanciclle sur les VCSPs binaires sont la suppression de valeurs, la 
projection et l’extension de coûts. Toutes ces opérations ne peuvent ajouter des 
contraintes binaires au problème sur lequel la technique de filtrage EDAC* est ap¬ 
pliquée. Donc le problème demeure binaire. Une autre propriété cruciale est que la 
technique de filtrage EDAC* calcule une borne inférieure pour le coût de la solution 
optimale et quand le VCSP est unaire, cette borne inférieure correspond exactement 
au coût de la solution optimale. 

L’algorithme le plus utilisé pour la résolution des VCSPs est l’algorithme qui 
maintient l’arc-consistance directionnelle existanciclle medac*. Ce dernier résulte 
de l’intégration de l’algorithme de filtrage EDAC* à un algorithme de recherche ar¬ 
borescente s’appuyant sur le mécanisme du retour sur trace chronologique. Nous 
présentons ci-dessous deux versions de l’algorithme medac* : 

- une version qui explore un arbre de recherche binaire, désigné par 2wmedac* 
(Voir Algorithme 14) ; 

- une version qui explore un arbre de recherche fc-aire, désigné par /cwmedac* 
(Voir Algorithme 15). 

L’ordre des variables est dynamique, c’est-à-dire qu’il peut varier pendant la 
recherche. Cet ordre est implicitement déterminé par la liste des variables non en¬ 
core instanciées Y, alors que les variables de X — Y précèdent ceux de Y dans l’ordre. 

A chaque nœud de l’arbre de recherche 2 wmedac* traite un sous-problème dé¬ 
noté par V , si le nœud traité est une feuilles, ce qui signifie que toutes les variables 
du problème initiale ont été traitées, le coût de la solution courante sera évalué et 
comparé au meilleur dernier coût trouvé pendant l’appel récursif ultérieur (ligne 1). 

Le plus faible entre les deux sera stocké dans c*. Sinon, si le nœud traité est un nœud 
interne, 2wmedac* sélectionne une variable x non encore affectée. Il sélectionne, en¬ 
suite, une valeur v G D x et réduit le domaine de x à v selon une heuristique de choix 
de valeur (ligne 4). Ensuite, l’algorithme EDAC* de filtrage est exécuté. Pour que l’al¬ 
gorithme puisse passer à un niveau plus profond dans l’arbre de recherche, le VCSP 
résultant ne doit pas contenir de domaine vide. Si l’exécution de l’algorithme de fil¬ 
trage ne conduit pas à un domaine vide, un appel récursif à 2 wmedac* est effectué 
dans le but de traiter les variables restantes. Dans une étape ultérieure, l’algorithme 
ignore l’effet de réduire D x à (x , v) et écarte l’élément (x , v) de D x . Il restaure encore 
une fois l’arc-consistance directionnelle existantiellc du sous problème résultant et 
teste à nouveau s’il n’y a aucun domaine vide. Si c’est le cas, l’algorithme exécute 
un deuxième appel récursif. 

La principale différence avec la version 2wmedac* est que avec la version fcWMEDAC*, 
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Algorithme 14: 2 wmedac*(R = (X,V,C,S),Y) : c* 

si Y = 0 alors 

si <p'(t) -< c* alors 

1 |_ c* <- 

sinon 

X t— SELECT(y) 

2 (x,v) u- Select (D x ) 

3 D x <- {(æ,v)} 

EDAC* (P) 

si 0 £ R alors 

c* t- 2 wmedac*(P, Y \ {x}) 

RESTAURE (R) 

D x <r~ D x \ {{x,v)} 

EDAC* (R) 
si 0 ^ T> alors 

C* U- 2WMEDAC *(V,Y) 

RESTAURE (R) 


le domaine de la variable courante est réduit, à tour de rôle, à chaque valeur du do¬ 
maine de cette variable, afin de préserver la complétude de l’algorithme. Ceci est 
assuré par la boucle tant que de la ligne 2.11 en résulte que l’arbre de recherche qui 
est explorée par /cwmedac* est un arbre fc-aire. 

2.6 Conclusion 

Dans ce chapitre nous avons présenté les principales techniques existantes de résolu¬ 
tion des Problèmes de Satisfaction de Contraintes Valuées binaires. L’idée fondamen¬ 
tale sur laquelle s’appuient ces techniques, est la réduction de l’espace de recherche 
à explorer. Les algorithmes dits de filtrage tirent profit du coût de la meilleure solu¬ 
tion déjà trouvée pour réduire l’espace à explorer par l’algorithme de résolution (un 
algorithme de retour sur trace chronologique, en général). Ainsi, le filtrage permet de 
guider, plus rapidement et plus efficacement, la recherche vers 1a, solution optimale. 
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Algorithme 15: fcWMEDAC*^, Y, c*) : c* 

// avec V = (X,V,C,S) 
si Y = 0 alors 

si 4>'(t) -< c* alors 
1 |_ c* <- (f>'(t) 


2 

3 

4 


sinon 

X <r- SELECT(Y) 

tant que true faire 

(x,v) <r- Select (D x ) 

D x <- {(x,x)} 

EDAC* (T) 

si 0 ^ V alors 

C* 4- fcWMEDAC* (P, Y \ {x}, C*) 


restaure(P) 

D x D x \ {(x,x)} 

EDAC*(P) 

si 0 £ V alors 

BREAK 


restaure(P) 
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Résumé 

Un Problème de Satisfaction de Contraintes Valuées (VCSP) est un problème d'optimisation NP- 
Difficile. Une partie importante de la recherche sur les VCSPs consiste à identifier des restrictions 
donnant lieu à des sous-problèmes solubles en temps polynomial. Ces sous-problèmes forment des 
classes polynomiales de VCSP. L’identification et l’utilisation de classes polynomiales spécifiques 
de VCSP a fait l'objet de plusieurs travaux de recherches dans ces dernières années. 

Dans ce chapitre, nous présentons trois grande classes de VCSPs polynomiaux. Nous utiliserons 
ces trois classes dans les chapitres suivants pour résoudre plus efficacement des VCSPs plus 
généraux. 


3.1 Introduction 

Une partie importante de la recherche sur les VCSPs consiste à identifier des 
restrictions sur les différentes composantes définissant les VCSPs qui garantissent 
une résolution en temps polynomial. On distingue deux types de restrictions : des 
restrictions qui touchent la structure des instances et des restrictions qui limitent la 
famille des fonctions de valuation autorisées dans le problème. Ces dernières sont des 
restrictions dites de langage. Plusieurs restrictions, limitant le langage des fonctions 
de valuation et garantissant la polynomialité, ont été identifiées dans la littérature. 
Nous citons principalement les travaux décris dans [16, 24]. 

Ce chapitre est structuré comme suit, dans la Section 3.2, nous étudions la com¬ 
plexité des VCSPs. Puis, dans la Section 3.3, nous décrivons la famille des fonctions 
de valuation sous-modulaires. Cette dernière est considérée comme la famille de res¬ 
triction de langage la plus importante pour résoudre les VCSPs. Ensuite, dans la 
Section 3.3.9, nous présentons la famille des fonctions modulaires. Enfin, dans la 
Section 3.4, nous présentons les VCSPs polynomiaux par substituabilité. 

3.2 Complexité d’un VCSP 

Pour pouvoir déterminer 1a, complexité d’un VCSP, nous présentons quelques 
notions de base sur la théorie de la complexité. Nous commençons par définir deux 
types de machines abstraites de calcul : les machines déterministes et les machines 
non déterministes. 

Définition 19. Une machine abstraite de calcul est dite déterministe si chaque 
action possible est unique, en d’autres termes l’action à effectuer est dictée de façon 
unique par l'état courant. Dans le cas où il peut y avoir plusieurs choix possibles 
d'actions à effectuer, la machine est dite non déterministe. 

A l’aide de ces deux types de machines, on définit la classe des problèmes poly¬ 
nomiaux. 
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Définition 20. Un problème de décision est dans la classe Polynomiale s ’il peut être 
décidé sur une machine déterministe en temps polynomial par rapport à la taille de 
la donnée du problème. Ce problème est alors qualifié de polynomial. Ce qui veut 
dire concrètement qu ’il existe donc un algorithme polynomial pour le résoudre. 

Nous définissons maintenant une classe, plus générale, de problèmes, à savoir la 
classe des problèmes Non-déterministe Polynomiaux (NP). 

Définition 21. Un problème de décision est dans une classe NP s ’il peut être décidé 
sur une machine non déterministe en temps polynomial par rapport à la taille de sa 
donnée. 

On qualifie alors le problème de Non-déterministe Polynomial. En d’autres ternie, 
un problème dans la classe NP est un problème qui peut être résolu par énuméra¬ 
tion de l’ensemble des solutions possibles dont chacune peut être testée à l’aide d’un 
algorithme polynomial. Il existe alors un algorithme polynomial non déterministe 
pour ce problème. 

Nous définissons à présent un problème plus difficile, à savoir un problème NP- 
complet : 

Définition 22. Un problème est NP-complet s’il vérifie les propriétés suivantes : 

1. Chaque solution peut être testée à l’aide d’un algorithme polynomial. 

2. Tous les problèmes de la classe NP peuvent se ramener à ce problème par une 
transformation polynomiale. En d’autres termes, un problème dans la classe 
NP-complet est au moins aussi difficile que tous les autres problèmes de la 
classe NP. 

Nous définissons, finalement, un problème encore plus difficile, à savoir un pro¬ 
blème NP-difficilc : 

Définition 23. Un problème est NP-difficile si tous les problèmes de la classe NP 
peuvent se ramener à ce problème par une transformation polynomiale. Notons que, 
chaque solution ne peut pas être testée à l’aide d’un algorithme polynomial déter¬ 
ministe. Il n’existe donc pas d ’algorithme polynomial non déterministe pour ce pro¬ 
blème. 

Enfin, nous pouvons déterminer la complexité d’un VCSP. 

Théorème 3. Un Problème de Satisfaction de Contraintes Valuées est un problème 
NP-difficile [82]. 

La Figure 3.1 illustre la localisation de la classe des VCSPs par rapport aux diffé¬ 
rentes classes de complexité. En conclusion, résoudre un VCSP est une tâche difficile 
du point de vue de la complexité puisqu’il s’agit de résoudre un problème NP-difficilc. 
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Figure 3.1 - Étude de complexité d’un VCSP : Diagramme d’Euler. 


3.3 Fonctions de valuations sous-modulaires 

Les fonctions sous-modulaires sont au centre de l’optimisation combinatoire et de 
l’intelligence artificielle. Dans cette thèse, nous nous intéressons à bénéficier des pro¬ 
priétés des fonctions de coût sous-modulaires afin d’accélérer le processus de résolu¬ 
tion des VCSPs. Dans un premier paragraphe, nous définissons la sous-modularité. 
Ensuite, nous présentons des exemples de fonctions sous-modulaires. Puis, nous pré¬ 
sentons un bref historique de la minimisation des fonctions sous-modulaires. A la fin 
de cette section, nous décrivons comment il est possible d’utiliser les fonctions de 
valuations sous-modulaires pour la résolution des VCSPs. 

3.3.1 Qu’est-ce que la sous-modularité? 

Dans [85] l’auteur a présenté une étude détaillée qui décrit la minimisation des 
fonctions sous-modulaires. Supposons qu’une usine est capable de produire tout 
sous-ensemble de E, où E est un ensemble donné de produits potentiels. Si l’usine 
décide de produire un sous-ensemble S Ç E de produits, alors elle doit amortir un 
coût d’installation c(S) afin que l’usine soit prête à produire S. Ce coût d’installation 
est un cas particulier d’une fonction de réglage : 

Etant donné un ensemble fini E (le motif prévu), la notation 2 e représente la 
famille de tous les sous-ensembles de E. Nous définissons, alors, une fonction scalaire 
/ : 2 e —y M. qui est appelée une fonction d’ensemble. Notons par f(S) les valeurs 
de / sur le sous ensemble S Ç E, nous utilisons n pour désigner \E |. 

Supposons que l’usine a décidé provisoirement de produire le sous-ensemble S, 
et que l’usine envisage la possibilité d’ajouter des produits e ^ S à sa gamme de 
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produits. Dans ce cas, le coût d’installation supplémentaire que l’usine aurait à payer 
est c(S U {e}) — c(S). Afin de simplifier les notions, nous utilisons la notation stan¬ 
dard S + e qui signifie S U {e}, S — e qui signifie S — {e}, et /(e) qui signifie /({e}). 
Dans cette notation le coût marginal d’ajouter e est c(S + e) — c(S). Nous utilisons 
S C T qui signifie S Ç T mais S ^ T. 

Notons que, les économistes suggèrent que, dans la plupart des situations du 
monde réel, ce coût supplémentaire est une fonction décroissante de S. Autrement 
dit, ajouter un nouveau produit e à un ensemble plus vaste devrait produire un 
coût marginal moins important que d’ajouter e à un ensemble plus petit. Pour une 
fonction / générale nous devrions avoir 

VScTcT + e, f(S + e) — f (S) > f (T + e) — f (T) (3.1) 

Quand une fonction / quelconque satisfait 3.1, alors on dit qu’elle est sous- 
modulaire. La connexion entre la sous-modularité et l’économie est très profonde. 
Beaucoup plus de détails à ce sujet sont disponibles dans le livre de D. M. Topkis 

[114]- 

3.3.2 Exemples de fonctions sous-modulaires 

Vu leur importance, les fonctions sous-modulaires ont été largement utilisées [85]. On 
peut citer les exemples des fonctions de rang des matroïdes [100], ou en théorie des 
graphes [41], la fonction qui associe à tout sous-ensemble de sommets d’un graphe 
la cardinalité de sa coupe [31]. 

D’autres exemples de fonctions sous-modulaires plus connues sont : 

- toutes les fonctions unaires sur un ensemble totalement ordonné, 

- les fonctions linéaires J2i a i x ;, 

- la fonction de distance Euclidienne ^ ( x 2 — x 1 ) 2 + ( y 2 — y 1 ) 2 . 

3.3.3 Historique de la minimisation des fonctions sous-modulaires (SFM) 

La minimisation des fonctions sous-modulaires (SFM) a été reconnue comme un 
problème important depuis la naissance de l’optimisation combinatoire, quand au 
début des années 1970 Edmonds [32] a établi deux résultats fondamentaux (voir 
Tableau 3.1) : 

- Optimisation gloutonne sur les Polyèdres sous-modulaires. 

- Les outils algorithmiques pour les Polyèdres sous-modulaires. 

Lorsque l’algorithme de l’Ellipsoïde a vu le jour en 1981, Grsotschcl, Lov’asz, et 
Schrijver [52] se sont rendus compte que c’est un outil fondamental pour trouver 
des algorithmes polynomiaux pour des problèmes tels que la SFM. En mathéma¬ 
tiques d’optimisation, la méthode de l’ellipsoïde est une méthode itérative utilisée 
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Edmonds [32] 

(1970) 

Optimisation gloutonne sur les Polyèdres sous-modulaires. 

Les outils algorithmiques pour les Polyèdres sous-modulaires. 

Grsotschel, Lov’asz, et 
Schrijver [52] (1981) 

Création d’un outil fondamental : 

L’algorithme polynomial de l’Ellipsoïde. 

Critiques 

L’algorithme d’Ellipsoîde n’est pas très pratique pour 
l’optimisation combinatoire. 


Tableau 3.1 - Travaux clés sur les SFMs 


pour minimiser des fonctions convexes. La méthode de l’ellipsoïde, conçue pour la 
résolution de problèmes d’optimisation linéaire avec des données rationnelles, est un 
algorithme qui trouve une solution optimale en un nombre fini d’itérations. Afin de 
converger vers un minimum d’une fonction convexe, la méthode de l’ellipsoïde génère 
une séquence d’ellipsoïdes dont le volume diminue uniformément à chaque itération. 

Toutefois, ce résultat est finalement peu satisfaisant, puisque l’algorithme Ellip¬ 
soïde n’est pas très pratique, pour l’optimisation combinatoire. La question est donc 
passée de “Est-ce que le problème SFM est polynomial?” (voir Tableau 3.1) à “Y 
a-t-il un algorithme combinatoire polynomial, non-Ellipsoïde, pour la SFM ?” (voir 
Tableau 3.2). 


Cunningham [7, 25, 26] 
(1984-1985) 

Présenter un algorithme pseudo-polynomial pour la SFM. 

Schriiver [1091 
(2000) 

Proposer un algorithme combinatoire : 

Algorithme de Schrijver. 

Iwata, Fleischer, et 
Fujishige [64] (2001) 

Proposer un autre algorithme combinatoire : 

Algorithme IFF. 

Orlin [97] 

(2007) 

Proposer un nouveau algorithme combinatoire : 
nettement plus rapide. 

Critiques 

Tous ces algorithmes utilisent des techniques d’algèbre linéaire. 

Iwata [65] 

(2009) 

Proposer un algorithme entièrement combinatoire hybride : 
basé sur l’algorithme IFF et l’algorithme de Schrijver à la fois. 

Orlin [98] 

(2009) 

Proposer un nouveau algorithme entièrement combinatoire : 
le plus rapide, (de complexité 0(n 6 + n 5 Ef )) 

Observations 

La plus part des algorithmes SFM utilisent la notion de flot 
maximal : l’intersection de matroïdes [33]. 


Tableau 3.2 - Algorithmes combinatoires (non-Ellipsoide) pour la SFM. 


En 1985, Cunningham [26] a déclaré que : “Trouver un algorithme pratique po¬ 
lynomial et combinatoire pour minimiser une fonction générale sous-modulaire est 
un problème ouvert”. Cunningham a donné ce qui s’est avéré être les principales 
contributions à ce problème dans le milieu des années 80, en utilisant un résultat de 
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programmation linéaire 1 de la dualité d’Edmonds 2 [32], il a proposé une solution 
algorithmique qui utilise l’algorithme du flot maximal pour la SFM. Les travaux de 
Cunningham présentent un algorithme pseudo-polynomial pour la SFM [7, 25, 26]. 
Puis, presque simultanément en 1999, deux travaux semblaient donner différents al¬ 
gorithmes combinatoires fortement polynomiaux 3 pour la SFM : 

- Schrijver [109] (publié en 2000). 

- Iwata, Fleischer, et Fujishige (IFF) [64] (publié en 2001). 

- En 2007, Orlin [97] a proposé un algorithme différent pour la SFM qui est 
nettement plus rapide. 

Notons que ces trois algorithmes sont basés sur les travaux de Cunningham. En 
plus, tous ces algorithmes utilisent un sous-progrannne nommé “l’algorithme Cara- 
théodory”. L’entrée de l’algorithme “Carathéodory” est une représentation d’un vec¬ 
teur y EW 1 comme une combinaison convexe de sommets c’est à dire y = AjîA 
La sortie est un ensemble d’au plus n des v l tels que l’enveloppe convexe de ces 
vecteurs contient toujours y. Ceci s’appuie sur des techniques d’algèbre linéaire. Ce 
qui a conduit Schrijver [109] à poser la question de savoir s’il existe un algorithme 
entièrement combinatoire 4 pour la SFM (voir tableau 3.2) ? 

Une version de cet algorithme entièrement combinatoire a été donnée par Iwata 
[65]. Puis améliorée par Orlin [98]. 

On peut remarquer que la plus part des algorithmes combinatoires SFM utilisent 
la notion de flot, dans la résolution des problèmes sous-modulaires. Ces algorithmes 
utilisent une forme généralisée de minimisation de flot en utilisant l’intersection de 
matroïdes. Cette dernière est développée par Edmonds et Giles [33]. Ces algorithmes 
sont dits algorithmes de flots sous-modulaires. 

1. La programmation linéaire est une branche de l’optimisation. Elle permet de résoudre de nombreux 
problèmes économiques et industriels. Un problème d’optimisation linéaire est un problème d’optimisation 
dans lequel on minimise une fonction linéaire sur un polyèdre convexe. La fonction-coût et les contraintes 
peuvent donc être décrites par des fonctions linéaires (on devrait dire affines), d’où vient le nom donné à 
ces problèmes. 

2. La dualisation est une technique. Elle est utilisée pour introduire à un problème d’optimisation 
initiale un problème dual équivalent. Cette technique prouve son utilité lorsque le problème dual introduit 
peut être plus facile à résoudre. 

3. Un algorithme est dit fortement polynomial s’il est polynomial et que de plus, lorsque les opération 
arithmétiques sont comptées pour 1 dans le temps de calcul mesuré pour des données numiriques, celui-ci 
ne dépend plus de la longueur d’écriture de ces données [84]. 

4. Un algorithme est dit entièrement combinatoire lorsque les configurations de collections d’objets 
sont finies, les combinaisons d’ensembles sont finis et les dénombrements sont finis, (utilisant seulement 
des additions, soustractions, comparaisons et des appels à un oracle) 
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3.3.4 Fonctions sous-modulaires et VCSP 

L’une des restrictions les plus importantes et les plus étudiées sur les contraintes 
valuées est la restriction des fonctions de valuations à des fonctions sous-modulaires. 
La classe des contraintes valuées ainsi obtenue est la seule classe polynomiale, non 
triviale, des problèmes d’optimisation dans le classement dichotomique des VCSPs 
booléen [16]. Elle est la seule classe polynomiale dans le classement dichotomique 
des Max-CSPs pour les domaines à trois éléments [69]. Pour les domaines arbitraires 
finis, la classe des fonctions sous-modulaire englobe des fonctions constantes, don¬ 
nant lieu à des contraintes à valeurs fixes [29]. La sous-modularité joue également 
un rôle crucial dans le problème Min-coûts dans un digraphe [54]. Ce problème est 
un cas particulier du cadre VCSP. 

La maximisation et la minimisation des fonctions sous-modulaires, éventuelle¬ 
ment sous certaines conditions supplémentaires, ont été largement traitées dans la 
littérature. Il a été montré que la maximisation des fonctions sous-modulaires est NP- 
difficilc [110]. En effet, elle généralise de nombreux problèmes NP-difficiles standards 
tel que le problème de la coupe maximale. En revanche, le problème de minimisation 
des fonctions sous-modulaires (SFM) peut être résolu efficacement, en utilisant l’un 
des algorithmes combinatoires qui ont été proposés dans [65, 98, 109]. La complexité 
temporelle de l’algorithme général le plus rapide [98] pour le problème SFM est de 
0(n 6 + n 5 Ef ), où n est le nombre de variables et Ef est le temps nécessaire pour 
évaluer la fonction f(S ) où S est un sous-ensemble d’éléments quelconque. 

3.3.5 Contraintes sous-modulaire et coupe minimale de graphe 

Les contraintes sous-modulaires jouent un rôle important dans la résolution des pro¬ 
blèmes de satisfaction de contraintes valuées (VCSPs). Il a, déjà, été montré en 
utilisant les résultats de la théorie de l’optimisation combinatoire, que les cas de 
VCSPs avec des contraintes sous-modulaires peuvent être résolus en temps poly¬ 
nomial. Cependant, l’algorithme général a une complexité de 0(n 6 ) et donc peu 
pratique. 

Classe de contraintes sous-modulaires d’arités binaires 

Une fonction binaire </> : D x D' —»■ E est sous-modulaire si et seulement si pour 
toutes u, v E D et pour toutes u',v r E D 'tel que u < v et u' < v', on a 

<j)(u,u') © 4>(v,v') ■< 4>(u,v r ) © 4>(v,v!) (3.2) 

à présent la définition des fonctions d’intervalle généralisé (GI). Ces fonctions forment 
une classe assez restreinte de fonctions binaires. Elles ont une structure très parti¬ 
culière. 
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Définition 24. Soit D un ensemble totalement ordonné, une fonction binaire <f> : 
D 2 -E E est dite fonction d’intervalle généralisé sur D si elle a la forme suivante : 


(p{x,y) 


0 si (x < a) V (y > b ) 
p sinon 


Pour toutes a, b E D et tout p E E cette fonction sera notée 77 ^ ^. 

L’ensemble de toutes les fonctions d’intervalle généralisé sur l’ensemble totalement 
ordonné D sera noté E ci- 


Dans [14], plus précisément le Théorème 4.7, les auteurs ont montré que pour 
tout ensemble de langage de contraintes, valuées, T, défini sur un ensemble fini to¬ 
talement ordonné D, si T ne contient que des fonctions sous-modulaires unaires ou 
binaires, la complexité en temps du VCSP(T) est de 0(n 3 d 3 ), ofi n est le nombre de 
variables et d est la taille du domaine de valeur le plus large. Ceci est possible en 
exprimant ces VCSPs comme des problèmes de calcul de la coupe minimale dans un 
graphe orienté et pondéré. Ce dernier est défini comme suit : 


Définition 25. Soit V = (X, D, C, S) est un VCSP binaire, dont les fonctions sont 
dans Tgi■ Un graphe orienté pondéré, dénoté par G-p associé à V, est défini comme 
suit : 


- Les nœuds de G-p sont : {s,t} U xe x {x v | n = 1,..., \D X \] 

- Les arcs orientés de Gp sont : 

- Pour chaque x E X, il existe un arc de s à x v pondéré par + 00 . 

- Pour chaque x EX, il existe un arc de xq à t pondéré par + 00 . 

- Pour chaque x E X, et chaque v E {1, 2,..., |_D X | — 2}, il existe un arc de 
x v à x v+ \ pondéré par + 00 . 

- Pour chaque contrainte ((x, y), rj? a M E C, il existe un arc de yb à x a -\ 
pondéré par p. Ces arcs sont appelés “arcs de contraintes”. 

Une coupe dans le graphe Gp est un ensemble d’arcs dont la suppression décon¬ 
necte la source s du puits t. Ce qui fait qu’il n’existe plus de chemin orienté qui part 
de la source s et qui arrive à la puits t. La capacité d’une coupe est l’agrégation 
des capacités de ces arcs. Une (s, t)-coupe se définit également par un sous-ensemble 
d’arcs A formant une partition de Gp en deux sous-ensemble de sommets S et T telle 
que s E S et t E T. E 11 d’autre termes, A est l’ensemble de tous les arcs (x, y) dont 
l’extrémité initiale x est dans S et l’extrémité terminale y est dans T. Le problème 
de recherche d’une Coupe Minimale dans le graphe Gp consiste à trouver une coupe 
entre S et T dont la capacité est minimale. 

Par ailleurs, un flot F dans un graphe Gp est un sous-ensemble a d’arcs de Gp 
dont les valuations sont dans M + . F assure la satisfaction de la loi des nœuds (lois 
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de Kirchhoff) en tout sommet de Gp sauf deux. On peut supposer que pour l’un de 
ces sommets, l’origine s, aucun flot ne rentre, tandis que pour l’autre, la destination 
du flot t, aucun flot ne sort. La valeur d’un flot est la somme totale du flot sortant 
de s. Cette somme est égale à celle du flot rentrant dans t par la conservation de 
flot découlant de la loi de Kirchhoff. 

Le théorème flot-max/coupe-min est un théorème bien connu de la théorie des 
graphes. Il révèle que le calcul d’une coupe minimale peut se ramener à un problème 
de flot maximum. 

Théorème 4. (indépendamment P. Elias, A. Feinstein et C. E. Shannon, et L. R. 
Ford, et D. R. Fulkerson, (1956)) Pour tout graphe, tout couple x,y de sommets du 
graphe, et pour toute pondération positive, la valeur maximale d'un flot F de s à t 
est égale à la capacité minimale d'une coupe séparant s de t. 

L’algorithme de FordFulkerson permet de calculer un flot maximal F dans un 
graphe orienté pondéré G. Ce flot maximal correspond à une coupe minimale dans 
ce graphe (voir Théorème 4). En résumé, une coupe minimale correspond au flot 
maximal Fp calculé dans le graphe Gp. L’algorithme de FordFulkerson [39] assure le 
calcule du flot maximal dans Gp et donc la coupe minimale dans ce même graphe. 
Cet algorithme est itératif. En effet, à chaque itération, la solution courante est 
un flot réalisable qui satisfait les contraintes de capacité et l’algorithme cherche à 
augmenter la valeur de ce flot. Considérons un flot possible dont les sous-ensembles 
sont les différents flots associés à chaque arc du graphe. 

On définit le graphe résiduel Gpf qui est un graphe avec capacité et sans flot. 

(j) f (u,v) = (/>(u,v) - f(u,v) 

Notons qu’il peut arriver qu’un flot de u hv est autorisé dans le graphe résiduel, s’il 
ne l’ai pas dans le graphe d’origine : si f(u,v ) > 0 et <f>(v,u) = 0 alors 

u ) = 0(u, u) - f(v, u) = f(u , v) > 0 

Notons par p un chemin de s à t, dans le graphe résiduel Gpf, tel que pour tout 
(u,v) G p, 4>f(u,v ) > 0. Un pseudo-code, de FordFulkerson est présenté dans 
l’Algorithme 16. 

Après une initialisation du flot maximal F* et de / (voir lignes 1 et 2), l’Al¬ 
gorithme de FordFulkerson, calcule, pour chaque chemin p , 4>f(p) : la capacité 
minimale de ce chemin (voir ligne 3). Ensuite, il envoie cette capacité minimale le 
long du chemin en mettant ainsi à jour / (voir ligne 4 et 5). Enfin, il met à jour le 
flot maximal F* (voir ligne 6 ), la valeur de ce dernier sera renvoyée lorsqu’il n’y a 
plus de chemin p. 

Le Théorème 4.8, de [14], montre que la classe polynomiale identifiée dans le 
Théorème 4.7 est maximale. En effet, étant donné T l’ensemble de toutes les fonctions 
sous-modulaires binaires sur un ensemble fini totalement ordonné D, tel que \D\ >2, 
pour toute fonction binaire tjj ^ T, VCSP(r U ÿ) est NP-difficilc. 
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Algorithme 16: FordFulkerson(Gp) : F* 
1 FVO 


pour u € D faire 
pour v € D' faire 

_ f(u,v) = 0 


3 


tant que p / 0 faire 
pour u G p faire 
pour v £ p faire 


_ 4>f{p) = min {<f> f (u,v)} 


4 

5 


pour u £ p faire 
pour v £ p faire 

f(u,v) <- f(u,v) + 4>f(p) 
_ f(v,u ) «- f(v,u) - (j)f(p) 


6 


F* «- F* + <t>f[p) 


3.3.6 Décomposition de contraintes sous-modulaires 

Les résultats de Stanislav Zivny, David A. Cohen, Peter G. Jeavons [122] ont donné 
lieu à plusieurs corollaires. Tout d’abord, ils ont caractérisé précisément quelles fonc¬ 
tions sous-modulaires d’arité quatre peuvent être exprimées par des fonctions bi¬ 
naires sous-modulaires. Ensuite, ils ont identifié une nouvelle classe de fonctions 
sous-modulaires d’arités arbitraires qui peuvent être exprimées par des fonctions bi¬ 
naires sous-modulaires, et donc minimiser efficacement en utilisant une réduction en 
sous-problèmes chacun pouvant être résolu par l’algorithme de la coupe minimale. 
Plus important encore, ces résultats ont impliqué des limites à ce genre de réduction 
et ont établi, pour la première fois, que cette technique ne peut pas être utilisée, en 
général, pour réduire les fonctions arbitraires sous-modulaires. 

Un sous-problème important et très étudié de la SFM est la minimisation loca¬ 
lement défini des fonctions sous-modulaires [19], ou des fonctions sous-modulaires 
avec représentation [34], Stanislav Zivny et Peter G. Jeavons ont appelé ce problème 
de minimisation de fonctions sous-modulaires bornées, SFMb. 

Dans ce scénario, la fonction sous-modulaire à minimiser est définie comme la 
somme d’un ensemble de fonctions qui ne dépendent que d’un nombre borné de va¬ 
riables. 

Ces Problèmes d’optimisation interviennent dans une variété de contextes : 

- Dans le contexte d’optimisation pseudo-booléenne, ces problèmes font interve¬ 
nir la minimisation des polynômes booléens de degrés bornés [8]. 

- Dans le cadre de la vision par ordinateur, ces problèmes permettent de for- 
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muler des problèmes de minimisation de l’énergie de Gibbs [47] ou champs de 
Markov (aussi Conditional Random Ficlds) [77]. 

- Dans le contexte de l’intelligence artificielle, ils ont été étudiés comme des 
Problèmes de Satisfaction de Contraintes Valuées [102] avec des contraintes 
d’arités bornées. 

Un algorithme général pour la SFM peut toujours être utilisé pour résoudre la 
SFMb, mais les particularités de ce problème plus restreint permettent parfois de 
proposer des algorithmes spéciaux plus efficaces. Par conséquent, certaines catégo¬ 
ries de SFMb sont connus pour être résolues plus efficacement que la SFM en général. 

Comme les VCSPs sous-modulaires sont équivalents à des SFMb, et les SFMb ont 
été étudiées dans plusieurs autres domaines de l’informatique, dans [123] les auteurs 
ont utilisé d’autres techniques, inspirées par d’autres domaines, pour résoudre les 
VCSPs sous-modulaires. 

Stanislav Zivnÿ et Peter G. Jeavons [123] focalisent sur les contraintes sous- 
modulaires sur un domaine Boolean {0,1}, qui correspondent précisément aux fonc¬ 
tions définies comme sous-modulaires [16]. Dans [123], les auteurs décrivent un al¬ 
gorithme basé sur une coupe de graphe qui peut être utilisée pour résoudre certains 
VCSPs à contraintes sous-modulaires sur un domaine booléen beaucoup plus efficace 
que le cas général. Certains de ces résultats sont étroitement liés à des cas connus 
efficaces pour les SFMb, et d’autres résultats antérieurs de différents domaines de 
l’informatique. En outre, ces auteurs identifient de nouvelles classes de VCSPs avec 
des contraintes sous-modulaires d’arités arbitraires qui peuvent être résolues effica¬ 
cement. 

En utilisant les résultats de la théorie de l’optimisation pseudo-booléenne 5 , on 
peut identifier plusieurs grandes catégories de contraintes sous-modulaires sur un 
domaine booléen. Ces catégories sont exprimables à l’aide des contraintes binaires 
sous-modulaires, et peuvent donc être minimisées en temps cubique [14]. Dans [122] 
les auteurs ont utilisé les techniques de l’optimisation pseudo-booléenne pour iden¬ 
tifier de nouvelles classes sous-modulaires de VCSPs qui peuvent être résolues effi¬ 
cacement. 

Les résultats de Stanislav Zivnÿ et Peter G. Jeavons [123] ont des conséquences 
directes pour les autres formalismes et des cadres dans lesquels le problème SFMb 
a été étudié. Par exemple, plusieurs problèmes de la vision par ordinateur peuvent 
être exprimés par la minimisation de l’énergie [71]. 

Cohen et al. ont montré que les contraintes de VCSPs sous-modulaires sur un 
domaine quelconque fini ne peuvent pas être réduites au problème de SFM sur une 
famille particulière d’ensembles connue sous le nom de famille de cycle [16]. Ce 
problème est équivalent alors au problème général de SFM [109] nécessitant ainsi 

5. Une fonction pseudo-booléenne est une fonction qui compare n valeurs booléennes à un réel. Le terme 
pseudo-booléen vient du fait que, bien que ces fonctions ne soient pas booléennes, elles font beaucoup penser 
aux fonctions booléennes et y empreintes quelques notions. 
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un algorithme d’ordre 0(n 6 + n 5 Ef). Ce résultat a depuis été généralisé à une large 
classe de contraintes valuées arbitraires à plus de domaines finis connu sous le nom 
de contraintes de tournoi par paire [17]. Une autre approche pour résoudre les cas 
VCSPs avec des contraintes sous-modulaires, basées sur la programmation linéaire, 
peut être trouvée dans [19]. 

3.3.7 Arc consistance flexible optimale (OSAC) 

L’algorithme d’arc consistance flexible optimale (Optimal Soft Arc Consistency ou 
OSAC [20]), est basé sur la programmation linéaire. Cet algorithme calcule un en¬ 
semble de transformations polynomiales qui préservent l’équivalence [18]. Il permet 
d’avoir la plus grande valeur possible de la borne inférieure fournie par cet algo¬ 
rithme, qui semble être trop coûteuse pour être maintenue durant la recherche. 

Dans son rapport d’habilitation Simon de Givry a considéré que la critique la 
plus importante de OSAC est la taille de son système linéaire 0(ed + n ) variables et 
0(ed r + nd) contraintes linéaires où r est l’arité maximale des fonctions. Il a aussi 
indiqué que depuis 1976 Schlesinger [107], dans le contexte du traitement d’images, 
a proposé une version équivalente à OSAC qui n’exploite pas les contraintes dures. 
Et en 2002, dans le cadre de formules pseudo-booléennes quadratiques (équivalent 
à Max-2SAT), l’auteur dans [8] a proposé une reformulation optimale sous forme 
d’un algorithme spécialisé de maximisation de flot. Finalement, notons qu’il est déjà 
connu que OSAC permet de résoudre les VCSPs avec des fonctions de coût sous- 
modulaire [19]. 

3.3.8 Arc Consistance Virtuel (VAC) 

Pour faire face à la grande taille du système linéaire de OSAC, Cooper et al. [22] 
ont proposé un algorithme de reformulation plus efficace. Il s’agit de l’arc-consistant 
virtuel (Virtual Arc Consistency (VAC)). VAC est une version affaiblie de OSAC. 
Il permet de détecter si, dans un VCSP, il existe une séquence de transformations 
locales qui puisse augmenter la borne inférieure (ft 0 du VCSP. Malheureusement, 
cette version plus rapide ne garanti pas l’optimalité de la borne inférieure trouvée. 

Le principe de cette nouvelle version est d’appliquer l’arc-consistance classique 
sur un CSP dérivé du VCSP d’origine, ce CSP est dénoté bool(V). Ce CSP est 
obtenu en remplaçant chaque fonction de coût autre que 0 0 par une fonction dûre 
ayant la même portée et obtemre en remplaçant les valuations supérieures à _L sont 
remplacées par T [21, 22] 

Définition 26. Un VCSP V est arc-consistant virtuel (Virtual Arc Consistency 
(VAC)) si le renforcement d’arc-consistance dans le CSP Bool(V) ne produit pas un 
domaine vide. 

Dans [22] les auteurs ont décri un exemple qui montre la limite de VAC puis¬ 
qu’il peut exister une infinité de séquences qui font augmenter la borne inférieure de 
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manière infinitésimale. Ils ont réussit à présenter un algorithme approché, dont la 
terminaison est assurée, pour les VCSPs binaires. Cet algorithme, dénoté VAC e , est 
relativement efficace en temps, 0(ed 2 k e) où e est le coût rationnel non nul minimal, 
et en espace 0(ed). 

Dans son rapport d’habilitation, Simon de Givry a indiqué qu’il a commencé 
à comparer ces travaux avec l’existant en recherche opérationnelle en prenant plu¬ 
sieurs benchmarks dont les problèmes WAREHOUSE et de couverture de minimum 
d’ensembles (TAGSNP et Steiner Triple Covering Problem (STEINER) [99]). Selon 
cette comparaison l’auteur a déduit que la qualité des bornes inférieures obtenues 
par VAC et OSAC est souvent proche de celle produite par la PLNE appliquée sur 
une formulation directe du problème, mais en un temps plus long. VAC s’appuie sur 
l’arc-consistance classique qui a de bonnes propriétés, telles que : 

- son indépendance de l’ordre des opérations. 

- ses classes polynomiales dont celle où la structure du réseau est un arbre et 
EDAC* et celle des contraintes max-closed correspondant à des fonctions sous- 
modulaires en optimisation. 

Quand un problème V est arc consistant virtuel, le problème Bool(V ) associé 
a un e fermeture par arc consistance n’ayant pas de domaine vide. Cela permet à 
VAC d’hériter différentes classes de problèmes polynomiaux qui sont résolus par 
arc-consistance dans le cadre CSP. Par exemple, VAC peut résoudre les problèmes 
de minimisation sous-modulaire, sur la structure de valuation Q + [16]. Où Q + = 
( Q + U {oo}, +, <) et Q + représente l’ensemble des nombres rationnels non-négatifs. 
Cooper et al [22] ont montré qu’un VCSP dont les fonctions de coût sont toutes 
sous-modulaires à une permutation près des domaines de valeurs, peut être aussi 
résolu par VAC, donc en temps polynomial. 

3.3.9 Exploitation de la (sous-)modularité 

La complexité algorithmique de la recherche de la solution optimale d’un VCSP a été 
étudiée dans plusieurs travaux de recherche. Et les VCSPs impliquant les fonctions 
de valuations sous-modulaires est une des classes polynomiales les plus étudiées. 
Néanmoins, les VCSPs résultant d’instances réelles ne sont pas limités à des fonc¬ 
tions sous-modulaires. 

Dans le cas où le problème implique des fonctions de valuations qui ne sont pas 
sous-modulaires, nous montrons comment un tel problème peut être décomposé en 
un ensemble de sous-problèmes formés, exclusivement, de fonctions de valuations 
sous-modulaires (voir chapitre 6). 

Dans [50], les auteurs ont proposé une théorie s’appuyant sur des permutations 
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des valeurs des domaines de variables dont l’objectif est de réduire les CSPs à des 
CSPs dans un langage de contraintes soluble en temps polynomial, à savoir le lan¬ 
gage de contraintes max-fermées [67]. Dans une troisième contribution notre objectif 
est d’étudier le problème de la découverte de permutations de domaines qui trans¬ 
forment les instances des VCSPs en des instances dans des classes solubles en temps 
polynomial, notamment, en des VCSPs sous-modulaires. 

Pour ce faire, nous montrons comment la théorie de la réduction de domaine 
par permutation proposée dans [50] peut être adaptée et étendue pour résoudre plus 
efficacement les VCSPs. 

La modularité est une condition plus forte que la sous-modularité, mais en com¬ 
pensation, exploiter la modularité pour résoudre les VCSPs, garanti que les VCSPs 
impliquant de telles fonctions peuvent être identifiés et résolus plus efficacement. 

Dans une deuxième contribution nous présentons un algorithme de décomposi¬ 
tion de domaine pour résoudre des VCSPs binaires qui ne sont pas limités à des 
fonctions de valuations sous-modulaires. Le processus de décomposition produit un 
ensemble de sous-problèmes dont les fonctions de valuations sont modulaires et dont 
la résolution est polynomiale. 

Plus précisément, pour une classe de VCSPs binaires où toutes les fonctions de valua¬ 
tion sont exclusivement modulaires, nous proposons un algorithme d’identification 
de complexité 0(ed 2 ) et un algorithme de résolution de complexité 0(ed) [93]. Nous 
proposons également une méthode de décomposition de VCSP s’appuyant sur cette 
classe polynomiale. 


3.4 Substituabilité directionnelle 

La notion de substituabilité de voisinage [42] a été proposée et utilisée comme un 
moyen de filtrage des domaines de valeurs des variables des CSPs binaires. Plusieurs 
variantes de cette notion ont été proposées, parmi lesquelles la notion de substitua¬ 
bilité directionnelle [91]. 

La substituabilité directionnelle est une forme faible de la substituabilité de voisinage 
qui, air départ, visait à améliorer la résolution des CSPs binaires. On part de l’idée 
que même si deux valeurs ne sont pas voisinage substituables, elles peuvent l’être si 
on restreint le voisinage en se référant à un ordre sur les variables. La substituabilité 
directionnelle ne peut pas être utilisée comme un moyen de ûltrage des domaines 
des variables comme c’est le cas pour la substituabilité de voisinage. Toutefois, on 
peut l’utiliser comme un moyen de décomposition des domaines des variables en 
de chaînes de valeurs directionnellement substituables qui peuvent être considérées 
simultanément lors de la résolution de problèmes. 
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Dans une première contribution, (voir Chapitre 4), nous proposons une méthode 
de décomposition de domaine en chaines de valeurs dites directionnellement substi¬ 
tuables pour résoudre les VCSPs en passant par leurs duales : les VCSPs binaires 
avec des Contraintes binaires Dures [92], 

Nous proposons une extension du travail présenté dans [91] qui consiste à géné¬ 
raliser davantage la notion de substituabilité directionnelle en considérant, comme 
référence, une orientation du graphe d’inconsistance au lieu d’un ordre total sur les 
variables. Pour garantir l’optimalité de la solution lors de la résolution des Problèmes 
de Satisfaction de Contraintes Valuées à Contraintes binaires Dures, nous introdui¬ 
sons des conditions supplémentaires de substituabilité. Nous utilisons la notion de 
substituabilité directionnelle proposée pour les VCSPs-CbC pour résoudre, de ma¬ 
nière indirecte, les VCSPs. 


3.5 Conclusion 

Motivés par le fait que le dual d’un VCSP est un VCSP binaire à Contraintes bi¬ 
naires Dures (VCSP-CbC), nous présentons dans le chapitre 4 une méthode de 
décomposition, basée sur la substituabilité directionnelle permettant de résoudre les 
VCSPs-CbC. Comme le dual d’un VCSP peut avoir un nombre de variables qui est le 
carré du nombre de variable dans le problème primai, il est plus utile de trouver des 
solutions algorithmiques de plus en plus puissantes et permettant de résoudre direc¬ 
tement le VCSP d’origine. Combiner les notions de modularité et de sous-modularité 
avec une méthode de décomposition pourra nous permettre de trouver de nouvelles 
solutions algorithmiques pour résoudre les VCSPs. 

Nous présentons dans les chapitres 5 et 6 deux méthodes de décomposition, la 
première s’appuie sur la modularité et la deuxième exploite la sous-modularité. Ces 
deux méthodes permettent de résoudre directement les VCSPs binaires. 
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substituabilité directionnelle 
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Résumé 

Motivés par le fait que le dual d'un VCSP est un VCSP binaire à Contraintes binaires Dures 
(VCSP-CbC), nous proposons, dans ce chapitre, des solutions efficaces pour résoudre ces derniers. 
Pour ce faire, nous généralisons la notion de substituabilité directionnelle, (proposée initialement 
pour les CSPs binaires), au cadre des VCSPs-CbC. 

Notre contribution, dans ce chapitre, consiste en une méthode de décomposition des VCSPs 
en passant par leurs duaux qui sont des VCSPs-CbC. Cette décomposition produit un ensemble 
de VCSP-CbC particuliers dont la résolution est polynomiale. 

Des résultats de simulation, sur plusieurs VCSPs-CbC, modélisant des problèmes d’ordon¬ 
nancement, montrent qu'une variante de l'algorithme du Branch-AND-Bound qui exploite la 
substituabilité directionnelle est souvent plus efficace que la version classique du Branch-AND- 
Bound. 

4.1 Introduction 

Ce chapitre a fait l’objet de la publication [92], 

Dans ce chapitre nous nous intéressons aux méthodes de résolution de VCSPs 
dites complètes : c’est à dire celles qui permettent de trouver une solution optimale 
si cette dernière existe. 

L’algorithme du BRANCH-AND-BOUND [78] est sans doute l’algorithme complet 
le plus connu et le plus utilisé pour résoudre les VCSPs-CbC. C’est une procédure 
d’exploration en profondeur d’abord munie d’une fonction heuristique servant à éva¬ 
luer le coût de la solution complète que l’on peut espérer atteindre par l’extension de 
la solution partielle courante. Le rôle de la fonction heuristique est d’éviter l’explora¬ 
tion de branches de l’arbre de recherche qui ne mènent pas à des solutions meilleures 
que celles déjà trouvées. 

Par ailleurs, la notion de substituabilité [42] constitue l’un des outils les plus 
efficaces qui permettent le filtrage des domaines des variables des CSPs binaires 
dans le but est d’accélérer leur résolution. 

Soient a et b deux valeurs appartenant à une variable x d’un CSP binaire. On 
dit que b est substituable à a si le fait de substituer b à a dans une solution du 
problème, produit une nouvelle solution qui, au moins, a un coût meilleur sinon 
identique. Dans une telle situation, il est possible de supprimer a du problème tout 
en étant sûre d’avoir gardé une solution qui est, au moins, de qualité équivalente. 

En se référant à [42], une valeur a d’une variable x est voisinage substituable 
(VS) à une valeur b de x si et seulement si l’ensemble des valeurs incompatibles avec 
a sont inclus dans l’ensemble des valeurs incompatibles avec b. 
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Plusieurs variantes de cette notion ont été proposées, parmi lesquelles la notion 
de substituabilité directionnelle [90]. Cette dernière est une forme faible de substi¬ 
tuabilité qui, au départ, visait l’amélioration la résolution des CSPs binaires. On 
part de l’idée que même si deux valeurs d’une variable x ne sont pas voisinage sub¬ 
stituables, elles peuvent l’être si on restreint le voisinage de x. La restriction du 
voisinage d’une variable peut être effectuer de plusieurs façons. Dans [90], cette res¬ 
triction a été effectuée en se référant à un ordre total sur les variables. Le voisinage 
la variable x est alors limité aux variables qui sont liées à x par une contrainte est 
qui, en plus, arrivent avant x dans l’ordre considéré. 

La substituabilité directionnelle 11 e peut pas être utilisée comme un moyen de 
filtrage des domaines des variables comme c’est le cas pour la substituabilité de voi¬ 
sinage. Toutefois, 011 peut l’utiliser comme un moyen de décomposition du problème, 
comme nous le montrerons dans la suite de ce chapitre. 

Notre contribution, dans ce chapitre, consiste en une méthode de décomposition 
des VCSPs-CbC qui s’appuie sur la substituabilité directionnelle. Nous proposons 
une généralisation de la notion de substituabilité directionnelle qui peut être appli¬ 
quée aux VCSPs-CbC et qui, en plus, restreint les voisinages des variables en utilisant 
une référence plus flexible qu’un ordre total sur les variables. Cette référence est une 
orientation d’un graphe lié au VCSP. Nous montrons enfin comment utiliser la notion 
de substituabilité directionnelle pour résoudre les VCSPs en passant par leurs duaux. 

Ce chapitre est organisé comme suit : dans la Section 4.2, nous introduisons les 
débilitions et les notations utilisées. Nous déhnissons, par la suite (Section 4.3), la 
notion de substituabilité directionnelle. Nous présentons, dans la même section, une 
version de l’algorithme du BRANCH-AND-BOUND qui utilise la substituabilité di¬ 
rectionnelle pour décomposer, en cours de résolution, des VCSPs-CbC. La Section 
4.5 présente un algorithme d’orientation du graphe d’inconsistance utilisée, comme 
référence, pour la détermination de la substituabilité directionnelle. Dans la Section 
4.6, nous montrons comment nous pouvons utiliser la notion de substituabilité di¬ 
rectionnelle pour résoudre les VCSPs. Enhn, dans la Section 4.7, nous reportons les 
résultats d’une expérimentation montrant l’apport de la notion de substituabilité 
directionnelle pour la résolution des VCSPs binaires. 

4.2 Définitions et notations 

Tout d’abord rappelons qu’un VCSP-CbC se déhnit par la donnée d’un ensemble de 
variables, chaque variable pouvant prendre une valeur parmi un ensemble de valeurs 
possibles appelé domaine. Puis, par un ensemble de contraintes binaires exprimant 
des restrictions sur les combinaisons de valeurs permises. Enhn, un ensemble de 
contraintes unaires exprimant des préférences sur les valeurs. 

Dans un VCSP-CbC, les contraintes binaires doivent être obligatoirement satis- 
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faites. Ce sont donc des contraintes binaires dures. Par contre, les contraintes unaires 
sont flexibles. Elles permettent d’exprimer des préférences entre les solutions qui sa¬ 
tisfont à toutes les contraintes binaires dûres. Ainsi, l’ensemble des valuations pour 
les contraintes binaires est la paire {_L,T} 12 . L’ensemble des valuations pour les 
contraintes unaires est un ensemble plus vaste totalement ordonné par A. L’opéra¬ 
teur d’agrégation est noté © comme d’habitude. Formellement, dans le cas binaire, 
un VCSP - CbC ce défini comme suit : 

Définition 27. Un Problème de Satisfaction de Contraintes Valuées à contraintes 
binaires dures (VCSP-CbC) est défini par un quadruplet (X, D,C, S) tel que : 

1. X est un ensemble de variables. 

2. D est un ensemble de domaines de valeurs. 

3. S = ( E , ©, A) est une structure de valuation. 

4. C = C 1 CC 2 où 

- ((x) , 0) G C 1 , =£> (j> x G {D x —> E}. 

- ({x,y) ,0) G C 2 , 4> xy G {D x x D y -* {_L,T}}. 


Rappelons que le but est de trouver un tuple t tel que <f-p(t ) = Y^ceC^cif 0 cr c ) 
est minimal. Rappelons aussi que, deux variables x et y reliées par une contrainte 
binaire, notée ((x,y) ,0), sont dites voisines. La valeur a G D x , est compatible avec 
b G D y si 0(a, b) — _L. Dans ce cas, on dit que b est un support de a. Rappelons que 
si chaque valeur du problème a, au moins, un support dans le domaine de chaque 
variable voisine alors le problème est dit arc-consistant. 

En se référant à [42], dans un CSP, une valeur a d’une variable x est voisinage 
substituable (VS) à une valeur b de x si et seulement si l’ensemble des valeurs 
incompatibles avec a sont inclus dans l’ensemble des incompatibles avec b. 

Définition 28. Soit V un CSP binaire, x une variable deV et a, b G D x . La valeur 
a est dite voisinage substituable à la valeur b (notation a fi b) si et seulement si, 
pour tout ((x, y) , 0) G C 2 et pour tout a' G D y , on a 0(a, a') fi 0(6, a') 

Cette définition, qui est valable uniquement pour les CSPs binaires, peut être 
étendue pour tenir compte des fonctions coût unaires de la manière suivante : 


Définition 29. Soit V un VCSP-CbC, x une variable deV et a,b G D x . a G D x est 
dite voisinage substituable à b G D x (notation a fi b) si et seulement si 4> x (a) fi <f> x (b) 
et, pour tout ((x, y) , 0) G C 2 et pour tout a' G D y , on a 0(a, a') fi 0(6, a') 

1. Rappelons que _L est l’élément neutre et T est l’élément absorbant 

2. Rappelons aussi que l’ordre total sur ces deux éléments est tel que IV. 
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Nous avons utilisé le même symbole (^) pour désigner la relation de substitua¬ 
bilité de voisinage et celle de l’ordre total sur les valuations, car la première dépond 
essentiellement de la deuxième. 

On peut aisément vérifier qu’à partir de toute solution dans laquelle x a pour 
valeur 6, on peut déduire une solution de qualité, au moins, égale si l’on substitue 
a à 5. Par conséquent, b peut être supprimée du domaine de x sans que l’on perde 
toutes les solutions optimales du problème. 

4.3 Substituabilité directionnelle pour les VCSPs-CbC 

La substituabilité directionnelle [90] est une forme faible de substituabilité de 
voisinage [42]. Au départ, cette notion a été définie pour les CSPs binaires. Dans 
le présent chapitre, nous généralisons la notion de substituabilité directionnelle de 
manière à tenir compte des fonctions coût unaires qui interviennent dans les VCSP- 
CbC. 

Définition 30. Le graphe d’inconsistance d’un VCSP-CbC V est un graphe simple 
GI(V ) dans lequel les sommets correspondent aux valeurs des variables et les arêtes 
relient les paires de sommets qui représentent des valeurs incompatibles. 


La substituabilité directionnelle est définie en utilisant, comme référence, une 
orientation du graphe d’inconsistance du VCSP-CbC. 

Définition 31. Une orientation A du graphe d’inconsistance d’un VCSPs-CbC est 
une affectation d’une orientation à chaque arête {a, b} du graphe donnant lieu, ou 
bien à l’arc (a, b) ou bien à l’arc (b, a). 

La notion de substituabilité directionnelle est alors une relation binaire définie 
comme suit : 

Définition 32. Soit V = ( X,D,C,S ) est un VCSP-CbC, x G X et a, b G D x . On 
dit que a est directionnellement substituable à b par rapport à A (notation a P a b) 
si (p x (a) P (f) x {b) et pour tout (( x , y ), 0) G C 2 et a' G D y , on a 

(a, a') G A =>■ {b, a') G A A 0(a, a') P 0(6, a') 

Exemple 4.3.1. Si l’on considère l’orientation du Graphe d’inconsistance indiquée 
par les arcs à la Figure 4.1 (a) alors les deux conditions de la Définition 32 sont 
vérifiées, d’où a Pa b. Dans la Figure 4.1(6) <f> xy (a,a') ^ <f> xy (b,a') et dans la Figure 
4.1(c) on a 4> x (a) ^ <fi y {b). 

Etant donnée une orientation A du graphe d’inconsistance, la relation P a définit 
un préordre sur le domaine de chaque variable. 

Propriété 5. -<a est un préordre sur D x . 
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Figure 4.1 - Cas où a Aa b 


Preuve. <a est trivialement réflexive. Montrons qne Aa est aussi transitive. À cette 
fin, supposons l’inverse et procédons pour obtenir une contradiction. 

Soit u, v, w G D x tel que u Aa v et v Aa u> mais u ^a w - Cela implique que, 
pour tout u' G D y , nous avons 


(f) x (u) A 4> x (v) A [(u, u') E A=ï (v, u') e A A 4> xy (u, u') A <t) xy (v , u')] (4.1) 

<t>x{v) A (j> x {w) A [(v,u') G A =>- (w,vf) G A A(j> xy (v,u') A (p xy (w,u')\ (4.2) 

comme u ^a w et A est un ordre total, il doit exister u' G D y tel que 

(j) x (u) y (j) x {w ) V [(w, u') e A A (w, u') £ A y <t> xy {u, u') y (p X y(w, u')} (4.3) 

A partir de (4.1) et (4.2) nous obtenons 

4> x {u) A (f> x (w) A [(w, ii')^4V ( w , u') G A A (j) xy (w, u') y (f> xy {u, u ')] 

qui est en contradiction avec (4.3). 

Ainsi, y a est transitive. 
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□ 

La relation binaire est définie sur D x comme suit : u ~a v si et seulement si 
u -<a n et n -<a u. 

Propriété 6. est une relation d’équivalence sur D x . 

Preuve. La preuve est découle de la Propriété 5. 


□ 

Ainsi, chaque domaine D x peut être subdivisé en des sous-ensembles D x = D x i U 
... U D xs tel que les éléments de chaque D x ^, k = 1 s sont tous deux à deux 
comparables, c’est-à-dire, que pour tout a, b G D X}k , on a a ^a b ou b ^a a. Chaque 
D x a est donc une chaine de valeurs totalement ordonnée par -<a- 

Dans chaque chaine D Xtk , on peut distinguer le sous-ensemble des éléments mi¬ 
nimaux. 


min {D Xtk ) = {a G D Xtk j V b G D x>k , a b} (4.4) 

4.4 Utilisation de la substituabilité directionnelle 

Rappelons que, nous nous intéressons aux méthodes de résolution de VCSPs- CbC 
dites complètes, c’est à dire celles qui permettent de trouver une solution optimale si 
cette dernière existe. L’algorithme du Branch-and-bound est l’algorithme com¬ 
plet le plus communément utilisé pour résoudre les VCSPs-CbC. Cet algorithme 
explore l’espace de recherche en effectuant une recherche arborescente en profon¬ 
deur d’abord 3 . La recherche arborescente en profondeur d’abord, par son caractère 
systématique, assure que les algorithmes qui l’utilisent sont complets. Les nœuds de 
l’arbre de recherche sont examinés dans un ordre particulier en utilisant des heuris¬ 
tiques de sélection de variable et de valeur. 

Les variables sont instanciées une par une, jusqu’à obtenir une solution ou jus¬ 
qu’à ce qu’une variable ait un domaine vide. Lorsque c’est le cas, on revient sur la 
dernière décision prise et on effectue un autre choix de valeur, si ceci est possible. 
Dans le cas contraire, on remonte plus haut dans l’arbre de recherche. Cette tech¬ 
nique est appelée Backtrack chronologique. Lorsque le nombre de solutions optimales 
est important, cette méthode permet de trouver une de ces solutions en des temps 
raisonnables. Par contre lorsque ce nombre est faible, cette approche devient coû¬ 
teuse. D’ou le besoin de l’optimiser. 


3. Une recherche arborescente dépend de la nature du problème et de l’espace des solutions : si dans 
un problème le nombre de solutions est très important, on cherche, selon un objectif donné, la meilleure 
de ces solutions. 
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Au niveau de chaque nœud de l’arbre de recherche, on peut procéder de deux 
manières différentes qui sont les suivantes : 

- La technique d’énumération : 

- Choisir une variable x ; 

- Choisir une valeur v G D x ; 

- Brancher avec x = v et x ^ v. 

Si toutes les variables ont une valeur unique (succès) alors une solution est 

trouvée. Si une variable admet un domaine vide (échec) alors un retour arrière 

(backtrack) est appliqué. 

- La technique de partitionnement consiste à : 

- Choisir une variable x ; 

- Choisir S Ç D x en vérifiant certaines conditions ; 

- Brancher sur v G S et v 0 S. 

Si toutes les variables ont une valeur G S (succès) =>- solution trouvée. Si une 

variable admet un domaine vide (échec) =>■ retour arrière (backtrack). 

Dans ce qui suit nous utilisons la notion de substituabilité directionnelle comme 
une technique de partitionnement de VCSP-CbC en plusieurs VCSPs-CbC polyno¬ 
miaux. 

La proposition suivante identifie une classe polynomiale de VCSPs-CbC qui s’ap¬ 
puie sur la relation de substituabilité directionnelle. 

Propriété 7. Soit V un VCSP-CbC et soit A est une orientation de GI{V). Si 
chacun des domaines de valeurs de V est une chaine non vide de valeurs direction- 
nellement substituables par rapport à A alors une solution optimale de V peut être 
trouvée en un temps polynomial. 


Preuve. On commence tout d’abord par rendre V arc-consistant. On montre, en¬ 
suite, qu’en sélectionnant un élément minimum (voir (4.4)) de chaque domaine 
de valeur, on obtient une solution optimale. Ce qui veut dire que tout n-uplet 
t G U xeX min(D x ) est une solution optimale. En se référant à la Définition 32 et 
air fait que la fonction 0 est calculable en temps polynomial, cette sélection peut 
être effectuée en temps polynomial. 

Supposons que t G II xe x min(Z? x ) n’est pas une solution optimale. Ceci implique 
que t est inconsistant ou que 0p(£) n’est pas minimum. 

Commençons par supposer que t est inconsistant. Donc t doit contenir, au moins, 
une paire de valeurs incompatibles. Soit a G min(D x ) et b G min(D y ) une telle 
paire. Puisque a et b sont incompatibles, {a, b} doit être une arrête de GI('P). Par 
conséquent, on doit avoir (a, b) & A ou (6, a) G A. 
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Supposons, sans perte de généralité, que (a, b) G A, (sinon on pourra raisonner 
sur b au lieu de a et obtenir le même résultat). Il s’ensuit que f xy (a,b) = T, et 
puisque a G min(D a .) alors pour tout a' G D x , on doit avoir <j) xy (a',b ) = T. Ce qui 
veut dire que b n’a pas de supports dans D x et donc que V n’est pas arc-consistant, 
d’où une contradiction. 

Supposons, à présent que (p-p(f) n’est pas minimum, donc, qu’il existe t! G 
Il xexD x tel que évif’) fvif). Puisque les valeurs de la fonction sont obte¬ 

nues à partir de celles de <f x& xif I x) en utilisant un opérateur monotone, il doit 
exister x G X tel que t f x = a At' l x = a' et f> x (a') -< <j> x (a). Il en résulte que 
a ^ min(H 3 .), d’où une contradiction. 


□ 

Ainsi, les sous problèmes considérés lelong d’un même chemin de l’arbre de re¬ 
cherche sont des réductions du problème initial que l’on peut définir de la manière 
suivante : 

Définition 33. Un problème V = (. X,D',C',S ) est une réduction d'un problème 
V = (X, D , C, S) (notation V’ ÇP ) si 

- D' est tel que D' x Ç D x , pour tout x G X, 

- C' = {(a, fi) | (a, fi) G C et fi est la restriction de fi à II xea D f x }. 

Exemple 4.4.1. Considérons le problème d'ordonnancement d’atelier (Job-Shop) 
classique décrit dans la Figure 4-2. Il s’agit de planifier l’exécution de deux jobs J\ 
et J 2 sur trois machines Mi,M 2 et M 3 . Chaque job Ji se compose de trois tâches 
Ti,i,T i: 2 et Tj )3 . Une tâche T tJ est dédiée à être exécutée sur la machine Mj. Pour 
cet exemple, on suppose que la durée de toutes les tâches est fixée à une unité de 
temps, saufT 2 ,2 qui consomme deux. Les problèmes du type job-shops font intervenir 
deux types de contraintes : des contraintes de précédence qui imposent que les tâches 
d’un même job soient exécutées les une après les autres ainsi que des contraintes 
de ressource qui empêchent qu’une machine exécute plus d’une tâche à la fois. Ré¬ 
soudre un tel problème revient à affecter un temps de début d'exécution à chacune 
des tâches de façon à satisfaire les contraintes de précédence et de ressource tout en 
minimisant le temps total d’exécution de toutes les tâches (makespan). 

Une modélisation possible des problèmes job-shops en termes de VCSP-CbC 
consiste à associer une variable à chacune des tâches. Ainsi, pour notre exemple, on 
obtient un problème impliquant six variables X\,... ,x 6 . Les domaines de valeurs des 
variables représentent une discrétisation du temps maximum impartie à l’exécution 
de toutes les tâches. Pour le présent exemple, on suppose que les différentes tâches 
peuvent commencer à des instants représentés par les entiers de 0 à 4. 

Les contraintes impliquées sont des contraintes binaires de précédence et de res¬ 
source. On en trouve en tout sept contraintes (voir Figure 4-%)- La fonction coût à 
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minimiser est celle définie par 

(fit) = maxla + duréeixA} 

l<i<6 

où t désigne une instanciation de toutes les variables du VCSP-CbC. En tenant 
compte des contraintes de précédence, la fonction <f> peut être calculée en prenant le 
maximum uniquement sur la paire {£ 3 , xfi\ au lieu de X. 


On en déduit une définition simplifiée des fonctions coût unaire 


, , , [ a + durée(xi) 

Ma) = { 0 


si i G {3, 6 } 
sinon 


(4.5) 


On obtient donc (fit) = max® = 1 (0j(t)), où rriax est bien un opérateur monotone. 
Après application d’un algorithme d’arc-consistance, on obtient le problème dont le 
graphe d’inconsistance est représenté dans la Figure f.3. Dans cette même Figure 
on peut également voir une orientation du graphe d ’inconsistance. 

En se référant à cette orientation, on obtient les ensembles des conflits direc¬ 
tionnels donnés dans le Tableau f.l. En examinant ce tableau, on constate que tous 
les domaines de valeurs sont des chaînes de valeurs DS sauf D 2 qui est composé 
de deux chaines : {(£ 2 ,1), (£ 2 , 2)} et {(£ 2 , 3 )}. En réduisant le domaine de x 2 à 
{(£ 2 ,1), {x 2 , 2)} puis à {(£ 2 , 3)}, on obtient deux sous-problèmes qui, d’après la pro¬ 
position 5, ne contiennent que des chaines de valeurs DS comme domaine. D’après 
la Propriété 1, ces deux sous-problèmes peuvent être résolus sans retour sur trace. 



contrainte de ressource Q 1 unité 

contrainte de précédence 2 unités 


Figure 4.2 - Un problème du type job-shop comportant deux jobs et six tâches. Toutes 
les tâches sont supposées avoir une durée d’une unité de temps, sauf T 2)2 qui en consomme 
deux. 


4.4.1 Algorithme de DS partition 

Le préordre définit par la relation peut être utilisé pour partitionner les 
domaines de valeurs des variables en de chaines de valeurs dir ect!o nn ellement sub- 
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x2 x3 



4 

3 

3 

4 


Figure 4.3 - Graphe d’inconsistance du job-shop décrit dans la Figure 4.2 après applica¬ 
tion d’un algorithme d’arc-consistance. On peut également voir une orientation du graphe 
d’inconsistance. 



0 

1 

2 

3 

4 

x\ 

(X 4 , 0) 

0 

0 



X 2 


(xi, 1) 

( xi, 2 ) 

(ah, 2) 

( x 3 , 3) 
(X 6 , 3) 


X3 




0 

0 

X 4 

0 

(xi, 1) 
(X5, 1) 




X5 


0 

i x 3, 3) 
(x 6 ,3) 



X 6 




0 

0 


Tableau 4.1 - Conflits directionnels des valeurs arc-consistantes du VCSP-CbC décrit dans 
la Figure 4.3. 

stituables (DS). 

Le préordre ^a et la relation d’équivalence ~ A induisent un ordre partiel <a 
sur l’ensemble D x \ ~a des classes d’équivalence de tel que [a] [6] si et 

seulement si a < A b, où [a] désigne la classe d’équivalence de la valeur a. 

En général, étant donné un ensemble partiellement ordonné E 1 , on peut avoir 
plusieurs partitions de E en chaines d’éléments totalement ordonnés. En théorie des 
ensembles, la partition optimale en chaines d’un ensemble partiellement ordonné 
est connue sous le nom de partition en chaines de Dilworth (DCP) [30]. C’est une 
partition qui contient le nombre minimum de chaines parmi toutes les partitions 
possibles. La taille d’une telle partition (i.e., le nombre de chaines de la partition), 
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définit la largeur de l’ordre partiel. Le problème qui consiste à trouver la DCP d’un 
ensemble partiellement ordonné peut être résolu en un temps polynomial en l’ex¬ 
primant comme un problème de recherche d’un couplage maximum dans un graphe 
bipartite [61]. 

Motivé par le fait qu’à chaque nœud de l’arbre de recherche, l’algorithme de 
résolution aura autant de choix que de chaines dans la DCP du domaine de la va¬ 
riable courante, nous proposons de calculer une DCP à chaque nœud de l’arbre de 
recherche. Par cette stratégie, on cherche à minimiser la taille de l’arbre de recherche 
que l’algorithme aura à explorer. 

La première étape du calcul de la DCP (voir la Fonction 17) consiste à déterminer 
la relation ~ A . Ceci peut être accompli en O (nd?) étapes en utilisant l’algorithme 
décrit dans [42] ou celui proposé dans [90]. On en déduit les classes d’équivalences 
D x j en 0(d) étapes. L’étape la plus coûteuse est celle du calcul de la relation 
<, 4 . Au pire des cas, d(d — l)/2 tests d’inclusion entre des paires d’ensembles de 
conflits directionnels sont nécessaires. Chaque test d’inclusion peut être accompli en 
0(nd ) puisque chaque ensemble de conflit directionnel contient au maximum d(n — 1) 
éléments. D’où une complexité de O (nd 3 ) pour cette étape. 



Xi 


x. 



Figure 4.4 - La partition en chaine de Dilworth (DCP) 


Ensuite, l’algorithme construit un graphe bipartite G = (V, U , E) tel que V = U = 
D x / et l’ensemble des arêtes E contient une arête ([a], [b]) si et seulement si 
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[a] <^4 [b], La construction de G nécessite O (d 2 ) étapes. Un algorithme de couplage 
maximum est alors appliqué à G. On utilise l’algorithme décrit dans [61] qui s’exé¬ 
cute en 0 (d 2 5 ) étapes. 

Les chaines de la DCP sont extraites du couplage maximum en incluant les élé¬ 
ments de [a] et ceux de [ 6 ] dans une même chaine chaque fois que l’arête ([a], [ 6 ]) 
fait partie du couplage maximum. Ceci demande O (d 2 ) étapes. D’où une complexité 
totale de O (nd 3 ). 

Notons par S x une chaîne de la DCP de D x obtenue selon le préordre U. 4 - Alors, 
nous avons ce qui suit 

Propriété 8. -<a est un préordre total sur ô x . 


Preuve. -<a est réflexive et transitive sur D x alors il est aussi réflexive et transitive 
sur chaque sous ensemble de D x , et en particulier sur ô x . De plus, ^a est total sur 
ô x du fait que ô x est une chaîne, et donc pour toutes u,v G S x , nous avons u -<a v 

OU V -<A U. 


□ 


Algorithme 17: DSPartition(D x , A, V) 


// avec V = (A, D, C, S) 

~A «— DirInterchangeable(ZU, A, V)\ 

C ExtractEqClass (~a,D x ,V) -, 

// C = D x j ~A 

<A<r- DirSubstituable(D x / ~ a ,V); 

G BlPARTITEGRAPH(Da;/ 

M a- CouplageMaximum(G); 

V x a- ExtractChaines(D 2; , ~a, M); 
renvoyer(P x ) ; 



4.4.2 Exemple 

Reprenons l’exemple du paragraphe 4.4.1. Comme nous l’avons déjà constaté, après 
application d’un algorithme d’arc-consistance, tous les domaines de valeurs se trouvent 
réduits à des chaines de valeurs DS sauf D 2 qui est composé de deux chaines : 
{(x 2 , 1), (x 2 ,2)} et {(x 2 ,3)}. L’application de l’algorithme BabDs à cet exemple se 
résume à instancier les variables aq, £ 3 , X 4 , X 5 et x§ par les seules chaines qui consti¬ 
tuent leurs domaines respectifs. Pour x 2 , il y a deux distanciations possibles qui 
correspondent aux deux chaines données ci-dessus. 
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Commençons par réduire D 2 à la chaîne {(x 2 ,l), (x 2 ,2)}. A ce stade, tous les 
domaines sont des chaines de valeurs DS. Le rétablissement de l’arc-consistance 
(ligne 2) ne change rien à ce fait d’après la Proposition 8 . Le graphe d’inconsistance 
du problème obtenu est illustré dans la Figure 4.5-(a). Par la suite, l’algorithme 
choisit un élément minimum de chacun des domaines (ligne 1). D’après le Tableau 
4.1 et la Figure 4.5-(a), le choix est unique pour chacune des variables et l’on obtient 
l’instanciation complète (1, 2, 3, 0,1, 3) qui est bien une solution qui satisfait toutes 
les contraintes et qui a un coût égal à 4 unités de temps. 

Enfin, BabDs réduit D 2 à {(x 2 , 3)}. En appliquant l’algorithme d’arc-consistance, 
puis la fonction heuristique h sur les domaines ainsi réduits, on déduit qu’il n’est 
pas possible de trouver une solution meilleure que celle déjà trouvée. En effet, les 
domaines de x 3 et de x§ sont tous les deux réduits à la seule valeur 4 (voir Figure 
4.5-(b)). Ce qui veut dire que la meilleure solution qu’on peut espérer trouver a un 
coût de 5 unités de temps. 



Figure 4.5 - (a) Graphe d’inconsistance du VCSP-CbC associé au job-shop décrit dans la 
Figure 4.2 après réduction de D 2 à la chaine {( x 2 ,1), ( x 2 , 2)} et application d’un algorithme 
d’arc-consistance, (b) Graphe d’inconsistance du même VCSP-CbC après réduction de D 2 
à la chaine {(x2,3)} et application d’un algorithme d’arc-consistance. 


4.5 Orientation du graphe d’inconsistance 

L’orientation du graphe d’inconsistance, utilisée comme référence (voir Définition 
32) pour déterminer les valeurs DS, a un grand impact sur l’occurrence des valeurs 
DS. Les meilleures résultats ont été obtenus lorsqu’on se référé à une orientation qui 
favorise l’occurrence des valeurs DS à des niveaux proches de la racine de l’arbre 
de recherche. Une telle orientation est obtenue dynamiquement en utilisant un algo- 
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rithme glouton (voir Fonction 18). 

Au départ du processus de recherche, l’orientation A est vide. A chaque étape, 
l’algorithme d’orientation ne considère que les arêtes qui relient les valeurs de la 
variable courante à ceux des variables déjà instanciées. Soit a une valeur quelconque 
de la variable courante qui est incompatible avec une valeur b d’une variable déjà 
instanciée. L’arête {a, b} donnera lieu à l’arc (a, b), qui sera ajouté à A. 

Cette stratégie favorise l’occurrence des valeurs DS à des niveaux proches de la 
racine. En effet, rappelons que seulement les arcs sortants d’une valeur sont pris en 
considération lors du calcul des conflits directionnels de cette valeur. En ne considé¬ 
rant que les arêtes qui vont vers les valeurs des variables instanciées, on augmente 
les chances des valeurs des variables traitées à des niveaux proches de la racine à 
être DS les une aux autres. La raison de ceci est que, proche de la racine, il y a peu 
de variables instanciées. 

Le fait de considérer de larges chaînes de valeurs DS à des niveaux peu pro¬ 
fonds de l’arbre de recherche est une stratégie souvent avantageuse car, à ce stade, 
l’heuristique d’ordre des valeurs ne dispose pas d’assez d’informations pour faire des 
choix corrects de valeurs. En prenant des décisions portant sur plusieurs valeurs, on 
peut éviter des erreurs qui, plus elles sont commises tôt, plus elles sont coûteuses en 
temps d’exploration. 


Algorithme 18: Orientation (x, Y, V) : A 

// avec V = (A, D, C, S) 
pour chaque y G X \ Y \ C xy G C faire 
pour chaque a G D x faire 

pour chaque b G D y \ (j) xy (a, b) = _L faire 

|_ A< — A U {(a, 6)}; 

renvoyer(A); 


4.5.1 Algorithme de résolution 

La notion de substituabilité directionnelle peut être utilisée pour améliorer la résolu¬ 
tion de VCSPs-CbC et ceci en l’intégrant à l’algorithme du Branch-and-Bound 
pour obtenir l’algorithme BabDs (voir Fonction 19). Pour avoir une méthode de 
résolution plus efficace, BabDs intègre également une procédure qui permet de 
maintenir l’arc-consistance durant la recherche de solution comme dans l’algorithme 
de résolution de CSP classique MAC [103]. 

BabDs prend comme paramètres le VCSP-CbC à résoudre (A, D, C,S ), (qui 
est supposé être arc-consistant), une orientation du graphe d’inconsistance du pro- 
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bicme A, l’ensemble des variables non encore instanciées Y et la meilleure solution 
déjà trouvée I* et procède comme suit : il sélectionne une variable non encore ins- 
tanciée, calcule la partition de son domaine de valeurs en des chaines de valeurs 
DS et décompose le problème courant en deux sous-problèmes. Le premier sous- 
problème est obtenu en réduisant le domaine de valeurs de la variable courante à 
une chaine de valeurs DS (notée D l k dans le pseudo-code). L’arc-consistance du pro¬ 
blème résultant est alors restauré en utilisant un algorithme d’arc-consistance. Puis, 
un appel récursif est effectué pour considérer les variables restantes. Cet appel per¬ 
met d’obtenir la meilleure solution du sous-problème qui limite les valeurs possibles 
de la variable courante aux éléments de D i k . Ensuite, un processus de restauration 
des domaines de valeurs est effectué et D i k est éliminée du domaine de la variable 
courante. Après restauration de l’arc-consistance du problème résultant, l’algorithme 
effectue un deuxième appel récursif pour considérer les autres chaines de valeurs DS. 

Les deux appels récursifs ne sont effectués que si une fonction heuristique (h) 
indique qu’il est possible d’obtenir des solutions de qualité meilleures que celles déjà 
trouvées. Si l’algorithme réussit à instancier toutes les variables alors il exécute un 
processus polynomial (ligne 1) pour extraire la meilleure solution à partir des chaines 
sélectionnées. 

4.6 Substituabilité directionnelle pour les VCSPs 

Dans cette section, nous utilisons la méthode de décomposition par la substituabilité 
directionnelle pour résoudre les VCSPs binaires en passant par leur duaux, qui sont, 
rappelons le, des VCSP-CbC. 

Le problème dual d’un VCSP est une reformulation du problème qui exprime 
chaque contrainte du problème original comme une variable. Les problèmes duaux 
contiennent uniquement des contraintes unaires et binaires et sont donc des pro¬ 
blèmes bianires. Par conséquent, on peut leur appliquer les algorithmes façonnés 
pour de tels problèmes. 

Le problème dual d’un VCSP est une reformulation de ce dernier qui considère 
chaque contrainte du problème original comme une variable. La contrainte unaire 
associée à une telle variable spécifie et les coûts unaires et les coûts binaires donnés 
par les contraintes unaires et binaires du problème original. 

Les contraintes binaires du problème dual expriment le fait que les variables 
communes à deux contraintes doivent avoir la même valeur. 

Définition 34. Soit V = (. X,T>,C , S) un VCSP bianire. Le dual de V, noté V* est 
défini par le quadruplet (X*,V*,C*, S*) tel que : 

- X* = C ; 

- V* est tel que D* = {t E D x x D y \ a(c) = (x , y)} ; 
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Algorithme 19: BabDs(7 :> , A, Y, /*) 


// avec V = (X, D, C, S) 
si Y = 0 alors 
1 renvoyer(min(D)) 


sinon 


2 


x Select(Y) 

5 X DS-PARTITION( J D æ ,A,P) 
D x ,k <r- Select(< 5 x ) 

-D x ^ 

AC{V) 

si 0 D et h{V) -< 0(7) alors 
I BabDs("P, Y \ {x}, A, /*) 
si 0(7) -< 0(7*) alors 

_ l 1 '^ 1 

Restaurer(R) 

Dx ^ D x \ D x ,k 

AC('P) 

si 0 D et h(V) -< 0(7*) alors 
7 BabDs("P, Y \ {x}, A, I *) 
si 0(1) -< 0(7*) alors 
LJ* <- I 

Restaurer(R) 
renvoyer(7*) 
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0c 

1 

2 

1 

a 

T 

2 

T 

/d 


Tableau 4.2 - La fonction de valuation de l’Exemple 4.6.1 


- S* = S; 

- C* = C* 1 U C* 2 où 

- C* 1 = {((c) ,ip)\ c= ((x, y) , 0) G C} avec <p(t) = <j> x (t {. x) © <f>(t) ; 

- C* 2 = {((ci, c 2 ) , ipi t2 ) | Ci, c 2 G C A cr(ci) fl cr(c 2 ) ^ 0} avec 

( t t r\ i -L si (t I cr(ci) D cr(c 2 )) = (t' | cr(ci) fl a(c 2 )) 

^ l!2 ' ’ | T sinon 

où cr(c) désigne la portée de la contrainte c. 

Comme on peut le constater de la définition précédente, toutes les contraintes 
binaires du dual d’un VCSP binaire sont des contraintes binaires dûres. D’où la 
propriété suivante : 

Propriété 9. Le dual d’un VCSP est un VCSP-CbC. 

Exemple 4.6.1. Soit V est un VCSP composé de trois variables x\, x 2 et x^. 
(Voir Figure 4.6 gauche). Les domaines D t de chaque variable Xi est formé de deux 
valeurs Di = { 1 , 2 }. Ce VCSP est composé de trois contraintes c\ = ((xi,x 2 ) ,(p c ), 
c 2 = ((x 2 , X 3 ) , 0 c) et C 3 = ((xi, X 3 ) , 0 c) dont la fonction de valuations 0 C est définie 
par le Tableau 4.2. 

En se référant à la Définition 34, le daul de V est un problème V* défini comme 
suit : 


- V* est un VCSP composé de trois variables (les contraintes de V) : ci, c 2 et 
C3. (Voir Figure 4.6 droite). 

- Le domaine Dj de chaque variable c 3 est formé de quatre tuples de valeurs 

Dj = {( 1 , 1 ), ( 1 , 2 ), ( 2 , 1 ), ( 2 , 2 )}. 

- Les contraintes binaires dans V deviennent des contraintes unaires dans V* : 

- pour chaque variable Cj, la valeur ( 1 , 1 ) a un coût unaire de a, la valeur ( 2 , 2 ) 
a un coût unaire de 0 et la valeur (2,1) a un coût unaire de T et elle sera 
filtrée par application d’un algorithme d’arc-consistance. 

- des contraintes binaires dans V doivent être ajoutées pour interdire de choisir 
deux valeurs appartenant à la même variable dans V. Ces contraintes binaires 
sont dures comme elles interdisent des solutions impossibles pour le problème 
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primai. 

Comme on peut le voir sur la Figure 4.6 V* est un VCSP-CbC. 



X3 


C3 


Figure 4.6 - Dual d’un VCSP. 

Il en découle de la Propriété 9 que l’on peut appliquer dual d’un VCSP binaire, 
les algorithmes façonnés pour les VCSP-CbC. Nous proposons donc appliquer, au 
dual d’un VCSP binaire, le schéma de décomposition basé sur la substituabilité di¬ 
rectionnelle. 


Propriété 10. Les VCSPs sont solubles par le schéma de décomposition par la 
substituabilité directionnelle. 


Preuve. En se référant à la Propriété 9, le problème dual d’un VCSP binaire est un 
VCSP-CbC. En se référant à la Définition 34 la transformation d’un problème pimal 
au dual correspondant se fait en temps polynomial car il faut créer On 2 variables 
ayant des domaines de taille de l’ordre de O(d 2 ). Enfin, le nombre de contraintes 
dans le dual est de l’ordre de ()(n :^ ). Donc, chaque solution du VCSP-CbC peut 
être transformée en une solution du VCSP original en temps polynomial. D’autre 
part, par la Propriétés 7 un VCSP-CbC arc-consistant dont chacun des domaines 
de valeurs est une chaîne non vide de valeurs dir ectio nn ellement substituables par 
rapport à une orientation A est résolu en temps polynomial. Selon la Propriété 8, le 
passage d’un problème à l’une de ses réductions préserve la substituabilité direction¬ 
nelle. L’application de l’Algorithme 19 sur le dual d’un VCSP permet d’exploiter la 
substituabilité directionnelle pour le décomposer en plusieurs sous-problcmes chacun 
est résolu en temps polynomial. 


□ 


74 









Chapitre IV : Décomposition par substituabilité directionnelle 


4.7 Résultats expérimentaux 

Pour montrer l’apport du schéma de décomposition proposé, nous avons mené une 
étude expérimentale impliquant des problèmes d’ordonnancement. Les problèmes 
test utilisés sont ceux référencés par js-taillard-15 disponibles sur [80]. Ce sont 
dix instances de job-sliops comportant 15 jobs et 15 machines, (donc 225 tâches), 
générés selon le modèle décrit dans [112]. Pour ces instances de job-shop, la discré¬ 
tisation du temps maximum impartie pour exécuter toutes les tâches donne plus de 
mille valeurs. Pour obtenir des solutions en des temps d’exécution raisonnables pour 
des VCSPs-CbC impliquant de si large domaines de valeurs, on a élargi davantage 
les domaines de valeurs pour obtenir des instances plus faciles. Ainsi, on a multiplié 
les tailles des domaines originaux par un facteur de 105%, 104% puis 103% pour 
obtenir des instances qui vont des plus faciles aux plus difficiles. 

L’algorithme d’arc-consistance employé est AC-3 [87]. L’heuristique de sélection 
de variables utilisée est min-doraain/wdeg [9], pour la sélection de valeurs, on a 
utilisé deux heuristiques : D’abord, celle qui privilégie les valeurs qui ont le moins de 
conflits [43] (le BAB sera noté dans ce cas par MC). Ensuite, celle qui privilège les 
valeurs les plus faibles (le BAB sera noté dans ce cas par MV). La fonction unaire 
4> x (a), qui est indispensable pour définir la substituabilité directionnelle, est déduite 
de <f>p(T) de manière analogue à celle spécifiée par l’équation (4.5). 

Les critères d’évaluation des performances sont la taille de l’arbre de recherche 
explorée, le temps de calcul et la qualité de la meilleure solution trouvée. Les deux 
algorithmes ont été implémentés en C++ et exécutés sur un PC à 2 GHZ de fré¬ 
quence et une mémoire vive de 4GB. 

Les résultats obtenus (voir Tableau 4.3, Tableau 4.4 et Tableau 4.5) montrent 
clairement que MC+DS et MV+DS ont trouvé des solutions pour plus d’instances 
que MC ou MV. Pour les instances auxquelles les deux algorithmes ont réussit à 
trouver des solutions, on constate que dans la plupart des cas, ou bien que la so¬ 
lution trouvée par BabDs a un coût meilleur que celle trouvée par BAB, ou bien 
que BabDs trouve une solution ayant le même coût mais en moins de temps. Sur les 
trente instances considérées, il y a uniquement six instances (MC-104-4, MC-104-7, 
MC-103-7, MV-105-6 et MV-104-9) sur lesquelles, BAB a été meilleur que BabDs. 
Enfin, la supériorité de BabDs ne semble pas dépendre de la difficulté des instances, 
ni de l’heuristique employée puisque BabDs est globalement meilleur sur les deux 
niveaux de difficulté. 


4.8 Conclusion 

Dans ce chapitre, nous avons proposé les extensions suivantes au travail présenté 
dans [90] : 
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temps coût temps coût 


instances 

MC 

MC+DS 

MC 

MC+DS 

MV 

MV+DS 

MV 

MV+DS 

105-0 

7065 

651 

1308 

1308 

- 

225 

- 

1290 

105-1 

90 

3847 

1326 

1325 

- 

1573 

- 

1304 

105-2 

11521 

748 

1295 

1295 

- 

1173 

- 

1285 

105-3 

- 

3743 

- 

1239 

- 

6193 

- 

1234 

105-4 

- 

5503 

- 

1297 

4491 

12943 

1291 

1285 

105-5 

- 

2532 

- 

1308 

- 

2249 

- 

1309 

105-6 

105 

6056 

1297 

1296 

2292 

9559 

1284 

1286 

105-7 

- 

376 

- 

1282 

- 

523 

- 

1282 

105-8 

3406 

74 

1353 

1353 

- 

2541 

- 

1343 

105-9 

3051 

4681 

1334 

1333 

739 

13850 

1334 

1316 


Tableau 4.3 - Résultats obtenus par MC, MC+DS, MV et MV+DS sur les instances js- 
taillard-15-105. Le temps CPU en secondes et le coût de la meilleure solution trouvée sont 
raportés. Pour le temps CPU, on indique les temps auxquels ont été trouvées les meilleures 
solutions. Le temps d’exécution a été limité à 4 heures par instance. 


temps coût temps coût 


instances 

MC 

MC+DS 

MC 

MC+DS 

MV 

MV+DS 

MV 

MV+DS 

104-0 

- 

4526 

- 

1296 

12157 

3366 

1292 

1279 

104-1 

- 

460 

- 

1313 

- 

12147 

- 

1302 

104-2 

- 

88 

- 

1283 

- 

13812 

- 

1269 

104-3 

646 

2716 

1228 

1227 

- 

11080 

- 

1228 

104-4 

1500 

12942 

1285 

1285 

- 

10504 

- 

1281 

104-5 

- 

643 

- 

1297 

- 

- 

- 

- 

104-6 

- 

1136 

- 

1284 

- 

766 

- 

1282 

104-7 

2606 

14028 

1270 

1270 

- 

- 

- 

- 

104-8 

61 

539 

1341 

1340 

- 

4137 

- 

1340 

104-9 

99 

79 

1321 

1321 

739 

- 

1318 

- 


Tableau 4.4 - Résultats obtenus par MC, MC+DS, MV et MV+DS sur les instances js- 
taillard-15-104. Le temps CPU en secondes et le coût de la meilleure solution trouvée sont 
raportés. Pour le temps CPU, on indique les temps auxquels ont été trouvées les meilleures 
solutions. Le temps d’exécution a été limité à 4 heures par instance. 
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temps 


coût 

temps 

coût 

instances 

MC MC+DS 

MC 

MC+DS 

MV MV+DS 

MV MV+DS 

103-0 

122 

- 

1284 

- 

- 

103-1 

885 

- 

1301 

3813 

1284 

103-2 

423 

- 

1271 

389 

1262 

103-3 

- 

- 

- 

9734 

1213 

103-4 

- 

- 

- 

218 

1270 

103-5 

- 

- 

- 

- 

- 

103-6 

1274 

- 

1273 

- 

- 

103-7 

2513 

1259 

- 

- 

- 

103-8 

- 

- 

- 

132 

1327 

103-9 

454 

- 

1309 

8926 

1287 


Tableau 4.5 - Résultats obtenus par MC, MC+DS, MV et MV+DS sur les instances js- 
taillard-15-103. Le temps CPU en secondes et le coût de la meilleure solution trouvée sont 
raportés. Pour le temps CPU, on indique les temps auxquels ont été trouvées les meilleures 
solutions. Le temps d’exécution a été limité à 4 heures par instance. 

Nous avons généralisé la notion de substituabilité directionnelle en considérant, 
comme référence, une orientation du graphe d’inconsistance au lieu d’un ordre sur les 
variables. Pour garantir l’optimalité de la solution lors de la résolution des Problèmes 
de Satisfaction de Contraintes Valuées à Contraintes binaires Dures (VCSPs-CbC), 
nous avons introduit des conditions supplémentaires de substituabilité. 

Enfin, nous avons exploité la notion de substituabilité directionnelle, proposée pour 
les VCSPs à Contraintes binaires Dures, pour résoudre les VCSPs 

Les résultats de simulation sur plusieurs VCSPs-CbC qui codent des problèmes 
d’ordonnancement du type job-shop ont montré qu’une variante de l’algorithme du 
Branch-and-Bound exploitant la substituabilité directionnelle est souvent plus 
efficace que l’algorithme original. 

Une extension possible de ce travail consiste à proposer une formulation encore 
plus générale que la substituabilité directionnelle afin de pouvoir l’appliquer aux 
VCSPs en profitant et en proposant des classes polynomiales plus importantes lors 
de la décomposition. 


77 









Chapitre 5 


Une méthode de décomposition 
modulaire pour les VCSPs binaires 

Sommaire 

5.1 Introduction . 79 

5.2 Définitions et Notations .80 

5.3 Les fonctions modulaires binaires.81 

5.3.1 Définitions et propriétés. 81 

5.3.2 Décomposition de fonctions modulaires binaires. 83 

5.3.3 Décomposition en fonctions modulaires binaires. 85 

5.4 Un schéma de décomposition pour les VCSPs.89 

5.4.1 VCSP avec fonctions modulaires binaires . 90 

5.4.2 La décomposition des domaines. 92 

5.4.3 La solution algorithmique. 93 

5.5 Résultats expérimentaux.94 

5.6 Conclusion .101 


78 
















Chapitre V : Une méthode de décomposition modulaire pour les VCSPs 


Résumé 

Dans ce chapitre nous étudions la résolution de certains problèmes combinatoires. Ce sont des 
problèmes qui peuvent être formulés en tant que problèmes de satisfaction de contraintes valuée 
(VCSPs). Dans le cadre des VCSPs, les contraintes sont définies par l'utilisation de fonctions de 
valuations afin de refléter plusieurs degrés de cohérence. 

Malgré la complexité des VCSPs dans le cas général (ce sont des problèmes NP-difficiles), 
plusieurs sous-problème polynomiaux peuvent être obtenues en se limitant à des fonctions de va¬ 
luation avec des propriétés mathématiques spécifiques. C’est le cas, par exemple, des VCSPs qui 
utilisent des fonctions de valuations sous-modulaires [14]. 

Dans ce chapitre, nous proposons un schéma de décomposition pour les problèmes de VCSPs 
binaires. Notre méthode utilise les fonctions de valuation modulaires, même lorsque le problème 
étudié ne se limite pas à cette classe de fonctions. Comparées aux fonctions sous-modulaires, 
les fonctions modulaires sont moins fréquentes, mais, en compensation, elles sont plus faciles à 
traiter. Le schéma de décomposition de problème proposé est intégré à un algorithme de résolution 
basé sur le retour-arrière chronologique et consiste à décomposer le problème à résoudre en un 
ensemble de sous-problèmes modulaires, dont la résolution est polynomiale. 


5.1 Introduction 

Ce Chapitre a fait l’objet des publications suivantes [56, 58]. 

La complexité algorithmique de la recherche de la solution optimale d’un VCSP 
a été étudiée dans plusieurs travaux de recherche. Plusieurs classes polynomiales de 
VCSPs ont été identifiées et résolues [24, 70, 17]. Les VCSPs impliquant les fonctions 
de valuations binaires sous-modulaires est une de ces classes polynomiales. En effet, 
en exprimant ces VCSPs comme des problèmes de calcul de la coupe minimale d’un 
graphe pondéré, il est possible de les résoudre en 0(n 3 cl 3 ) étapes où n est le nombre 
de variables et d est la taille du domaine de valeur le plus large [14]. 

Néanmoins, les VCSPs résultant de cas réels ne sont pas limités à des fonc¬ 
tions sous-modulaires. Dans de tels cas, pouvons nous procéder de manière plus 
efficace qu’une recherche exhaustive tout en garantissant l’optimalité de la solution 
proposée ? Y a-t-il un moyen d’exploiter la sous-modularité dans un contexte plus 
restreint ? Ce chapitre contribue à fournir des réponses positives à ces questions. 
Motivé par les résultats obtenus sur les CSPs [91], nous présentons une méthode de 
décomposition de domaine pour résoudre les VCSPs binaires. Elle s’appuie sur la 
classe des fonctions modulaires [93]. Le processus de décomposition produit des sous 
problèmes dont les fonctions de valuations sont exclusivement modulaires. 

Pour les VCSPs binaires avec de telles fonctions de valuations, nous proposons 
un algorithme d’identification de complexité 0(ed 2 ) et un algorithme de résolution 
de complexité 0(ed). 
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Le chapitre est organisé comme suit : la section suivante présente quelques défini¬ 
tions et notations. La Section 5.3 est consacrée aux fonctions binaires modulaires et 
leur décomposition. Dans la Section 5.4, nous présentons le schéma de décomposition 
et des solutions algorithmiques pour les VCSPs binaires. Des résultats expérimen¬ 
taux sont présentés dans la Section 5.5. Finalement, la Section 5.6 est une brève 
conclusion. 

5.2 Définitions et Notations 

Rappelons qu’une structure de valuation est un triplet S = (E,©, tQ tel que : 

- E est l’ensemble de valuations contenant un élément minimal _L et un élément 
maximal T ; 

- 7 < est un ordre total sur E ; 

- © est un opérateur commutatif, associatif et monotone. 

Dans la suite de ce manuscrit, nous ne traiterons que des structures de valuations 
propres et strictement monotone (voir Chapitre 1). 

Exemple 5.2.1. Considérons la structure de valuation donnée par S'® = ([0,1], 0 , < 
), ou l’opérateur d’agrégation 0 est défini de la sorte que a 0 /3 = a + (3 — aj3. 
Notons que a © /3 est dans [0,1], pour toutes a,/3 dans [0,1]. En plus, pour toutes 
a, (3, 7 G [0,1] tel que a < (3 e t 7 ^ 1, nous avons a(l — 7 ) + 7 < (3(1 — 7 ) + 7 , et 
donc a + 7 — « 7</3 + 7 — /hy, qui est équivalent àa® 7 <, 8 ® 7 . Ainsi, S 1 ® est 
strictement monotone. 

Avec la structure de valuation S 1 ®, si a, (3 G [0,1], a A /3, alors la différence entre 
[3 et a, notée par /3 0 a, existe, elle est dans [0,1] et elle est définie par 

a f j—^ si a 7 ^ 1 

i? 0 a= r 

\ 1 smon 

Il est facile de vérifier que a 0 (/3 0 a) = (3. D’ofi S 1 ® est une structure de valuation 
propre et strictement monotone. 

Les structures de valuations propres et strictement monotone vérifient les deux 
propriétés suivantes, qui sont prouvés dans [18], (Lennne 7 et Théorème 38). 

Propriété 11. Soit S = (E, ©, A) une structure de valuation et strictement mono¬ 
tone. Alors pour tout a, /3 ,7 G E tel que 7 A (3, on a (/ 3 © 7 ) 7 ! (3 et (a® 7 )©(/?© 7 ) = 
a © j3. 

Propriété 12. Soit S = ( E , ffi, 7 /) une structure de valuation propre et strictement 
monotone. Alors pour tout a, /3 ,7 G E tel que 7 7 ! /3, on a (affi/3) ©7 = a© (^© 7 ). 

Par ailleurs et comme nous l’avons mentionné précédemment, nous nous intéres¬ 
sons aux algorithmes complets qui fonctionnent par recherche en profondeur d’abord 
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(DFS) dans un arbre de recherche jusqu’à ce qu’une solution optimale soit trouvée 
ou qu’il soit prouvé qu’aucune solution n’existe. 

A chaque nœud de l’arbre de recherche, les algorithmes de résolution du type, 
DFS, réduisent le problème courant en affectant une valeur à une variable et en en¬ 
levant les valeurs qui sont désormais incohérentes avec la solution partielle courante 
ou ceux menant à des solutions sous-optimales. 

Pour des raisons de simplicité, nous supposons que les domaines de valeurs sont 
des sous-ensembles de nombres entiers. Chaque variable x devrait être affectée par 
une valeur de son domaine. L 1 arité d’une contrainte valuée est la taille de sa por¬ 
tée. Dans ce chapitre, nous traitons les VCSPs binaires qui sont des VCSPs avec 
des contraintes unaires et binaires seulement. Etant donné un ordre —> sur les va¬ 
riables du problème, nous définissons le sous-ensemble de contraintes binaires dont 
les portées sont des séquences ordonnées de variables comme suit : 

c = {((a, y) , <f>zy) eC \ x ^y} 

Chaque fonction de valuation </> : D x D' —>• E qui sera examinée lors de l’exposi¬ 
tion de notre approche est supposée être arc-consistante. En d’autres termes, chaque 
valeur doit avoir au moins un support, ce qui se traduit par la formule suivante : 

Vu G D , 3v r G D', tel que (j>(v,v r ) -< T (5.1) 

Ceci ne présente en aucun cas une limitation de notre approche, puisque toute 
fonction de valuation non arc-consistante peut être transformée en une autre fonc¬ 
tion qui soit arc-consistante simplement en supprimant les valeurs qui ne remplissent 
pas la condition (5.1). Dans le cadre VCSP, appliquer une telle transformation sur 
chaque fonction de valuation est une tâche polynomiale standard. Cette tâche se 
traduit par la réduction des domaines de valeurs sans changer l’ensemble de toutes 
les solutions possibles au problème. En effet, les valeurs éliminées ne peuvent parti¬ 
ciper à aucune solution optimale. Le problème qui résulte des suppressions décrites 
ci-dessus est dit arc-consistant. 


5.3 Les fonctions modulaires binaires 

Notre schéma de décomposition s’appuie essentiellement sur la classe des fonctions 
modulaires. Dans cette section, nous définissons la fonction modulaire binaire, met¬ 
tons en évidence certaines de ses propriétés et montrons qu’elle peut être exprimée 
comme la somme de deux fonctions de valuation unaires. 

5.3.1 Définitions et propriétés 

Considérons l’ensemble de toutes les fonctions binaires (f> : D x D* —ï E, oti D et 
D' sont deux ensembles finis totalement ordonnés et E est l’ensemble de valuations 
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possibles d’une structure de valuation propre et strictement monotone (E, ©, ^). 

Définition 35. Soit S = (E, ©, P) une structure de valuations propre. Une fonction 
binaire <f> : D x D' —»■ E est modulaire si, pour toutes u, v G D, u', v' G D', on a 

0(w, u') © 0(u, v') = 0(w, v') © 0(u, u') (5.2) 

Notons que si D est totalement ordonné alors toutes les fonctions unaires D —* E 
sont modulaires. 

Désignons par A4 l’ensemble de toutes les fonctions unaires et binaires modu¬ 
laires. 

Soit 0 et <f> T deux fonctions binaires telles que pour toutes v G D et v' G D', on 
a cf> T (v,v') = (f>(v',v). Nous pouvons facilement déduire de (5.2) que 0 est dans A4 
si et seulement si 0 T est dans A4. 


Propriété 13. Pour toute 0 dans AA et toutes u',v' dans D', on a 

3 ueD, 0(m, u') -< 0(w, v') =r- Vu G D, 0(n, u') P 0(w, v') (5.3) 


Preuve. Nous supposons l’inverse et procédons pour obtenir une contradiction. 
Supposons donc qu’il existe u,v G D et u', v' G D' tel que 0(w, u') -< 0(w, v') et 
4>(v,v') -< <f>{v,u'). Le fait que © est strictement monotone implique que 0(w, «') © 
4>(v,v') -< 4>(u,v') © 0(n,n / ), ce qui contredit le fait que 0 est modulaire. 


□ 

La Propriété 13 suggère que la comparaison entre 0(w, u ') et 0(w, v') pour toutes 
v!,v' G D' ne dépend pas de u. En outre, puisque A est un ordre total sur E , il 
induit un préordre total sur les éléments de D ', qui sera noté par Ce préordre 
est défini comme suit 


u'&def 0(m, u') A v') (5.4) 

où m = min D. 

Propriété 14. Pour toute 0 G Al, si 0 et 0 T sont deux fonctions arc-consistantes 
alors <f>{u,u') -< T pour tout (u,u') G D x D'. 


Preuve. Nous supposons l’inverse, et procédons pour obtenir une contradiction. 
Soit u,u' G D x D' tel que 0(-u,-u / ) = T. Puisque 0 et 0 T sont arc-consistantes 
alors de (5.1) nous déduisons qu’il doit exister une valeur » G D et une valeur 
v' G D' tel que <f>{u,v') -< T et <f>{v,u') -< T. Nous obtenons 0(u,n') © 0(n,u / ) -< 
0('U,'U / ) © 0(n,n / ) = T, qui est une contradiction avec le fait que 0 est modulaire. 
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□ 

Propriété 15. Pour toute 0 G A4, si 0 et f> T sont deux fonctions arc-consistantes, 
alors pour toutes u,v G D et u',v' G D' tel que 4>{u,v') P f>(u,u'), on a 

(f>(u, u') © 0(w, v') = f(v,u') © 0(u, v') 


Preuve. 0 est dans A4 alors, pour toutes u,v G D, u',v' G D' , on a 


f(u,u') © = f{u,v') © f{v,u') (5.5) 

D’après la Propriété 14, on a f{v,v') -< T. On obtient donc 


0(w, U') = (0(w, v') © 001, U ')) © 0(v, V ') 

D’autre part, on doit avoir f(v,v') P car sinon, par la stricte monotonicité 

de ©, (5.5) ne pourra pas être vraie. Il s’en suit, par la Propriété 12 et le fait que © 
est commutative, que 


0(w, U') = (001, U ') © 001, V ')) © 0(w, V') 

En plus, puisque <f{u,v') P <f(u,u'), on obtient 

4>(u,u') © 0(w, v') = (4>(v,u') © 4>(v,v') © f{u,v')) © 4>(u,v') 

Encore d’après la Propriété 14, on a f>(u,v') -< T. Il en résulte d’après la Propriété 
12 que 

0(u, u ') © 0(u, v') = 0(w, u') © 0(w, v') 


□ 


5.3.2 Décomposition de fonctions modulaires binaires 

Dans ce qui suit, nous montrons qu’une fonction modulaire binaire peut être ex¬ 
primée comme la somme de deux fonctions unaires. Par conséquent, dans ce para¬ 
graphe, 0 désignera une fonction modulaire arc-consistante. Désignons par m* l’un 
des éléments minimaux de D' selon On a donc 

m* = arg min <f(y,v') (5-6) 

v'&D' 

ofi v est un élément quelconque de D. Puis considérons les fonctions unaires 0i et 
02 définies respectivement sur D et D' avec des valuations dans E comme suit 


0i (v) = <t>(v,m*) (5.7) 

020 /) = 4>(u,v') © 4>(u,m*) (5.8) 

où u est un élément quelconque de D. Notons que 0 2 (v 1 ) existe pour toute v' G D' 
comme par (5.6), nous avons © <f{u,v') pour toute u E D. 

En plus, selon la Propriété 15, 0 2 ('î/) ne dépend pas de u. 
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Propriété 16. Soit 0 : D x D' —y E est une fonction binaire alors 0 est modulaire 
si et seulement si = 0i(w) © 0 2 (u , ) J pour tout (u,u') E D x D'. 


Preuve. 0 est modulaire et arc-consistante, alors d’après la Propriété 14, toutes les 
valuations obtenues par 0 sont plus petites que T. De plus, pour toutes u, v E D , u' 
E D ', on a 

0 (w, u') © 0(u, m*) = 0 (w, m*) © <f>(v, u') 

0(w,«') = 0i( m) © (0(u, w') © 0(u,m*)) 

4>(u,u') = 0i(m)©0 2 (m') (5.9) 

La seconde égalité est obtenue par (5.7) et la Propriété 12, et la dernière est obtenue 
par (5.8). 

Dans le sens inverse, si 0(w, w') = 0i(w) © 0 2 (u'), pour toutes (u,u') E D x D ', 
alors nous avons 

0(M,M') © 0(u,'l/) = 01 («) © 02 («') © 01 (ü) © 02 K) 

= 0 (w, m*) © ( 0 (v, m') © 0 (w, m*)) © 0 (n, m*) 

© (0(m,u') © 4>(u,m*)) (5.10) 

= 0 (w,u') © 4>(v,u') 

La seconde égalité est obtenue par (5.7) et (5.8). La dernière égalité, qui est obtenue 
par la Propriété 11, implique que 0 est modulaire. 

□ 

Ci-dessous est présenté le pseudo-code d’un algorithme qui prend comme para¬ 
mètres un couple de D x D ', une fonction modulaire binaire 0 : D x D' —> E et un 
élément u de D et calcule les fonctions unaires 0i et 0 2 selon (5.7) et (5.8). 


Algorithme 20: Decompose(0, D, D') 

m* <— argmhvgxr 4>(v,v') 
pour u E D faire 
|_ 01 (u) E- 0(tt, m*) 
pour u' E D 1 faire 

02 ( 11 ') E- f(v, u ') 0 (f>(v, m*) 

RENVOYER 0i, 0 2 


Propriété 17. Soit <f E M. alors l’algorithme Décomposé calcule deux fonctions 
unaires vérifiant la Propriété 16 en 0(max(|D|, |D'|)) étapes. 


Preuve. Ceci peut être facilement déduit du pseudo-code de Décomposé. 
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: o 0.5 ' 

I ® © j 


® © 

0 0.5 | 


(a) 


(b) 


Figure 5.1 - Illustration de l’algorithme Décomposé 


0 

0 

1 

0 

0 

0.5 

1 

0.5 

0.75 


Tableau 5.1 - Définition de la fonction binaire fi 


□ 

Propriété 18. Une fonction binaire peut être identifiée comme modulaire en 0(|D|.|.D / |) 
étapes. 


Preuve. D’après la Propriété 16, il suffit d’appliquer l’algorithme Décomposé à 
la fonction en question, qui calcule fi 1 et 0 2 en 0(max(\D\, \D '\)) étapes, puis de 
vérifier les |D|.|.D , | équations données par la Propriété 16. 


□ 

Exemple 5.3.1. Considérons encore la structure de valuation S ® = ([0,1], 0 , <) 
introduite dans l’Exemple 5.2.1. Soit (j) : {0, l} 2 —y [0,1] la fonction binaire définie 
par le Tableau 5.1 (voir la Figure 5.1(a)). fi est modulaire puisque 0(0, 0) <8)0(1,1) = 
0(0,1) 0 0(1, 0) = 0.75. L’algorithme Décomposé calcule m qui est égale à 0 et m! 
qui est aussi égale à 0. Puis, il calcule la fonction unaire fi\. Selon (5.7), fi\ est tel 
que 0i(0) = 0(0,0) = 0 et 0i(l) = 0(1,0) = 0.5. Ensuite, l’algorithme DECOMPOSE 
calcule la fonction unaire 0 2 selon (5.8) pour obtenir 0 2 (O) = 0(0,0) 0 0(0,0) = 0 
et 0 2 (1) = 0(0,1) 0 0(0, 0) = 0.5 (voir la Figure 5.1(b)). 


5.3.3 Décomposition en fonctions modulaires binaires 

Dans ce paragraphe, nous montrons que, étant donné une fonction binaire fi : D x 
D' —y E qui n’est pas nécessairement modulaire, il est possible de partitionner D 
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en de sous-ensembles de sorte qu’en limitant le premier argument de 0 à un de ces 
sous-ensembles, nous obtenons une famille de fonctions binaires modulaires. 

Définition 36. Soit 0 : D x D' —* E une fonction binaire et soit u et v deux 
éléments de D. Nous disons que u et v sont modulaires par rapport à 0 (notation 
u v) si et seulement si la restriction de 0 à {u, v } x D' est modulaire, c ’est-à-dire, 

U V ^def Vu',v'eD', 0(«, u') © <f>(v, v') = 0(lt, v') © 0(f, u') 

Notons que 


0 G Ai u v, pour tout u,veD (5.11) 

Propriété 19. est une relation d’équivalence sur D. 

Preuve. La relation ^ rp est trivialement réflexive et symétrique. La transitivité de 
est obtenue comme suit. Soit u,v,w G D sont tel que u v et v w. Ceci 
implique que, pour tout u',v' G D' , on a 

0(w, u') © 0(u, v') = 4>{u,v') © 4>{v,u') (5.12) 

<p(v,u') © <j)(w,v') = <f>(v,v') © cp(w,u') (5.13) 

Nous distinguons les quatre cas suivants : 

- 0(ryu') -< T A 4>(v,v') © T : 

Ce qui implique que 

0(u, u') © 0(u, u') = 0(y, v') © 0(u, v') = _L (5.14) 

En agrégeant les membres gauches et les membres droites de (5.12) et (5.13), 
on obtient 

0 (w, u') © 0(u, v') © 0(u, U ') © 0(u>, V ') = 0(w, ?/) © 0(v, «') © 0(u, A) © 0(lü, «') 

En soustrayant 0(n, «'), et par la suite 0(w, u'), des deux membres de la dernière 
égalité, et puis, en utilisant la Propriété 12 et (5.14), on obtient 

0 (w,w') © 0 (iu,î/) = 4>(u,v') © (p(w,u') 
ce qui signifie que u w. 

- 0(ryu') ©TA 4>(v,v r ) = T : 

Pour avoir (5.12) et (5.13), nous devons avoir <f>{u,v') = 0(w >,v r ) = T. Il en 
résulte que 


<j)(u,u') © 4>(vj,v') = 4>(u,v') © 4>(vj,u') 
ce qui signifie que u 
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- 4>(v,u') = T A 0(u,u') -< T : 

Ce cas est similaire au cas précédant. 

- (j)(v,u') = T A 4>{v,v') = T : 

Comme 0 est arc-consistante, il doit exister w' G D 'tel que 0(u, w') -< T. Cela 
implique que 4>{u,u') = <f)(w,u') = T, car sinon on aura 

0(w, u') © (j)(v, w') -< 0(w, w') © (j)(v, u') 

ou alors 

cj)(v,w') © -< 0(n,w') © <f>(w,w') 

Ce qui contredit le fait que u v et v w. Il en résulte que 

(j>(u,u r ) © 4>(w,v') = 4>(u,v') © <f){w,u') 
cela signifie que u w. 


□ 

L’algorithme EquivClass, donné ci-dessous, prend une fonction binaire 0 : D x 
D' —ï E et un élément v de D comme arguments et calcule la classe d’équivalence 
de v, (dénotée par v dans le pseudo-code), conformément à la Définition 36. 


Algorithme 21: EquivClass(0, D, D', v) 

(0i,02) g- Décomposé(0, D,D') 
v 0 

pour u € D faire 

equiv true 
pour u' G D' faire 

si (j}(u,u') / 0i (u) © 02(rO alors 
equiv t— false 
break 

si equiv alors 

v <— v U {«} 

RENVOYER V 


Propriété 20. Pour une fonction binaire 0 : D x D’ —^ E et une valeur v G D, 
Valgorithme EquivClass calcule la classe d’équivalence de v selon la définition 36 
dans 0(\D\.\D'\) étapes. 


Preuve. Les fonctions unaires 0i et 0 2 sont obtenues à la suite d’un appelle de 
Décomposé avec les paramètres 0 et v, nous avons alors 
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0iO) = 00 ,m* v ) 

0200 = 00 ) u ') © 00 ) m*) 
tel que m* = argmnveD'O, O). 

D’autre part, on peut déduire du pseudo-code EquivClass qu’un élément u est 
inséré dans v si et seulement si 00 , 0 ) = 0iO) ©02 00 P our tout O £ D'. D’abord, 
nous vérifions que v lui même est dans v. Ceci est vrai puisque, pour tout O G D', 
on a 


0iO)©0aO / ) = 00) m v) © (00 ) v> ) © 00) m v)) 

= 00, O) (5.15) 

où la seconde égalité est obtenue par la Propriété 15 et la troisième par la Propriété 
11. 


Nous prouvons à présent que si u v alors u G u. Puisque u u, alors, pour 
tout u' G D' , on a 


00 , u ') © 

00, «') 
00,0) 
00,0) 


(j){u,m* v ) © 00,0) 

(010) ® 00) O)) © 00, m*) 
0iO) © (00,0) © 4>(v, ml)) 
010 ) © 0200 


L’arc-consistance de (p implique que (p(v, m') -< T. Donc, 0(u, m')Q(p(v, m') = _L. 
La seconde égalité est obtenue en utilisant (5.7). La troisième égalité résulte de la 
Propriété 12 et la quatrième est déduite de 1a, Propriété 15. La cinquième égalité est 
dérivée de (5.7) et (5.8). Elle implique que u G v. 


Pour le deuxième sens, puisque u,v G v, alors, pour tout u' G D ', nous avons 
( p(u,u') = <pi(u) ©02(0) et 0(u,O) = 0i(u) ©02(0). Par (5.7) et (5.8) et pour tout 
0,0 G D’ , on obtient ce qui suit 


00,0) © 00,0) = 0iO) © 0 2 (O) © 0iO) © 0 2 (O) 

= 00,0) © 02(O) © 0iO) 

= 00,0) © (00,0) © 00, m*)) © 00, m*) 

= 00 , 0 ) © 00 , 0 ) 

(5.16) 

La seconde égalité est dérivée de (5.8), 1a, troisième est obtenue de la Propriété 15. 
Ensuite, (5.7) est utilisée pour obtenir la quatrième égalité. La dernière égalité, qui 
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implique que u v, est obtenue de la Propriété 11. 

Enfin, nous pouvons facilement voir à partir du pseudo-code de l’algorithme 
EquivClass qu’il nécessite 0(\D\.\D'\) étapes. 


□ 

Propriété 21. Soit 0 est une fonction et soit D i,..., Dk la partition de D induite 
par Alors la restriction de 0 à D s x D est modulaire pour s = 1, ..., k. 


Preuve. La Propriété résulte immédiatement de la Définition 36 et le fait que chaque 
D s est une classe d’équivalence de 


□ 

Exemple 5.3.2. Considérons encore la structure de valuation S ® = ([0,1], ®,<) 
introduite dans l’Exemple 5.2.1. Soit 0 : {A, C, G, T} 2 —> [0,1] la fonction binaire 
définie comme suit : 

0 ( „, „') = (» si 6 U °} 2 u (c, r} 2 

' I 0.5 smon 

En constatant que <j)(A, u') = u ') et 0(C, u ') = 0(T, u ') pour tout u' G {A, C,G,T}, 

nous déduisons que 


0(4, u ') ® 0(G, v') = cj)(A, v') ® 0(G, u') 
et 

0(C, u') ® 0(T, u') = 0(C, u') ® 0(T, u') 

pour tout u',v' G {A, C, G, T}. D’où, les deux équivalences suivantes : A G et 
C T. Les deux fonctions obtenues par la restriction du premier argument de 0, 
à {A, G} puis {C, T} sont donc modulaires. 

5.4 Un schéma de décomposition pour les VCSPs 

Dans ce qui suit, nous désignons par VCSP(U) la classe de VCSPs binaires n’impli¬ 
quant que des fonctions de valuations dans l’ensemble de fonction C. 

Dénotons par U l’ensemble de toutes les fonctions unaires D —> E. 

Propriété 22. La complexité temporelle de VCSP{U ) est 0(nd), où n est le nombre 
de variables et d est la taille du domaine le plus large du VCSP. 


Preuve. Un algorithme qui résout des instances de VCSP (Z//) procède en exploitant 
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la consistance de nœud, qui garantit qu’il n’y a aucun domaine vide. Il retourne un 
n-tuples t qui vérifie : 


11 x = arg min f x (v ), x E X 

veD x 

Ceci peut être effectué en 0(nd ) étapes. De plus, on a (f>-p(t) = © æe x 00 I x) qui 
est minimal puisque l’opérateur d’agrégation © est monotone. 


□ 

Dans ce qui suit, nous montrons comment réduire une instance de VCSP(Al) à 
une instance de VCSP(W). Ensuite, nous décrivons une méthode pour décomposer les 
instances de VCSPs binaires quelconques en ensembles d’instances dans VCSP(Al). 

5.4.1 VCSP avec fonctions modulaires binaires 

Soit c xy est une contrainte binaire dont la fonction de valuation <p xy est dans Ai. 
Alors, d’après la Section 5.3, c xy peut être exprimée comme la conjonction de deux 
contraintes unaires, dont les portées sont x et y et dont les fonctions de valuations, 
notées (p x et ip y , sont définies conformément à (5.7) et (5.8). De plus 


f X y Px © Py (5.17) 

Définition 37. Soit V = (X,V,C, S) une instance de VCSP(Xt). La réduction 
unaire de V se référant à l’ordre total —y sur X est le VCSP V = (X, V, C', S) tel 
que 

C' = {(0, (f' x ) | x e X et <j/ x = 4> x © 0 ip x ® 0 ç p y } 

P-xy £ ^ Cyx £ C 

Notons que, par définition, chaque réduction unaire est nécessairement dans VCSPifA). 

Les réductions unaires dépendent de l’ordre total selon lequel elles sont définie. 
En effet, différents ordres impliquent différents ensembles (7, qui peuvent conduire 
à leur tour à différentes réductions unaires, comme il est illustré dans l’exemple 
suivant. 

Exemple 5.4.1. Considérons un VCSP défini sur deux variables x et y dont les 
domaines sont bi-valués. Le VCSP implique deux contraintes unaires c x et c y , une 
contrainte binaire c xy et son transposé c yx . Les divers fonctions de valuation sont : 



' 1 ' 


' 2 ’ 


’ 1 2 ’ 


’ 1 2 " 


2 

, fy ~ 

1 

7 4>xy 

2 3 

5 4*yx 

2 3 


Si nous considérons l’ordre x —>■ y, alors C = { c x , c y , c xy } et nous obtenons 



' 1 ' 


' 0 ’ 


' 2 ’ 


' 2 " 

^Px 

2 

; Py 

1 

5 4*X ® ^Px 

4 

1 0 y 4>y ® ip y 

2 


90 



















Chapitre V : Une méthode de décomposition modulaire pour les VCSPs 


Maintenant, si nous considérons l’ordre inverse, alors C = { c x , c y , c yx } et nous obte¬ 
nons 



' 1 ' 


' 0 ’ 


' 1 ' 


' 3 " 

Py 

2 

) Px 

1 

1 Ûx (fix © Px 

3 

) © Py 

3 


Ainsi, les réductions obtenues différent suivant l’ordre choisi sur les variables. 

Propriété 23. Soit P une instance de VCSP(JPt) alors la réduction unaire de V 
peut être calculée en 0(ed), où e est le nombre de contraintes et d est la taille du 
domaine des valeurs le plus large. 


Preuve. Il suffit de substituer deux contraintes unaires à chacune des contraintes 
binaires de P comme suggérée par la Propriété 16. Cette tâche peut être faite en 
0(d), et comme il y a e contraintes dans P, on a besoin de 0(ed) étapes. 


□ 

Propriété 24. Chaque instance de VCSP(Ai) est équivalente à sa réduction unaire. 


Preuve. Nous prouvons que la valuation d’une instantiation t d’un ensemble de 
variables f Ç I est la même par V et par sa réduction unaire P', C’est-à-dire 

D’après (1.1),(5.17) et la Définition 37, nous obtenons ce qui suit : 




0 </>x{tix) 

© 0 (p X y(t ix,t IV) 

xev 

c xy ec,x,yev 

0 (t>x(tlx) 

© 0 Px(t l x) ® ip y (t l y) 

xev 

c xy ec,x,yev 

0 <t> x {t\.x) 

© 0 Px{tlx)® 0 <Py(tly) 

xev 

c xy eC,x,yeV c xy eC,x,yeV 

0 <l>x(t±x) 

© 0 Px{tix)@ 0 (p x (tlx) 

xev 

c xy eC,x,yeV c yx eC,x,yeV 


<t> x (t ix)@ 0 <p x (t | x) 

c xy eC,yeV 

0 

xev 


© 

xGV 


0 (p x (t ; X) 
Cyx£-C >yeV 


□ 
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5.4.2 La décomposition des domaines 

Notre stratégie de décomposition s’appuie essentiellement sur les paires de valeurs 
modulaires (voir Définition 36). La notion de valeurs modulaires est transposée sur 
les VCSPs comme suit : 

Définition 38. Soit V = (V, D , C, S) un VCSP binaire muni d ’un ordre total sur les 
variables. Les valeurs u et v G D x sont directionnellement modulaires selon l’ordre 
considéré, (notation u ~ x v), si et seulement si pour tout c xy G C, on a u ~<f, xy v. 

Propriété 25. ~ x est une relation d’équivalence sur D x , pour tout x G X. 


Preuve. Tout d’abord, notons que la relation est l’intersection des ~ t j >sc tel que 
((x,y), (p xy ) G C. D’autre part, il est connu que l’intersection d’un ensemble de 
relations d’équivalence est une relation d’équivalence, d’où le résultat recherché. 


□ 

Propriété 26. Soit V un VCSP binaire muni d’un ordre total sur ses variables. Si 
chaque produit une seule classe d’équivalence sur D x alors V est dans VCSP(Ai). 


Preuve. Si chaque D x , x G X, est partitionné en une seule classe d’équivalence par 
la relation alors les valeurs de chaque D x sont toutes deux à deux modulaires 
pour tout x G X. Il en résulte d’après (5.11) que toutes les (p xy tel que x —> y sont 
modulaires. Et donc si (p xy G A4 alors (p yx = 4>ly £ -M- D’où toutes les fonctions de 

valuation de V sont modulaires et donc V est dans VCSP (Al). 


□ 

Toute fois, les conditions requises par la Propriété 26 sont trop fortes pour être 
satisfaites si on considère des VCSPs binaires quelconques. Pour cette raison, nous 
proposons de tirer profit de la modularité en décomposant les VCSPs à résoudre. 
Notre schéma de décomposition repose sur la propriété suivante : 

Propriété 27. Soit V un VCSP binaire, —> un ordre total sur ses variables et V 
un VCSP binaire obtenu en réduisant le domaine de chaque variable x de V à une 
classe d’équivalence de ~ x , alors V est dans VCSP(A4). 


Preuve. Ceci est obtenu par le fait que dans chaque classe d’équivalence de ~ x , x £ 
X , les valeurs sont toutes deux à deux modulaires par rapport à —K Le résultat est 
obtenu en appliquant la Propriété 26. 


□ 
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5.4.3 La solution algorithmique 

L’algorithme mmedac* décrit ci-dessous est une variante de l’agorithme medac* 
(voir Chapitre 1). Le dernier est un algorithme basé sur une recherche avec retour 
en arrière chronologique qui utilise l’algorithme de filtrage par arc-consistance di¬ 
rectionnel existentiel (edac*) [76]. mmedac* décompose progressivement le VCSP 
initial en sous-problèmes dans VCSP (.Ad). 

L’ordre des variables selon lequel les paires de valeurs modulaires sont identifiées 
est dynamique, c’est à dire, qu’il peut varier au cours de la recherche. Cet ordre est 
implicitement déterminé par la liste des variables non encore instanciées Y dont les 
éléments viennent après ceux de X \ Y dans l’ordre. 

La fonction principale dans le processus de décomposition est modequivclass 
qui est appelée à chaque nœud de l’arbre de recherche. Elle retourne la classe d’équi¬ 
valence D XjV de la valeur sélectionnée (x,v) (voir ligne 2 dans le pseudo-code). Elle 
est basée sur EQUIVCLASS (voir Section 5.3) qui est appelée pour chaque fonction de 
valuations (p xy de V telle que x est la variable courante et y est une variable déjà ins- 
tanciée, i.e., une variable dans X\Y. Le domaine de x est donc réduit à D x v (ligne 3). 

Ensuite, le sous-problème résultant est de nouveau transformé en vue d’obtenir 
progressivement un sous problème dans la classe VCSP(W). Cette étape est effectuée 
par la Procédure moddecompose (ligne 4). Cette dernière est également exécutée 
à chaque nœud de l’arbre de recherche. Elle fait appel à DECOMPOSE (voir section 
5.3) pour chaque fonction de valuations binaire é xy tel que x est la variable courante 
et y est une variable déjà instanciée. Le rôle de cette procédure est de remplacer (j> xy 
par deux fonctions de valuation unaires. L’arc-consistance directionnel existentiel 
est ensuite rétabli. Pour que l’algorithme puisse considérer les variables restantes, le 
VCSP arc-consistant directionnel existentiel ne doit pas contenir un domaine vide. 
Un appel récursif est par la suite effectué pour envisager les variables restantes (ligne 
5). Cet appel retourne le coût c* de la meilleure solution obtenue en réduisant le 
domaine de x à D XjV . Ensuite, l’algorithme annule l’effet de la réduction de D x à D x ^ v 
et élimine les éléments de D x>v de D x . 11 restaure à nouveau l’arc-consistance direc¬ 
tionnel existentiel sur le sous-problème résultant et vérifie s’il n’y a pas de domaine 
vide. Si c’est le cas, l’algorithme effectue un deuxième appel récursif qui retourne le 
coût de la meilleure solution du problème reçu en paramètre (ligne 6). 

A chaque fois que l’algorithme parvient à instancier toutes les variables, on ob¬ 
tient un VCSP dont les contraintes sont toutes un-aires, i.e., dans la classe VCSP(W). 
Le dernier appel de EDAC* le long du chemin courant fournit, par conséquent, le coût 
de la solution actuelle en temps polynomial. Le coût de cette dernière est enregistré 
dans c* au cours du dernier appel récursif (ligne 1). 
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Algorithme 22: mmedac*(7 7 ,Y,lb) : c* 


// avec V = (X,V,C,S) 
si Y = 0 alors 

1 c* -<— Z6 


sinon 

X 4- SELECT(y) 

V 4- select(U x ) 

2 v 4- MODEQUIVCLASS(’P, Y, D x ,v) 
D x 4 — v 

3 lb 4- ED AC* ("P) 
si 0 V alors 

C* MMEDAC*(P, Y \ {x}, lb) 


4 


restore(P) 

D x t D x \v 

lb 4- ED AC* (P) 
si 0 ^ V alors 

C* 4- MMEDAC* (P, F, lb) 


RESTORE(P) 


5.5 Résultats expérimentaux 

Les problèmes étudiés dans notre expérimentation sont, les problèmes Max-CSPs 
aléatoires binaires, les problèmes de coloration des graphes, les problèmes d’affec¬ 
tation de liens de radio-fréquences (RLFAP), les problèmes de planification et le 
problème de sélection génétique TagSNP. 

Nous comparons l’algorithme MMEDAC* avec medac*. L’algorithme arc-consistance 
directionnel existentiel employé pour ces deux algorithmes est edac* [74]. Les bornes 
supérieures initiales utilisées par les deux solveurs sont ceux données dans les diffé¬ 
rents fichiers tests disponibles sur [80]. 

L’heuristique d’ordre des variables est min-domain/deg [9]. Pour celle de l’ordre 
des valeurs, les deux algorithmes choisissent en priorité les valeurs qui ont le plus 
petit coût unaire projetés. Ce coût est automatiquement calculé par la procédure 
EDAC*. Les critères de valuations sont le nombre de nœuds développés, le temps 
de calcul en secondes et le coût de la meilleure solution trouvée. Nous avons égale¬ 
ment reporté, pour MMEDAC*, le nombre de nœuds-multiples développés, qui sont 
des nœuds impliquant des classes d’équivalence contenant plus d’une seule valeur. 
Tous les algorithmes ont été implémentés en C++ et exécutés sur un PC à 2 GHZ 
de fréquence et une mémoire vive de 4 Gb. 

Max-CSPs aléatoires binaires : Pour les Max-CSPs aléatoires binaires, nous 
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avons expérimenté sur les instances décrites dans [75] et générées en fonction du 
modèle des quatre paramètres [104], L’objectif est de trouver une affectation avec 
un nombre minimum de contraintes non satisfaites pour un CSP surcontraint. 

Le problème est facilement exprimé en terme de VCSP en utilisant un coût bi¬ 
naire nul pour les paires de valeurs compatibles et un coût binaire de un pour les 
paires incompatibles. 

Nous avons expérimenté sur quelques uns des instances disponibles dans [80]. 
Ces instances sont obtenues en considérant trois niveaux de densité de graphe de 
contraintes : problèmes peu dense “sparse” (S), dense (D) et complet (C) et deux 
niveaux de densité des contraintes : contraintes lâches (L) et contraintes dures (T). 
La combinaison de ces classes donne six classes de VCSPs : SL, ST, DL, DT, CL, CT. 

Nous considérons exactement deux instances de chacune des six classes décrite 
ci-dessus. La plupart des instances sont difficiles à résoudre (voir la Tableau 5.2). 
Ils ont en commun un domaine de valeur uniforme de taille dix. Le nombre de va¬ 
riables qui est indiqué immédiatement après les initiales de la classe dans le nom de 
l’instance, varie de 15 à 80. 

Le Tableau 5.2 fournit les résultats obtenus par mmedac* et medac*. Line heure 
de temps d’exécution a été consacré à chaque instance. Pour ces instances, non struc¬ 
turées, aucun des deux algorithmes ne surpasse systématiquement l’autre, mmedac* 
domine medac* sur les instances peu dense sparse et lâches. Ces problèmes sont les 
plus difficiles. Le dernier algorithme domine légèrement le premier sur les instances 
complet et lâches. Les instances restantes sont très faciles à résoudre et aucun des 
deux algorithmes ne domine l’autre. 

Coloration de graphe : Ce problème est une généralisation du problème de la 
coloration de graphe présenté dans [68]. Étant donné un graphe simple et un nombre 
de couleurs d, le but est de colorier les sommets du graphe de telle sorte que deux 
sommets adjacents Xi et x 3 n’aient pas la même couleur. 

Nous considérons précisément les graphes géométriques générés par M. Trick 
[115]. Les sommets de ces graphes sont générés dans une grille 10.000 par 10.000. 
Deux sommet sont reliés par une arête s’ils sont suffisamment proches. Les valuations 
sont inversement proportionnelles à la distance (d t . 3 ) entre les sommets. 

Nous avons précisément expérimenté sur neuf instances de ceux disponibles sur 
[80], quatre instances sont extraites du problème GEOM30a et notés GEOM30a-3 à 
6 . 

Cinq autres instances sont extraites du problème GEOM40 et sont notées par 
GEOM40-2 à 6 (voir Tableau 5.3). Les instances GEOM30a-d impliquent chacune 
30 variables et d valeurs par domaine. Les instances de GEOM40-d impliquent cha¬ 
cune 40 variables et d valeurs par domaine. Le Tableau 5.3 montre que MMEDAC* 
est plus performant que MEDAC* sur les deux instances qui sont difficiles à résoudre 
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instance 

nombre de nœuds 

temps ( 

en sec.) 


coût 

MEDAC* 

MMEDAC* 

MEDAC* 

MMEDAC* 

MEDAC 

* MMEDAC* 

SL80-10-6-60-1 

2611 

1859 

2246 

1723 

4 

4 

SL80-10-6-60-2 

2229 

1411 

2209 

1362 

2 

2 

DL30-10-25-48-1 

2.18 

9.67 

2.75 

8.37 

2 

2 

DL30-10-25-48-2 

20.51 

17.28 

10.2 

6.9 

2 

2 

CL25-10-100-18-1 

154 

262 

252 

252 

2 

2 

CL25-10-100-18-2 

103 

467 

115 

370 

2 

2 

ST25-10-21-85-1 

5.00 

6.93 

2.67 

4.38 

19 

19 

ST25-10-21-85-2 

1.60 

4.49 

2.01 

6.03 

20 

20 

DT25-10-25-87-1 

8.93 

13.11 

4.53 

5.24 

28 

28 

DT25-10-25-87-2 

16.87 

14.80 

7.12 

6.94 

30 

30 

CT15-10-100-93-1 

55.76 

65.92 

32.42 

36.37 

73 

73 

CT15-10-100-93-2 

71.32 

73.59 

43.04 

38.60 

74 

74 


instance 

nœuds-multi 

moyenne des tailles de domaine 

SL80-10-6-60-1 

569 

5.49 

SL80-10-6-60-2 

375 

5.84 

DL30-10-25-48-1 

0.48 

5.12 

DL30-10-25-48-2 

3.67 

5.11 

CL25-10-100-18-1 

49.78 

2.25 

CL25-10-100-18-2 

71.26 

2.23 

ST25-10-21-85-1 

2.18 

6.47 

ST25-10-21-85-2 

0.47 

5.34 

DT25-10-25-87-1 

2.08 

5.08 

DT25-10-25-87-2 

2.19 

2.79 

CT15-10-100-93-1 

7.80 

2.64 

CT15-10-100-93-2 

13 

3.21 


Tableau 5.2 - Résultats obtenus par MEDAC* et MMEDAC* sur les instances Max-CSPs 
aléatoires binaires. Le nombre de nœuds et nœuds multiples explorés (en milliers), le 
temps CPU en secondes et le coût de la meilleure solution trouvée sont raportés. Le temps 
d’exécution a été limité à une heure par instance. Un tiret indique que le programme n’a 
pas terminé au bout d’une heure de temps CPU. 


GEOM40-4 et GEOM40-5. Sur la dernière instance seul mmedac* a réussi à explo¬ 
rer tout l’arbre de recherche en moins d’une heure. 

Problème d’affectation de fréquence radio (RLFAP) : Ce sont des instances 
provenant d’un problème du monde réel, à savoir le problème cclar RLFAP [11] qui 
a la particularité d’être peu contraint. Le Problème d’affectation de fréquence ra¬ 
dio (RLFAP) consiste à fournir des moyens de communication entre des ressources 
spectrales (radios) tout en assurant un minimum de perturbations subies par ceux 
souhaitant communiquer dans un réseau de communication par fréquences radio. Ce 
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nombre de nœuds 

temps (en sec.) 

coût 

instance 

MEDAC* 

MMEDAC* 

MEDAC* 

MMEDAC* 

MEDAC* 

MMEDAC* 

GEOM30a-3 

6836 

80.58 

5.76 

7.79 

11 

11 

GEOM30a-4 

23.43 

21.76 

1.27 

1.16 

4 

4 

GEOM30a-5 

0.28 

0.53 

0.02 

0.04 

1 

1 

GEOM30a-6 

0.03 

0.03 

0.00 

0.00 

0 

0 

GEOM40-2 

25.35 

38.56 

1.22 

1.80 

22 

22 

GEOM40-3 

71.84 

115 

4.08 

7.15 

7 

7 

GEOM40-4 

23033 

889 

1117 

40.15 

3 

3 

GEOM40-5 

- 

10897 

- 

571 

1 

1 

GEOM40-6 

0.04 

0.04 

0.06 

0.05 

0 

0 


instance 

nœuds-multi 

moyenne des tailles de domaine 

GEOM30a-3 

12.31 

2.07 

GEOM30a-4 

1.39 

3.14 

GEOM30a-5 

0.13 

2.18 

GEOM30a-6 

0.06 

6 

GEOM40-2 

8.99 

2 

GEOM40-3 

7.99 

2.41 

GEOM40-4 

59.93 

2.00 

GEOM40-5 

125 

4.27 

GEOM40-6 

0.01 

6 


Tableau 5.3 - Résultats obtenus par MEDAC* et mmedmac* sur les instances de coloration 
de graphe. Le nombre de nœuds et nœuds multiples explorés (en milliers), le temps CPU 
en secondes et le coût de la meilleure solution trouvée sont raportés. Le temps d’exécution 
a été limité à une heure par instance. Un tiret indique que le programme n’a pas terminés 
au bout d’une heure de temps CPU. 
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problème est un problème d’optimisation combinatoire (NP-difficilc). Nous avons 
précisément expérimenté sur dix instances de ceux disponibles dans [80]. Cinq de 
ces instances sont extraites du problème celar6 et notés celar6-sub0 à 4. Les cinq 
autres instances sont extraites du problème celar7 et sont notées celar7-sub0 à 4 
(voir Tableau 5.4). Chacune de ces instances comprend 16 variables et une taille 
maximale de domaine de 44 valeurs. 

En examinant le Tableau 5.4, nous remarquons que lorsque les deux algorithmes 
terminent d’explorer tout l’arbre de recherche, mmedac* est toujours plus rapide que 
medac*. Sur l’instance celar6-sub2 seul mmedac* a réussi d’explorer toute l’arbre 
de recherche. Sur les autres instances, qui sont plus difficiles à résoudre, seulement 
sur la dernière instance, medac* a trouvé un coût plus faible. Globalement, on peut 
dire que, contrairement aux instances aléatoires, MMEDAC* surpasse MEDAC* sur 
toutes ces instances structurées sauf une : celar7-sub4. 

Problèmes de planification : Nous considérons précisément le problème de trans¬ 
port DriverLog. Ce problème implique un ensemble de conducteurs pour conduire 
des camions transportant des objets d’un emplacement de départ vers une destina¬ 
tion finale. Plus de détails sur ce problème se trouve dans [51]. 

Les dix instances sur lesquels nous avons expérimenté sont modélisés comme des 
WCSPs de densité moyenne puisque environ la moitié des contraintes possibles sont 
présentes. Chacune des instances implique des centaines de variables et une taille de 
domaine d’au plus onze valeurs. 

Le Tableau 5.5 montre que MMEDAC* est globalement meilleure que MEDAC*. En 
effet, sur les instances pour lesquelles MEDAC* domine MMEDAC*, le temps d’exécu¬ 
tion des deux algorithmes est très proche. 

Le problème de sélection TagSNP se pose dans le domaine de la génétique. 
L’objectif est de sélectionner un sous-ensemble restreint de polymorphisme nucléo¬ 
tide unique (SNP) qui capte l’information génétique critique des individus d’une 
même espèce. Le problème est étroitement liée à un problème NP-difficile, celui de 
couverture par des sommets de graphes pondérés, et peut être exprimé comme un 
VCSP binaire tel qu’il est présenté dans [1]. Le problème résultant est plutôt struc¬ 
turé comme le graphe sous-jasent à chacune de ses instances est construit sur la base 
d’une distance entre chaque paire de SNP. Nous expérimentons précisément sur plu¬ 
sieurs instances de tagsnp_r08 qui sont valables dans [2], Le nombre de variables de 
ces instances varie de 33 à 93 les tailles de leurs domaines varient de 15 to 74. 

La Figure 5.2 compare les résultats obtenus par les algorithmes MMEDAC* et 
MEDAC*. Pour ces instances structurées, MMEDAC* domine MEDAC* sur toutes les 
instances sauf une. L’instance pour laquelle MEDAC* domine MMEDAC* est (25023), 
elle est très facile à résoudre, comme elle est résolue en moins d’une seconde par les 
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instance 

nombre de nœuds 

temps ( 

en sec.) 


coût 

MEDAC* 

MMEDAC* 

MEDAC* 

MMEDAC* 

MEDAC* 

MMEDAC* 

celar6-sub0 

0.11 

0.13 

0.38 

0.24 

159 

159 

celar6-subl 

13235 

2253 

2353 

1582 

2669 

2669 

celar6-sub2 

- 

448 

- 

1949 

2746 

2746 

celar6-sub3 

- 

- 

- 

- 

3398 

3177 

celar6-sub4 

- 

- 

- 

- 

6307 

5643 

celar7-sub0 

0.48 

0.20 

0.40 

0.18 

10310 

10310 

celar7-subl 

66.42 

45.46 

87.65 

45.68 

142640 

142640 

celar7-sub2 

753 

355 

1422 

541 

173252 

173252 

celar7-sub3 

- 

- 

- 

- 

242842 

223460 

celar7-sub4 

- 

- 

- 

- 

292639 

372448 


instance 

nœuds-multi 

moyenne des tailles de domaine 

celar6-sub0 

0.04 

8.33 

celar6-subl 

326 

5.64 

celar6-sub2 

35.18 

3.09 

celar6-sub3 

- 

- 

celar6-sub4 

- 

- 

celar7-sub0 

0.08 

7.07 

celar7-subl 

9.86 

3.99 

celar7-sub2 

50.27 

3.66 

celar7-sub3 

- 

- 

celar7-sub4 

- 

- 


Tableau 5.4 - Résultats obtenus par medac* et mmedac* sur les instances extraites des 
problèmes réels : celarb et celar7. Le nombre de nœuds et nœuds multiples explorés (en 
milliers), le temps CPU en secondes et le coût de la meilleure solution trouvée sont rapor- 
tés. Le temps d’exécution a été limité à une heure par instance. Un tiret indique que le 
programme n’a pas terminé au bout d’une heure de temps d’exécution. 
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instance 

nombre de nœuds 

temps (en sec.) 

coût 

MEDAC* 

MMEDAC* 

MEDAC* 

MMEDAC* 

MEDAC* 

MMEDAC* 

driverlog02ac 

14.96 

14.05 

20.29 

19.86 

2010 

2010 

driverlog02bc 

11.28 

10.92 

16.35 

15.83 

2085 

2085 

driverlog02cc 

5.63 

6.93 

9.20 

12.72 

2428 

2428 

driverlog03 

2.14 

4.59 

3.69 

6.51 

1225 

1225 

driverlog04ac 

3.78 

8.71 

5.31 

12.37 

1790 

1790 

driverlog04bc 

3.83 

7.49 

5.41 

11.46 

1921 

1921 

driverlog04cc 

3.45 

2.94 

5.03 

4.71 

2932 

2932 

driverlog05ac 

3.33 

3.14 

14.17 

10.50 

1350 

1350 

driverlog05bc 

3.34 

3.13 

14.30 

11.07 

1344 

1344 

driverlog05cc 

3.43 

2.72 

14.23 

9.18 

2659 

2659 

driverlogsOG 

294 

22.12 

271 

20.63 

1055 

1055 

driverlog08ac 

37.51 

7.48 

318 

50.89 

2220 

2220 

driverlog08bc 

37.92 

7.23 

273 

50.27 

2282 

2282 

driverlog08cc 

36.02 

7.10 

281 

60.06 

4052 

4052 


instance 

nœuds-multi 

moyenne des tailles de domaine 

driverlog02ac 

5.70 

2.09 

driverlog02bc 

4.26 

2.12 

driverlog02cc 

2.60 

2.12 

driverlogs03 

1.46 

2.06 

driverlog04ac 

3.89 

2.31 

driverlog04bc 

3.12 

2.28 

driverlog04cc 

1.24 

2.17 

driverlog05ac 

1.55 

2.12 

driverlog05bc 

1.52 

2.12 

driverlog05cc 

1.32 

2.29 

driverlogsOG 

14.79 

2.09 

driverlog08ac 

2.14 

2.01 

driverlog08bc 

2.08 

2.01 

driverlog08cc 

2.09 

2.00 


Tableau 5.5 - Résultats obtenus par MED AC* et mmedac* sur les instances extraites du 
problème de planification driverlog. Le nombre de nœuds et nœuds multiples explorés 
(en milliers), le temps CPU en secondes et le coût de la meilleure solution trouvée sont 
raportés. Le temps d’exécution a été limité à une heure par instance. Un tiret indique que 
le programme n’a pas terminé au bout d’une heure de temps CPU. 
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deux algorithmes. En plus, sur les deux premières instances, mmedac* est deux fois 
plus rapide que MEDAC*. Finalement, MMEDAC* a résolu plus d’instances en quatre 
heures de temps CPU. 

5.6 Conclusion 

Dans ce chapitre, nous avons présenté un algorithme de décomposition de VCSPs 
binaires qui vise à accélérer la résolution de ces derniers. Le processus de décom¬ 
position proposé instancie les variables en réduisant leurs domaines respectifs à des 
sous-ensembles de valeurs modulaires au lieu de les réduire à des simple valeurs 
comme c’est le cas avec les algorithme de résolution classiques. 

L’algorithme proposé produit, à chaque feuille de l’arbre de recherche, un sous- 
problème dont les fonctions de valuations sont toutes modulaires. Ces sous-problèmes 
sont résolus en temps polynomial. 

Une extension naturelle de ce travail est de dériver un schéma de décomposition 
similaire pour les VCSPs n-aires puisque la propriété 16 serait facilement généra¬ 
lisable au cas n-aire. Un schéma de décomposition basé sur la classe des fonctions 
sous-modulaires sera considéré dans le chapitre suivant. 
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instance 

nombre de nœuds 

temps (en sec.) 

coût 

MEDAC* 

MMEDAC* 

MEDAC* 

MMEDAC* 

MEDAC* 

MMEDAC* 

10734 

229 

8.03 

157 

6.01 

150538 

150538 

11053 

40.53 

0.41 

180 

1.23 

140279 

140279 

12929 

- 

4.46 

- 

121 

1165000 

666797 

13106 

1.96 

0.11 

9.68 

1.06 

278934 

278934 

14481 

- 

42.96 

- 

295 

827112 

518137 

16864 

1.47 

0.15 

3.37 

0.29 

120356 

120356 

17289 

10.16 

0.18 

2.50 

0.04 

72638 

72638 

18827 

- 

0.38 

- 

23.70 

1611272 

920935 

19684 

0.06 

0.04 

0.09 

0.06 

96463 

96463 

20927 

8.50 

0.32 

3.76 

0.12 

82339 

82339 

21758 

0.27 

0.07 

5.77 

1.42 

243734 

243734 

22707 

- 

8.5 

- 

46.20 

287211 

287211 

25023 

0.05 

0.22 

0.08 

0.21 

132066 

132066 

26578 

- 

488 

- 

3459 

852467 

662633 

27116 

39.40 

2.66 

101 

4.41 

154277 

154277 

28188 

- 

1.10 

- 

4.39 

258747 

258747 


instance 

nœuds-multi 

moyenne des tailles de domaine 

10734 

7.33 

6.68 

11053 

0.37 

11.73 

12929 

4.17 

20.38 

13106 

0.10 

13.19 

14481 

40.26 

10.30 

16864 

0.40 

11.67 

17289 

0.12 

4.88 

18827 

0.34 

23.70 

19684 

0.03 

8.15 

20927 

0.30 

5.04 

21758 

0.06 

18.19 

22707 

8.5 

14.45 

25023 

0.2 

9.62 

26578 

453 

9.68 

27116 

2.17 

8.32 

28188 

1.06 

10.21 


Figure 5.2 - Résultats obtenus par MEDAC* et mmedac* sur les instances extraites du 
problème réel TagSnP. Le nombre de nœuds et nœuds multiples explorés (en milliers), le 
temps CPU en secondes et le coût de la meilleure solution trouvée sont raportés. Le temps 
d’exécution a été limité à une heure par instance. Un tiret indique que le programme n’a 
pas terminé au bout d’une heure de temps CPU. 
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Résumé 

Plusieurs problèmes combinatoires peuvent être formulés comme des problèmes de Satisfaction 
de Contraintes Valuées (VCSPs), où les contraintes sont définies par des fonctions de valuations 
pour tenir compte de divers degrés de cohérence. Ce problème est, en général, NP-difficile. Toute¬ 
fois, des versions polynomiales peuvent être obtenues en se limitant à des fonctions de valuations 
avec des propriétés spécifiques. C'est le cas des instances impliquant des fonctions binaires et 
sous-modulaires. 

Dans ce chapitre, nous identifions une nouvelle classe sous-modulaire de VCSPs binaires qui 
étend la classe polynomiale de VCSP décrite dans [14]. Cette nouvelle classe, que nous appelons 
VCSPs binaires sous-modulaires par permutations, est obtenue en appliquant indépendamment 
des permutations sur chaque domaine de valeur. Nous montrons que ces sous-problèmes peuvent 
être identifiés en 0(n 2 d 4 ) étapes et résolus en 0(n 3 d 3 ) étapes. 

Nous proposons également un algorithme de décomposition d’un VCSP n’appartenant pas à 
cette super-classe en des sous-problèmes de VCSPs de cette classe. 


6.1 Introduction 

Ce chapitre a fait l’objet des publications [56, 57, 59]. 

Nous commençons par rappeler qu’en exprimant les VCSPs sous-modulaires en 
termes de recherche d’une coupe minimale d’un graphe orienté et pondéré il est pos¬ 
sible de les résoudre en 0{n 3 d 3 ) étapes où n est le nombre de variables et d est la 
taille du domaine de valeur le plus large [14]. Plus de détails sur ce point ont été 
fournis dans le chapitre 3. 

Néanmoins, les VCSPs résultant de cas réels ne se limitent pas à des fonctions 
sous-modulaires. Dans de tels cas, pouvons nous procéder de manière plus efficace 
qu’une recherche exhaustive tout en garantissant l’optimalité ? Y a-t-il un moyen 
d’exploiter plus efficacement la sous-modularité ? 

Dans [108], Schlesinger a identifié une super-classe de VCSPs binaires sous- 
modulaires qui peut être identifiée en temps polynomial. Les instances, dans cette 
super-classe qui utilise l’addition comme opérateur d’agrégation, sont obtenues des 
instances sous-modulaires en permutant séparément les valeurs. L’auteur a, éga¬ 
lement ,indiqué que l’algorithme fonctionne seulement lorsque les valuations sont 
finies. Cependant, l’utilisation des valuations infinies est essentielle parce qu’elles 
fournissent le moyen le plus naturel pour exprimer l’insatisfiabilité. Cette dernière 
est essentielle dans la modélisation des problèmes découlant du monde réel. 

La classe des VCSPs binaires sous-modulaires par permutations est, évidement, 
plus large que la classe obtenue “sans-permutation”. Cependant, elle ne peut pas 
aussi encoder toutes les situations découlant du monde réel. Par exemple, plusieurs 
problèmes de vision par ordinateur impliquent des fonctions d’énergie non sous- 
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modulaires qui ne peuvent pas être transformées en des fonctions sous-modulaires 
par le biais de permutations des valeurs des domaines [3]. 

En outre, la décomposition du problème est une stratégie largement utilisée en 
intelligence artificielle, qui vise à améliorer la résolution des problèmes en les dé¬ 
composant en un ensemble de sous-problèmes plus faciles à résoudre. Motivé par 
les résultats obtenus sur les CSPs classiques et les VCSPs [91, 93], nous présentons, 
dans ce chapitre, un schéma de décomposition qui tire profit des fonctions sous- 
modulaires, même si le problème à résoudre n’est pas définie uniquement à l’aide 
de ces fonctions. Le schéma de décomposition est implémenté par une recherche 
arborescente qui permet de résoudre, à chaque feuille de l’arbre de recherche, un 
sous-problème dont les fonctions de valuations sont toutes sous-modulaires. De tels 
sous-problèmes peuvent être résolus en temps polynomial [14]. 

Le chapitre est organisé comme suit : la Section suivante présente les définitions 
et les notations nécessaires à la présentation de notre approche. La Section 6.3, 
présente le schéma de décomposition proposé pour la résolution de notre problème 
et l’algorithme correspondant. Divers résultats expérimentaux sont présentés dans 
la Section 6.5 et quelques travaux connexes sont décrits dans la Section 6 . 6 . Nous 
concluons dans la Section 6.7. 

6.2 Définitions et notations 

Rappelons, tout d’abord, qu’un opérateur d’agrégation © est dit strictement mo¬ 
notone si pour tout a, /3 ,7 dans E tel que a -< (3 et 7 7 ^ T, nous avons «©7 -< /3© 7 . 

Considérons la structure de valuation S = (E,©, t<). Si a, (3 G E sont tel que 
a 7 ! f3 et il existe une valuation 7 E E qui vérifie a © 7 = /3, alors 7 est connue 
comme la différence entre (3 et a. Rappelons que, la structure de valuation S est 
propre si pour chaque paire de valuations a, (3 G E , tel que a 7! /3, il existe une 
différence maximale entre /3 et a. Cette différence maximale unique de /3 et a est 
dénotée par (3 © a. 

Dans ce chapitre, nous traitons les VCSPs binaires, qui sont, des VCSPs avec 
des contraintes unaires et binaires seulement. 

Nous supposons également que toutes les fonctions de valuation <p : D x D' —»• E 
utilisées lors de l’exposition de notre approche sont supposées être arc-consistantes, 
ce qui signifie que 


Vu G D , 3 u' G D' , tel que <p(u,u') -< T (6.1) 

Ceci ne présente, en aucun cas, une limitation de notre approche, puisque toute 
fonction de valuation non arc-consistante peut se transformer en une autre qui est 
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arc-consistante simplement en supprimant les valeurs qui ne remplissent pas la condi¬ 
tion (6.1) ci-dessus. Dans le cadre VCSP, l’application d’une telle transformation sur 
chaque fonction de valuation se fait en temps polynomial. Cette tâche se traduit par 
la réduction des domaines de valeurs sans changer l’ensemble de toutes les solutions 
possibles au problème. En effet, les valeurs rejetées ne peuvent participer à aucune 
solution optimale. Le problème qui en résulte est appelé arc-consistant. 

Dans le cas binaire, les fonctions sous-modulaires peuvent être définies comme 
suit : 

Définition 39. Une fonction binaire (f> : D x D' —> E est sous-modulaire si et 
seulement si pour toutes u,v G D et pour toutes u', v' G D' tel que u < v et v! < v', 
nous avons 

<t>(u,v!) © 4>(v,v r ) S 4>(u,v') © 4>{v,u') (6.2) 

Exemple 6.2.1. Considérons la structure de valuation S = (N,+,<). Pour tout 
entier À > — 2 , la famille des fonctions binaires <j)\ est définie sur n’importe quel 
sous-ensemble fini et non vide D de N par v) = u 2 + v 2 + A uv. Nous montrons 
que 4 >a est sous-modulaire si et seulement si A < 0. 

Soit u, v, u', v' E D tel que u < v, u' < v'. En observant que Xuu'+Xvv' < Xuv'+Xvu' 
pour toute A < 0, nous obtenons 

4>\(u, u') + 4>\(v, v') = u 2 + u' 2 + Xuv! + v 2 + v' 2 + Xvv' 

< u 2 + v' 2 + A uv' + v 2 + u' 2 + A vu' 

< (j)\(u,v') + <j>\{v,v!) 

D’autre part, si A > 0 alors nous avons </>a( 0, 0) + 1) = 2 + A qui est supérieure 

à 0 a(O, 1) + 0a( 1, 0) = 2. Par conséquent, si A > 0 alors <f>\ n’est pas sous-modulaire. 

Les fonctions sous-modulaires sont largement étudiées comme elles sont très uti¬ 
lisées dans de nombreux problèmes d’optimisation discrète [94, 96, 109, 114, 121]. En 
particulier, les VCSPs qui n’utilisent que des fonctions de valuations sous-modulaires 
et impliquant l’opérateur d’agrégation + peuvent être résolus en temps polynomial. 
Dans le cas binaire, par exemple, de tel problèmes peuvent être résolus en les ex¬ 
primant comme un problème de recherche de la coupe de poids minimal dans un 
graphe orienté et pondéré [14]. Dans ce qui suit, un VCSP qui ne fait intervenir que 
des fonctions de valuations sous-modulaires sera appelé VCSP sous-modulaire. 

6.3 Décomposition de VCSPs 

L’approche décrite dans [108] peut échouer à trouver un ensemble de permutation 
de domaine qui transforme un VCSP binaire donné en un VCSP dont toutes les fonc¬ 
tions de valuations sont sous-modulaires. Cela se produit tout simplement lorsque le 
problème n’est pas dans la classe de VCSPs binaires sous-modulaires par permuta¬ 
tions. Cependant, l’idée de permutations de domaines peut être utilisée pour dériver 
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un schéma de décomposition de problème dans le but est d’accélérer la résolution 
des VCSPs binaires. 

Notre approche consiste à décomposer le problème à résoudre en un ensemble de 
sous-problèmes qui ne font intervenir que des fonctions de valuations sous-modulaires. 
Ces problèmes sont polynomiaux lorsque l’opérateur d’agrégation est strictement 
monotone [14]. Nous avons déjà proposé un schéma de décomposition similaire dans 
le précédant chapitre, sauf que la classe des fonctions de valuations cible était celle 
des fonctions de valuations modulaires. Intuitivement, un schéma basé sur la décom¬ 
position utilisant la super-classe des fonctions sous-modulaires est plus générale. 

Le schéma proposé s’intégre au sein d’un algorithme de recherche en profondeur 
d’abord qui explore un arbre binaire de recherche. Il consiste à appliquer l’une des 
deux actions suivantes à chaque nœud de l’arbre de recherche : 

1 . A chaque nœud interne, le domaine de la variable courante est partitionné en 
deux : une chaîne 1 et un résidu qui est égale au domaine privé de la chaîne. 
La chaîne peut être égale à tout le domaine. Dans ce cas, le résidu est le 
vide. L’algorithme réduit le domaine de la variable courante à la chaîne et 
explore le sous-arbre résultant de cette action, puis réduit le domaine de la 
variable courante au résidu et explore le sous-arbre de recherche résultant de 
ce deuxième choix. 

Les valeurs dans la chaîne sont éventuellement permutées. L’objectif de cette 
permutation est de transformer chaque fonction de valuation, liant la chaîne 
de la variable courante à une chaîne d’une variable déjà instanciée, en une 
fonction sous-modulaire (voir la Définition 39). C’est la clé de notre schéma de 
décomposition. 

Les permutations des valeurs de la chaîne courante sont (en partie) déduites de 
celles des chaînes affectées aux variables passées. Idéalement, une permutation 
unique du domaine courant est déduite. Dans ce cas, la chaîne correspondra 
au domaine courant et le résidu sera vide. Il peut aussi arriver que plusieurs 
permutations sont possibles. Dans une telle situation, une permutation unique 
devrait être (arbitrairement) choisie. Inversement, il peut arriver qu’aucune 
permutation n’est compatible avec les chaînes affectées aux variables passées. 
Dans ce cas, la chaîne ne pourra pas inclure toutes les valeurs du domaine, et 
ce dernier doit être partitionné. Dans le cas extrême, la chaîne sera réduite à 
une seule valeur. 

2. A chaque nœud feuille, un algorithme de calcul de la coupe minimale d’un 
graphe orienté et pondéré est exécuté afin de résoudre le problème sous-modulaire 
disponible au niveau de ce nœud. 

L’exécution de l’action 1 à chaque nœud interne de l’arbre de recherche produit, 
à chaque feuille de l’arbre de recherche, un sous-problème sous-modulaire. De tels 

1. L’utilisation du terme chaîne sera justifiée plus tard 


107 




Chapitre VI : Permutation de domaine pour les VCSPs binaires 


sous-problèmes peuvent être résolus par l’algorithme proposé dans [14] en temps 
polynomial. 

La description informelle du schéma de décomposition donnée ci-dessus montre 
que le processus qui sera effectué au niveau des nœuds internes de l’arbre de recherche 
est composé de deux phases : une phase de partitionnement de domaine et une phase 
de permutation de valeurs. Ces deux phases seront détaillées ci-dessous. 

6.3.1 Partitionnement de domaines 

La relation binaire suivante permettra de tirer profit des fonctions sous-modulaires 
même si les fonctions de valuations intervenant dans le problème initial ne sont pas 
sous-modulaires par permutation. 

Définition 40. Soit V un VCSP binaire muni d'un ordre total sur ses variables et 
soit x une variable de V. La valeur u G D x est dite directionnelle sous-modulaires 
à v G D x (notation u -< x v) si pour toute contrainte ((x,y) ,(p) G C et pour tout 
u’,v' G D y , u' < v', on a 

4>(u,u') © (f)(v,v') A 4>(u,v') © 4>(v,u') (6.3) 

Dans ce qui suit, nous montrons que -< x est un préordre sur D x . A cette fin, nous 
utilisons la relation d’équivalence (voir chapitre 5) définie sur D x par : u ~ x v si 
et seulement si pour toute contrainte (( x,y ) , (p) G C, on a 

</>(«,«') © 4>(v,v') = 4>(u,v') © 4>(v,u') 

Notons que u ~ x v si et seulement si u -< x v et v -< x u. 

Nous avons déjà montré, dans le chapitre précédant, (Propriété 19) que ~ x est 
une relation d’équivalence sur D x . 

Propriété 28. ~< x est un préordre sur D x . 


Preuve. < x est trivialement réflexive. Montrons que cette relation est aussi transi¬ 
tive. À cette fin, supposons l’inverse et procédons pour obtenir une contradiction. 

Soit u,v,w G D x tel que u -< x v et v P x w mais u -Ji x w. Cela implique que, pour 
tout u’,v’ G D y tel que u 1 < v', nous avons 

cp(u, u’) © cj)(v, v') A 4>(u,v r ) © cp(v,u') (6.4) 

4>(v, u') © (p(w, v') P 4>{v,v') © 4>(w,u') (6.5) 

et, comme © est un ordre total, il doit exister u',v' G D y , u' < v' tel que 

4>(u,u') © 4>(w,v') >- 4>(u,v') © 4>(w,u') (6.6) 
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Nous montrons que (6.4) et (6.5) ne peuvent tenir que si on a une égalité. En effet, si 
le membre à gauche de (6.4) est plus petit que le membre à droite alors l’agrégation 
du membre à droite de (6.6) avec le membre à gauche de (6.4) et vice, versa nous 
donne 


0 (u,w') © (f>(w,v') >- (j){v,v') © (j>{w,v!) 

qui est en contradiction avec (6.5). Une contradiction similaire peut être dérivée si 
nous supposons que les deux membres de (6.5) ne sont pas égaux. Ainsi, pour tous 
les u',v' G D y tel que u' < v', (6.4) et (6.5) peuvent tenir si on a une égalité. Par 
symétrie, nous pouvons déduire que (6.4) et (6.5) tiennent aussi si on a une égalité 
pour tous u',v' G D y tel que u' > v'. Ceci signifie que u ~ x v et v ~ x w. Par la 
Propriété 19, nous pouvons déduire que u ~ x w, puis pour tous u', v ' G D y , nous 
devons avoir </>(«, u') © 4>(w, v') = </>(«, v') © u') qui est en contradiction avec 

(6.6). Ainsi, A x est transitive. 


□ 


6.3.2 La plus longue partition en chaîne 

L’idée est d’essayer simultanément le maximum de valeurs possibles le long de 
chaque branche de l’arbre de recherche. En ce qui suit, nous montrons que cela est 
possible pour n’importe quel sous-ensemble de paires de valeurs sous-modulaires. 

Le préordre A x peut être utilisé pour partitionner D x en deux sous ensembles : 
une chaîne de valeurs qui sont totalement ordonnées par < x et le reste de D x . Plus 
précisément, -< x induit un ordre partiel < x sur l’ensemble quotient de D x par ~ x , 
qui est D x j Cet ordre partiel est défini par 

[u] < x [u] si et seulement si u -< x v 

ofi [u] désigne la classe d’équivalence de u. Pour identifier le plus large sous ensemble 
de paire de valeurs sous-modulaires, nous avons à déterminer la plus longue chaîne 
de l’ordre partiel < x . Trouver une telle chaîne peut être exprimée comme le problème 
de trouver le plus long chemin dans un graphe acyclique ordonné. Ce problème peut 
être résolu en temps polynomial en utilisant l’algorithme de plus court chemin avec 
des poids négatifs appropriés. D’où, à chaque nœud interne de l’arbre de recherche, 
un sous-programme est intégré à l’algorithme afin de calculer la plus longue chaîne de 
valeurs qui sont totalement ordonnées par < x . La première étape dans le calcule de 
la plus longue chaîne consiste à calculer la relation et l’ensemble quotient de D x 
par ~ x , D x j X • Ceci peut être accompli en 0[nd 3 ) étapes en utilisant l’algorithme 
proposé dans [93] (voir la Propriété 5). Le facteur n dans 0{nd 3 ) est du au fait que 
chaque variable peut être impliquée dans, au maximum, n — 1 contraintes. L’étape 
la plus longue dans l’algorithme de partitionnement est le calcul de l’ordre partiel 
< x . Selon la Définition 40, ceci nécessite, dans le plus pire cas, à effectuer 0{nd A ) 
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testes de sous-modularité, par exemple les testes donnés par (6.3). Supposons le 
dernier teste est achevé en temps constant, nous obtenons un temps de complexité 
de 0{nd A ) pour cette étape. 

Ensuite, l’algorithme construit un graphe acyclique orienté G = (V,E), où l’en¬ 
semble de nœuds V contient les éléments de D x / et l’ensemble des arcs E contient 
un arc ([«], [u]) si et seulement si [u] ^ [u] et [w] < x [u]. Les arcs sortants du nœud 
correspondant à [w] ont le poids de —card[u], où card[u ] est le nombre de valeurs 
dans [«]. Ce poids favorise la plus large classe d’équivalence. Nous devons aussi ajou¬ 
ter un nœud de départ s et un nœud terminal t. Le nœud de départ s est connecté 
à tous les nœuds dont les degrés entrant valent zéro, et les arcs assurant ces connec¬ 
tions sont pondérés par 0. Le nœud terminal est connecté à chaque nœud [u] dont 
les degrés sortant valent zéro, et les arcs assurant ces connections sont pondérés par 
—card[u]. La construction de G peut être réalisée en 0(d 2 ) étapes. L’algorithme de 
Bellman-Ford [5] est, donc, appliqué à G. Cet algorithme nécessite 0(d 3 ) étapes, où 
d est la borne supérieure pour le nombre de nœuds dans G. La plus longue chaîne 
est extraite de G en incluant les éléments d’une classe [u] dans la plus longue chaîne 
lorsque [u] est dans le plus long chemin obtenu par l’algorithme de Bellman-Ford. Le 
coût de la dernière étape est 0(d 2 ). Le temps de complexité global de l’algorithme 
calculant la plus longue partition en chaîne est, donc, 0(nd A ). 


6.3.3 La partition en chaîne de Dilworth 

Le préordre -< x peut être utilisé pour partitionner D x en chaines de valeur qui 
sont totalement ordonnées par -< x . Plus précisément, -< x induit un ordre partiel < x 
sur l’ensemble de quotient D x par ~ x , qui est D x j ~ x . Cet ordre partiel est défini 
par 


M <* b] ssi u< x v 

ofi [u] désigne 1a, classe d’équivalence de u. Pour un seul ensemble partiellement or¬ 
donné, on peut avoir plusieurs partitions en chaines. En théorie d’ordre, la partition 
en chaine optimale d’un ensemble partiellement ordonné est connue par la partition 
en chaine de Dilworth (DCP) [30]. C’est une partition en chaîne qui comporte le 
nombre minimum de chaînes sur l’ensemble des partitions en chaînes possibles. La 
taille d’une telle partition définit 1a, largeur de l’ordre partiel. Le problème de trouver 
la DCP d’un ensemble partiellement ordonné peut être résolu en temps polynomial 
en l’expriment comme un problème du couplage maximum dans un graphe biparti 

[40]. 

Motivé par le fait que l’algorithme doit avoir, à chaque nœud non-feuille, autant 
de choix qu’il y a des chaînes dans la DCP du domaine actuel, nous proposons de 
calculer un DCP à chaque nœud non-feuille de l’arbre de recherche. Le but est de 
minimiser le nombre de nœuds qui seront développées lors de la recherche. 
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La première étape dans le calcule de la DCP consiste à calculer la relation et 
l’ensemble de quotient de D x obtenu par D X J Ceci peut être accompli en 
0(nd 3 ) étapes, où n est le nombre de variables et d est la taille du domaine le plus 
large, en utilisant l’algorithme que nous avons proposé dans [93] (voir la Propriété 
5 de [93]). Le facteur n dans 0{nd 3 ) est du au fait que chaque variable peut être 
impliquée dans, au plus, n — 1 contraintes. L’étape la plus longue dans l’algorithme 
de partitionnement est le calcule de l’ordre partiel < x . En s’accordant à la Définition 
39, à calculer dans le pire cas O (nd 4 ) tests de sous-modularité, ces tests sont données 
par (6.3). En supposant que le dernier test réalisable en un temps constant, nous 
obtenons un temps de complexité de 0(nd 4 ) pour cette étape. Ensuite, l’algorithme 
construit un graphe biparti G = (V, U, E) où V — U — D x / et l’ensemble des 
arêtes E contient une arrête [u],[u] si et seulement si [n] ^ [n] et [n] < x [n]. La 
construction de G peut se faire en 0(d 2 ) étapes. Un algorithme de couplage maxi¬ 
mum est donc appliquée à G. Ceci peut être accompli par l’algorithme décrit dans 
[61] qui consomme 0(d 2 ' 5 ) étapes. Les chaînes de la partition sont extraites du cou¬ 
plage maximum en incluant les éléments d’une classe [ u] et ceux de la classe [n] dans 
une même chaîne chaque fois qu’il y a une arête [u], [n] dans le couplage maximum. 
Les valeurs représentées par des sommets qui ne sont pas membres du couplage 
maximum donnent lieu à des chaînes à valeur unique. Le coût de la dernière étape 
est 0(d 2 ). La complexité en temps de l’ensemble de la phase de partitionnement en 
chaîne est, par conséquent, 0{nd A ). 

Notons par ô x la plus longue chaîne ou une chaîne de la partition en chaîne de 
Dilworth de D x obtenue selon le préordre -< x . Alors, nous avons ce qui suit 

Propriété 29. ~< x est un préordre total sur S x . 


Preuve. -< x est réflexive et transitive sur D x alors il est aussi réflexive et transitive 
sur chaque sous ensemble de D X: et en particulier sur 5 X . -< x est un ordre total sur 
5 X du fait que ô x est une chaîne, et donc pour toutes u,v G S x , nous avons u -< x v 
ou v -< x u. 


□ 


6.3.4 Permutation de valeurs 

La phase de permutation de valeurs prend comme entrée une chaîne S x du do¬ 
maine en question, (D x ). Son rôle est de permuter les valeurs de S x , afin que, pour 
chaque contrainte ((x,y) ,(p) tel que x -> y, la restriction de (j) à la chaîne permutée 
produit une fonction sous-modulaire. Plus précisément, la tâche est de trouver une 
permutation n des éléments de ô x telle que la fonction <p, obtenue à partir de (j) en 
restreignant D x à ô x et en appliquant tt à ô x , est sous-modulaire. La fonction ip est, 
donc, définie pour toutes ( u , v) G 8 X x D y par 

<p(n(u),v) = <j)(u,v) 
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ce qui implique que 


cp(u,v) = (j>( 7r 1 (m),u) (6.7) 

Nous montrons que si 7r est choisie de tel manière qu’elle vérifies l’implication sui¬ 
vante alors tp est sous-modulaire. 

Si u -< x v => n(u) < 7 r(v) (6.8) 

Soit 7r est une permutation de ô x qui vérifies (6.8). Supposons que n’est pas 
sous-modulaire. Alors il doit exister u, v G S x et u', v' G D y tel que u < v, u' < v' et 

ip(u,u') © ip(v,v') ^ tp(u,v') © tp(v,v!) 

Par (6.7), nous obtenons 

0(7r _1 (w),w') © 0(7r _1 (u), v') ^ 0(7T _1 (m), v') © 0(7T _1 (u), u') 

Ceci implique que 7r _1 (w) 7r^ 1 (u) et depuis -< x est un préordre total sur S x nous 

déduisons que 7 r _1 (u) -< x tx _1 (m). Par (6.8) ceci implique que v < u, qui est en 
contradiction avec la condition sur u et v. D’ou, est sous-modulaire. 

Etant donné ô x , il est toujours possible de trouver une permutation 7r qui vérifies 
(6.8) pour toutes u, u G S x . U 11 simple algorithme qui calcule une telle permutation 
peut être décrit comme suit : 

1. Trier les éléments de la chaîne S x dans un ordre croissant selon pour obtenir 
une chaîne ordonnée ô x . 

2. Trier les éléments de la chaîne 5 X dans l’ordre croissant naturel pour obtenir 
une chaîne ordonnée 5 X . 

3. Puis choisir 7r pour quelle soit une permutation qui associe, à chaque élément 
de ô x , l’élément qui occupe la même position dans S x . 

L’application d’une telle permutation sur chaque chaîne extraite au niveau des 
nœuds internes produit un sous-problèmes sous-modulaire aux niveau des feuilles 
de l’arbre de recherche. La dernière étape dans le calcul d’une solution optimale 
(temporaire) consiste à appliquer, au niveau des feuilles, les permutations inverses 
de celles appliquées au niveau des nœuds internes afin de récupérer les solutions du 
problème initiale (sans permutation). 

6.4 L’algorithme de résolution 

Nous proposons deux algorithmes de résolution, le premier se base sur la plus 
longue partition en chaîne et le second se base sur la partition en chaîne de Dilworth. 
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6.4.1 Algorithme utilisant la plus longue partition en chaîne 

Notre algorithme de décomposition lcpmac* (pour Longest Chain MAC*), dont 
le pseudo-code est donné par l’Algorithme 23, est obtenu en intégrant la stratégie 
de décomposition décrite ci-dessus dans l’algorithme standard, nommée MAC* [75]. 
lcpmac* traite, à chaque nœud de l’arbre de recherche, un sous-problème dénoté 
par V. Si le nœud traité est un nœud feuille, le sous-problème courant est alors 
sous-modulaire. L’algorithme appelle, par la suite, la fonction MINCUT, (voir ligne 
1 de lcpmac*). Rappelons que, selon [14], le problème est exprimé en problème de 
recherche de coupe minimale dans un graphe pondéré. Il est résolu et une solution 
optimale s est renvoyée. Le coût de la dernière solution déjà trouvée est comparée 
à la solution courante. S’il y a une amélioration, alors les permutations inverses 
utilisées au niveau des nœuds internes, sont appliquées à s et la solution obtenue est 
stockée dans s*. 

Si, toutefois, le nœud traité est un nœud interne alors l’algorithme sélectionne 
une variable non encore instanciée, x, et calcule le préordre ^ x comme nous l’avons 
déjà expliqué (voir ligne 2). Ensuite l’algorithme appelle la fonction LONGCHAIN. Le 
rôle de cette dernière consiste à calculer, pour le domaine courant, la plus longue 
chaîne formée de paires de valeurs sous-modulaires. Ensuite, l’algorithme calcule 
une permutation de cette chaîne. Au lieu de réduire le domaine courant à une va¬ 
leur unique, comme fait l’algorithme de résolution classique, le domaine courant est, 
plutôt, réduit à la plus longue chaîne de valeurs. Puis, un algorithme d’étoile arc- 
consistance (ac*) est exécuté et si cela ne donne pas un domaine vide, un appel 
récursif à lcpmac* est effectué dans le but de traiter les autre variables non encore 
traitées. Un point crucial qui doit être souligné ici est que la structure de la chaîne 
des domaines transformés est conservée par le filtrage du domaine. Autrement dit, 
lorsque certaines valeurs sont supprimées d’un domaine déjà traité, par l’effet de 
l’étoile arc-consistance, les valeurs restantes forment, toujours, une chaîne de paires 
de valeurs sous-modulaires. Evidemment, le domaine courant doit être réduit, à son 
tour, à des valeurs qui ne sont pas dans la chaîne, afin de préserver un algorithme 
complet. Ainsi, le deuxième appel récursif de LCPMAC*, est destiné à traiter le reste 
du domaine courant. 


6.4.2 algorithme utilisant la partition en chaîne de Dilworth 

Notre algorithme de décomposition DCPMEDAC* (pour Partition en Chaîne de 
Dilworth medac*), dont le pseudo-code est donné par l’algorithme 24, est obtenu 
en intégrant la stratégie de décomposition, dans l’algorithme standard, nommée 
medac* [76]. Pour dcpmedac*, le domaine de la variable courante est réduit, à 
tour de rôle, à chaque valeur du domaine de cette variable. Et ceci afin de préserver 
la complétude de l’algorithme. Dans ce cas, la taille de l’arbre de recherche qui sera 
exploré par dcpmedac* dépend du nombre de branches consacré à chaque nœud 
interne. 
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Algorithme 23: LCPMAC* (V,Y,s*,tt) : s* 


// avec V = (X,V,C,S) 
si Y = 0 alors 

1 s <— mincut('P) 

si COSt(s) -< COSt(s*) alors s* invert(s, - k ) 


sinon 

X <r- SELECT(y) 

2 rfb-S—PREORDER(LL) 

3 (S x , n x ) LONGCHAIN(.D x ) 

V 


4 


D x t T^x(Sx) 

V <- AC* (P 7 ) 

si 0 ^ V alors 

s* LCPMAC* (P 7 , Y \ {x}, s*, 7r) 


P 7 ^P 


D' X ^D' X \5 X 
V t- AC* (P 7 ) 
si 0 ^ P 7 alors 

S* «- LCPMAC* (P 7 , Y, S*,7r) 


Exemple 6.4.1. Rappelons que le Problème de Satisfaction Maximale de Contraintes 
(Max-CSP) est un VCSP destiné à calculer des solutions approchées pour les CSP 
sur-contraint. Lors de la résolution d’un Max-CSP, un certain nombre de contraintes 
sont autorisées à être violées. L’objectif est de réduire le nombre de contraintes vio¬ 
lées. En d’autres termes, un Max-CSP est un VCSP avec une structure de valuation 
ne contenant que deux valuations : une valuation nulle, qui est habituellement as¬ 
sociée avec les tuples consistants dans le CSP original, et une valuation de 1 pour 
les tuples incompatibles dans le CSP original. L’opérateur d’agrégation utilisé dans 
le réseau de Max-CSP est simplement la somme et l’opérateur de comparaison est 
la comparaison de nombres entiers. La structure de valuation associée au Max-CSP 
est, donc, Sm = (N, +, <)• 

Considérons le Max-CSP P défini sur quatre variables (voir Figure 6.1). Les do¬ 
maines des variables x\ et contiennent trois valeurs : a, b et c, tandis que celles de 
X 2 et X 3 contiennent seulement a et b. Un trait continu indique une incohérence entre 
les valeurs liées, puis une utilisation simultanée de ces valeurs entraîne un coût de 
1. Les valeurs non reliées sont compatibles et peuvent, par conséquent, être utilisées 
en même temps sans engendrer un coût supplémentaire. 

L’instance n’est pas dans la classe des VCSPs binaires sous-modulaires par per¬ 
mutations définis dans [108], parce qu’il n’y a pas de permutations de domaine 
qui peuvent rendre les trois fonctions connectant Xi,X 2 et X 4 simultanément sous- 
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Algorithme 24: DCPMEDAC*(P, Y, s*) : s* 

// avec P = (X,V,C,S) 
si Y = 0 alors 

1 s 4- mincut(P) 

si COSt(s) -< COSt(s*) alors s* <— inverse(s, tt) 

sinon 

X <r- SELECT(y) 

^X<r- PREORDER^z) 

G- DCP {D x ,±x) 
tant que true faire 
5 ^-select(A x ) 
v Tx <- PRM (5, ±x) 

V' 4-V 

D'x 4- 7T x (â) 

V' 4- ED AC* (P') 

si 0 ^ V alors 

S* 4- DCPMEDAC*(P , J Y \ {x}, S*) 

D x D x \ 5 

V ^ T 

V' 4- ED AC* (P') 
si 0 G P' alors 
break 
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modulaires. 

Si nous considérons les variables le long de l’ordre donnée par (aq, x 2 , x 3 , x 4 ) alors 
appliquer le schéma de décomposition proposées à cette exemple, résulte les étapes 
suivantes : 

- A la racine de l’arbre de recherche et en se référant à la Définition 40, on 
en déduit que a -< X1 b, b -< X1 c et a -< X1 c. Donc D 1 n’est pas partitionné : 
Ai = {Di} et ses valeurs sont identiquement permutées : 7r 4 = (a, b, c). 

- Ensuite, nous passons à la variable x 2 . En se référant à la Définition 40, nous ob¬ 
tenons b -< X2 a. Ainsi, a et b doivent être permutées. La permutation tt 2 = ( b , a) 
est donc appliquée pour transformer la fonction de valuation entre X\ et x 2 en 
une fonction sous-modulaire. Une situation similaire se produit pour x 3 , où D 3 
n’est pas partitionné, mais ses valeurs sont permutées (voir le graphe à droite 
de la Figure 6.1). 

- La situation est plus complexe pour x 4 . La permutation choisie pour D\ exige 
que vt 4 (6) < 7r 4 (a) afin d’avoir une fonction sous-modulaire entre X\ et x 4 , tan¬ 
dis que la permutation appliquée à D 2 et D 3 exige l’inverse. Par conséquent, 
nous 11 ’avons ni a A X4 b ni b -< XA a. Cela signifie que a et b ne peuvent pas être 
placés dans la même chaîne lors de la partition de D 4 . D’autre part, la per¬ 
mutation appliquée à D\ ne produit aucune condition sur la position relative 
de a et c comme ces deux valeurs ne sont pas connectées à une valeur de D\. 
Cependant, tt 2 et 7r 3 impliquent toutes les deux que 7r 4 (a) < 7 r 4 (c). Ceci signifie 
que a et c peuvent être placées dans la même chaîne lors de la partition de 
D 4 à condition qu’elles ne sont pas permutées. Par conséquent, D 4 doit être 
partitionné en deux chaînes : A 4 = {{a,c}, {£>}}. 

La Figure 6.2 résume ces étapes en montrant les chemins explorés par un algo¬ 
rithme qui applique notre schéma de décomposition (voir l’arbre de droite) comparé 
à ceux explorés par un algorithme classique (voir l’arbre à la gauche de la figure). 
Notons que, dans l’arbre à droite, seiil le dernier domaine est divisé en deux chaînes, 
ce qui donne deux sous-problèmes sous-modulaire. En d’autre termes, l’algorithme 
que nous avons proposé ne prend aucune décision sur les trois premières variables. 
Ceci peut l’empêcher de faire de mauvais choix. 


6.5 Résultats expérimentaux 

Les problèmes étudiés dans notre expérimentation sont, les problèmes Max-CSPs 
aléatoires binaires, les problèmes de coloration des graphes, les problèmes d’affec¬ 
tation de liens de radio-fréquences (RLFAP), les problèmes de planification et le 
problème de sélection génétique TagSNP. 
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Figure 6.1 - Un Max-CSP (gauche) et sa décomposition en sous-problèmes (arbre) sou- 
modulaires (droite). 



Figure 6.2 - Arbre de recherche du Max-CSP de la Figure 6.1 exploré par un algorithme 
classique (gauche) et Arbre de recherche exploré par un algorithme utilisant notre schéma 
de décomposition (droite). 


Dans une première expérimentation, nous comparons, l’algorithme LCPMAC* avec 
MAC*. L’algorithme arc-consistance soft employé pour ces deux algorithmes est AC* 
[73]. Les limites supérieures initiales utilisées par les deux solveurs sont ceux données 
dans les différents fichiers contenant le codage des instances. 

L’heuristique d’ordre de variables employée est min-domain/deg [9]. Pour l’ordre 
des valeurs, les deux algorithmes choisissent les valeurs ayant le coût unaire minimal 
après projection. Ces coûts sont automatiquement calculés par application de AC*. 
Les critères d’évaluations sont le nombre de nœuds développés, le temps de calcul en 
secondes et le coût de la meilleure solution trouvée. Nous avons également reporté, 
pour lcpmac*, le nombre de nœuds multiples explorés. Ces nœuds sont ceux au 
niveau desquels la variable courante a été instanciée par plus d’une valeur. 

Dans une deuxième expérimentation, nous comparons l’algorithme DCPMEDAC* 
avec MEDAC*. L’algorithme arc-consistance flexible employé pour ces deux algo- 
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rithmes est edac* [74]. Les bornes supérieures initiales utilisées par les deux solveurs 
sont ceux données dans les différents fichiers tests disponibles sur [80]. 

L’heuristique d’ordre des variables est min-domain/deg [9] combinée avec la der¬ 
nière variable causant conflit [79]. Pour celle de l’ordre des valeurs, les deux algo¬ 
rithmes choisissent en priorité les valeurs qui ont le plus petit coût unaire projeté. 
Ce coût est automatiquement calculé par la procédure edac*. Les critères d’évalua¬ 
tions sont le nombre de nœuds développés, le temps de calcul en secondes et le coût 
de la meilleure solution trouvée. Nous avons également reporté, pour dcpmedac*, 
le nombre de nœuds multiples explorés. Ces nœuds sont ceux au niveau desquels la 
variable courante a été instanciée par plus d’une valeur. 

Tous les algorithmes ont été implémentés en C++ et exécutés sur un PC à 2 
GHZ de fréquence et une mémoire vive de 4 Gb. 

Max-CSPs aléatoires binaires : Pour les Max-CSPs aléatoires binaires, nous 
avons expérimenté sur les instances décrites dans [75] et générées en fonction du 
modèle des quatre paramètres [104], L’objectif est de trouver une affectation avec 
un nombre minimum de contraintes non satisfaites pour un CSP sur-contraint. 

Le problème est facilement exprimé comme un VCSP en utilisant, pour les coûts 
binaires, un coût nul pour les paires de valeurs compatibles et 1 pour les paires 
incompatibles. 

Nous avons expérimenté sur quelques un des instances disponibles dans [80]. 
Ces instances sont obtenus en considérant trois niveaux de densité de graphe de 
contraintes : problèmes peu dense (S), dense (D) et complet (C) et deux niveaux de 
densité des contraintes : contraintes lâches (L) et contraintes dures (T). La combi¬ 
naison de ces classes donne six classes de VCSPs : SL, ST, DL, DT, CL, CT. 

Pour la première expérimentation, nous avons utilisé dix instances appartenant 
à quatre de ces classes. La plupart de ces instances sont plutôt difficiles à résoudre 
(voir le Tableau 6.1). Ces instances ont en commun un domaine de valeur uniforme 
de taille dix. Le nombre de variables, qui est indiqué immédiatement après les ini¬ 
tiales de la classe dans le nom de l’instance, varie de 25 à 50. 

Le Tableau 6.1 fournit les résultats obtenus par lcpmac* et MAC*. Quatre heures 
de temps d’exécution ont été consacrées à chaque instance. Nous avons rapporté le 
nombre de nœuds développés, le temps d’exploration de l’arbre de recherche et le 
coût de la meilleure solution trouvée. Pour ces instances non structurées, aucun des 
deux algorithmes ne surpasse systématiquement l’autre. Néanmoins, on remarque 
que LCPMAC* est plus concurrentiel sur les contraintes lâches et les instances peu 
dense. 

Pour la deuxième expérimentation, nous considérons exactement deux instances 
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de chacune des six classes décrites ci-dessus. Ces instances ont en commun un do¬ 
maine de valeur uniforme de taille dix. Le nombre de variables qui est indiqué im¬ 
médiatement après les initiales de la classe dans le nom de l’instance, varie de 15 à 80. 

Le Tableau 6.2 fournit les résultats obtenus par dcpmedac* et medac*. Une heure 
de temps d’exécution a été consacrée à chaque instance. Nous avons rapporté le 
nombre de nœuds développés, le temps d’exploration de l’arbre de recherche et le 
coût de la meilleure solution trouvée. Pour ces instances, non structurées, globale¬ 
ment aucun des deux algorithmes ne surpasse systématiquement l’autre. Néanmoins, 
sur les deux premières instances, qui sont les plus difficiles, dcpmedac* domine me¬ 
dac*. En plus, DCPMEDAC* a réussi à terminer la preuve d’optimalité de l’instance 
SL80-10-6-60-1, contrairement à MEDAC*. En ce qui concerne les instances qui sont 
faciles à résoudre medac* domine dcpmedac* . Cependant, si nous considérons le 
temps total pour résoudre toutes ces instances, nous pouvons dire que dcpmedac* 
est globalement plus rapide que medac*. 


instance 

nombre de nœuds 

nœuds-multi temps (en sec.) 

cost 

MAC* 

LCPMAC* 

LCPMAC* 

MAC* 

LCPMAC* 

MAC* 

LCPMAC* 

SL50-10-10-60-1 

11431 

6369 

1636 

- 

6809 

5 

4 

SL50-10-10-60-2 

10377 

7728 

1250 

12982 

8472 

3 

3 

SL40-10-13-60-1 

1059 

762 

142 

1078 

719 

3 

3 

SL40-10-13-60-2 

2261 

1837 

385 

2314 

1660 

3 

3 

DL30-10-25-48-1 

67 

215 

16 

64 

204 

2 

2 

DL30-10-25-48-2 

915 

181 

37 

867 

251 

2 

2 

ST25-10-21-85-1 

527 

1206 

89 

314 

770 

19 

19 

ST25-10-21-85-2 

604 

3852 

394 

370 

2321 

20 

20 

DT25-10-25-87-1 

1871 

16581 

2405 

1266 

11332 

28 

28 

DT25-10-25-87-2 

4332 

9833 

511 

2822 

7498 

30 

30 


Tableau 6.1 - Résultats obtenus par MAC* et lcpmac* sur les instances Max-CSPs aléa¬ 
toires binaires. Le nombre des nœuds et nœuds multi explorés (en milliers), le temps CPU 
en secondes et le coût de la meilleure solution trouvée sont raportés. Le temps d’exécution 
a été limité à 4 heures par instance. Un tiret indique que le programme n’a pas mis fin au 
bout de 4 heures de temps d’exécution. 

Coloration de graphe : Ce problème est une généralisation au problème de colo¬ 
ration de graphe présenté dans [68]. Etant donné un graphe simple et un nombre de 
couleurs d, le but est de colorier les nœus du graphe de telle sorte que deux nœuds 
adjacents V{ et Vj n’aient pas la même couleur. Deux nœuds sont adjacents si la 
distance est la plus petite. 

Pour les deux expérimentations étudiées, nous considérons précisément les graphes 
géométriques générés par M. Trick [115]. Les nœuds sont générés dans une grille de 
10.000 par 10.000 et sont reliés par une arête si elles sont suffisamment proches. Les 
évaluations sont inversement proportionnelles à la distance dij entre les nœuds et 


119 







Chapitre VI : Permutation de domaine pour les VCSPs binaires 


instance 

nombre de nœuds 

temps (en sec.) 

coût 

MEDAC* 

DCPMEDAC* 

MEDAC* 

DCPMEDAC* 

MEDAC* 

DCPMEDAC* 

SL80-10-6-60-1 

- 

2184 

- 

1374 

4 

4 

SL80-10-6-60-2 

2296 

1862 

1834 

1170 

2 

2 

DL30-10-25-48-1 

2.17 

10 

0.97 

4.50 

2 

2 

DL30-10-25-48-2 

21 

20 

9.41 

7.08 

2 

2 

CL25-10-100-18-1 

163 

568 

122 

333 

2 

2 

CL25-10-100-18-2 

110 

807 

79 

314 

2 

2 

ST25-10-21-85-1 

5 

11 

1.80 

3.32 

19 

19 

ST25-10-21-85-2 

1.62 

5.33 

0.69 

1.94 

20 

20 

DT25-10-25-87-1 

9.18 

11.50 

4.73 

4.74 

28 

28 

DT25-10-25-87-2 

17 

19 

7.44 

8.86 

30 

30 

CT15-10-100-93-1 

57 

78 

34.38 

43.09 

73 

73 

CT15-10-100-93-2 

72 

80 

45.40 

44.04 

74 

74 


Tableau 6.2 - Résultats obtenus par MEDAC* et dcpmedac* sur les instances Max-CSPs 
aléatoires binaires. Le nombre des nœuds (en milliers), le temps CPU en secondes et le 
coût de la meilleure solution trouvée sont raportés. Le temps d’exécution a été limité à 
une heure par instance. Un tiret indique que le programme n’a pas mis fin au bout d’une 
heure de temps d’exécution. 


elles sont uniformément générées. 

Nous avons précisément expérimenté sur neuf instances de ceux disponibles dans 
[80], quatre instances sont extraites du problème GEOM30a et notées GEOM30a-3 
à 6. Cinq autres instances sont extraites du problème GEOM40 et sont notées par 
GEOM40-2 à 6 (voir Tableau 6.3). Chacune de ces instances de GEOM30a-d im¬ 
plique 30 variables et un domaine de taille maximale de 6 valeurs. Et chacune des 
instances de GEOM40-d implique 40 variables et un domaine de taille maximale de 
6 valeurs. 

Les résultats des deux expérimentations sont résumés dans les Tableaux 6.3 et 
6.4. Nous pouvons voir que LCPMAC* et DCPMEDAC* sont plus performants que 
MAC* et medac* sur les deux instances les plus difficiles de ce problème, à savoir, 
les instances GEOM-4 et GEOM-5. lcpmac* et dcpmedac* sont, clairement, plus 
rapides que MAC* et medac* sur ces instances. En plus, lcpmac* et dcpmedac* 
ont réussi, tous les deux, à terminer la preuve d’optimalité de l’instance GEOM40-5, 
contrairement à MAC* et MEDAC*. En ce qui concerne les instances qui sont faciles 
à résoudre, MAC* et medac* ont été plus rapides que lcpmac* et dcpmedac*. Ce¬ 
pendant, si nous considérons le temps total pour résoudre tous les instances GEOM, 
nous pouvons dire que LCPMAC* et DCPMEDAC* sont globalement plus rapides que 
MAC* et MEDAC*. 

Problème d’affectation de lien de radio-fréquence (RLFAP) : Ce sont des 
instances provenant d’un problème du monde réel, à savoir le problème celar RL- 
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instance 

nombre de nœuds 

nœuds-multi 

temps (en sec.) 


cost 

MAC* 

LCPMAC* 

LCPMAC* 

MAC* 

LCPMAC* 

MAC* 

LCPMAC* 

GEOM30a-3 

321 

782 

85 

55 

133 

11 

11 

GEOM30a-4 

73 

396 

54 

18 

75 

4 

4 

GEOM30a-5 

0 

1 

<1 

1 

1 

1 

1 

GEOM30a-6 

0 

<1 

<1 

1 

1 

0 

0 

GEOM40-2 

590 

480 

60 

63 

54 

22 

22 

GEOM40-3 

759 

768 

152 

117 

127 

7 

7 

GEOM40-4 

29948 

933 

68 

6122 

199 

3 

3 

GEOM40-5 

60426 

10400 

117 

- 

2992 

1 

1 

GEOM40-6 

0 

<1 

<1 

1 

1 

0 

0 


Tableau 6.3 - Résultats obtenus par MAC* et LCPMAC* sur les instances de coloration de 
graphe. Le nombre des nœuds et nœuds multi explorés (en milliers), le temps CPU en 
secondes et le coût de la meilleure solution trouvée sont raportés. Le temps d’exécution a 
été limité à 4 heures par instance. Un tiret indique que le programme n’a pas mis fin au 
bout de 4 heures de temps d’exécution. 



nombre de nœuds 

temps (en sec.) 

coût 

instance 

MEDAC* 

DCPMEDAC* 

MEDAC* 

DCPMEDAC* 

MEDAC* 

DCPMEDAC* 

GEOM30a-3 

90 

261 

3.40 

9.26 

11 

11 

GEOM30a-4 

28 

89 

1.37 

4.04 

4 

4 

GEOM30a-5 

0.37 

0.80 

0.02 

0.04 

1 

1 

GEOM30a-6 

0.03 

0.03 

0.01 

0.01 

0 

0 

GEOM40-2 

36 

47 

1.39 

1.69 

22 

22 

GEOM40-3 

87 

216 

4.34 

10.97 

7 

7 

GEOM40-4 

28833 

597 

1222 

26.74 

3 

3 

GEOM40-5 

- 

3357 

- 

148.19 

1 

1 

GEOM40-6 

0.04 

0.04 

0.01 

0.01 

0 

0 


Tableau 6.4 - Résultats obtenus par MEDAC* et dcpmedac* sur les instances de coloration 
de graphe. Le nombre des nœuds explorés (en milliers), le temps CPU en secondes et le 
coût de la meilleure solution trouvée sont raportés. Le temps d’exécution a été limité à 
une heure par instance. Un tiret indique que le programme n’a pas mis fin au bout d’une 
heure de temps d’exécution. 
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instance 

nombre de nœuds 

nœuds-multi 

temps (en sec.) 

cost 

MAC* 

LCPMAC* 

LCPMAC* 

MAC* 

LCPMAC* 

MAC* 

LCPMAC* 

celar6-sub0 

1184 

425 

75 

1129 

499 

159 

159 

celar6-subl 

9586 

1853 

357 

8616 

3618 

2669 

2669 

celar6-sub2 

12015 

2142 

277 

- 

8550 

3163 

2746 

celar6-sub3 

13259 

5989 

1015 

- 

- 

4821 

3177 

celar6-sub4 

7443 

3530 

123 

- 

- 

6247 

5605 

celar7-sub0 

11 

1 

<1 

48 

6 

10310 

10310 

celar7-subl 

2530 

425 

74 

4439 

733 

142640 

142640 

celar7-sub2 

7836 

6391 

149 

- 

- 

193253 

212436 

celar7-sub3 

8776 

5364 

238 

- 

- 

292840 

242843 

celar7-sub4 

5759 

3936 

369 

- 

- 

383254 

292640 


Tableau 6.5 - Résultats obtenus par MAC* et lcpmac* sur les instances extraites des pro¬ 
blèmes réels : celar6 et celar7. Le nombre des nœuds et nœuds multi explorés (en milliers), 
le temps CPU en secondes et le coût de la meilleure solution trouvée sont raportés. Le 
temps d’exécution a été limité à 4 heures par instance. Un tiret indique que le programme 
n’a pas mis fin au bout de 4 heures de temps d’exécution. 

FAP [11] qui a la particularité d’être peu contraint. Le Problème d’affectation de 
fréquence radio (RLFAP) consiste à affecter des fréquences entre plusieurs ressources 
(radio) tout en assurant un minimum de perturbations subies par ces fréquences ra¬ 
dio. Ce problème est un problème d’optimisation combinatoire NP-difficile. 

Pour les deux expérimentations effectuées, nous avons précisément expérimenté 
sur dix instances de ceux disponibles dans [80], dont cinq sont extraites du problème 
celarô et notés celar6-sub0 à 4. Les cinq autres instances sont extraites du problème 
celar7 et sont notées celar7-sub0 à 4 (voir Tableau 6.5). A Chacune de ces instances 
comprend 16 variables et une taille maximale de domaine de 44 valeurs. 

Pour la première expérimentation, en examinant le Tableau 6.5, nous remar¬ 
quons que lorsque les deux algorithmes trouvent la solution optimale dans un temps 
ne dépassant pas quatre heures, lcpmac* est toujours plus rapide que MAC*. En 
outre, sur les autres instances, qui sont plus difficiles à résoudre, LCPMAC* a systé¬ 
matiquement trouvé une solution qui est, au moins, aussi bonne que celle trouvée 
par MAC* sur toutes les instances sauf une (celar7-sub2). Globalement, on peut dire 
que, contrairement aux instances aléatoires, LCPMAC* surpasse MAC* sur toutes ces 
instances sauf une. 

Pour la deuxième expérimentation, en examinant le Tableau 6.6, nous notons 
que sur les six instances résolues par les deux algorithmes, dcpmedac* est souvent 
plus rapide que medac*. Une autre instance a été seulement résolue par dcpme¬ 
dac*. Dans les instances qui sont difficiles à résoudre, DCPMEDAC* domine MEDAC* 
en terme de coût trouvé. Globalement, en comparaison avec les problèmes Max- 
CSPs aléatoires binaires, nous pouvons dire avec plus de confiance que dcpmedac* 
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instance 

nombre de nœuds 

temps (en sec.) 

coût 

MEDAC* 

DCPMEDAC* 

MEDAC* 

DCPMEDAC* 

MEDAC* 

DCPMEDAC* 

celar6-sub0 

0.250 

0.221 

0.75 

0.29 

159 

159 

celar6-subl 

481 

371 

1529 

1124 

2669 

2669 

celar6-sub2 

224 

161 

782 

595 

2746 

2746 

celar6-sub3 

- 

744 

- 

3114 

3079 

3079 

celar6-sub4 

- 

- 

- 

- 

3428 

3230 

celar7-sub0 

0.480 

0.183 

0.66 

0.24 

10310 

10310 

celar7-subl 

240 

291 

62 

71 

142640 

142640 

celar7-sub2 

522 

321 

901 

495 

173252 

173252 

celar7-sub3 

- 

- 

- 

- 

223460 

223460 

celar7-sub4 

- 

- 

- 

- 

252441 

242443 


Tableau 6.6 - Résultats obtenus par MEDAC* et dcpmedac* sur les instances extraites 
des problèmes réels : celar6 et celar7. Le nombre des nœuds explorés (en milliers), le 
temps CPU en secondes et le coût de la meilleure solution trouvée sont raportés. Le temps 
d’exécution a été limité à 4 heures par instance. Un tiret indique que le programme n’a 
pas mis fin au bout de 4 heures de temps d’exécution. 

domine MEDAC* sur ses instances provenant de problèmes réels. 

Problèmes de planification : 

Nous considérons précisément le problème de transport DriverLog. Ce problème 
implique un ensemble de conducteurs pour conduire les camions transportant des 
objets d’un emplacement de départ vers une destination finale d’une manière ef¬ 
ficace. Plus de détails sur ce problème se trouve dans [51]. Les quatorze instances 
sur lesquelles nous avons expérimenté sont modélisées comme des WCSPs de densité 
moyenne puisque environ la moitié des contraintes possibles sont présentes. Chacune 
des instances implique des centaines de variables et une taille de domaine d’au plus 
onze valeurs. 

Pour la première expérimentation, le Tableau 6.7 montre que LCPMAC* domine 
MAC* sur toutes, sauf une (driverlog04bc) des dix instances instances simulées. 

Pour la deuxième expérimentation, les résultats présentés sur le Tableau 6.8, 
montrent que DCPMEDAC* est plus performant que MEDAC* sur onze instances des 
quatorze instances simulées. En plus, les instances sur les quelles MEDAC* domine 
dcpmedac* sont faciles. 

Pour les deux expérimentations étudiées, nous pouvons alors conclure que les 
deux algorithmes (DCPMEDAC* et LCPMAC*), qui partitionnent chaque domaine en 
exploitant les fonctions sous-modulaires, sont dominants sur toutes les instances 
structurées difficiles à résoudre. 
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instance 

nombre de nœuds 

nœuds-multi temps (en sec.) 

coût 

MAC* 

LCPMAC* 

LCPMAC* 

MAC* 

LCPMAC* 

MAC* 

LCPMAC* 

driverlog04ac 

7 

1 

<1 

26 

5 

85 

85 

driverlog04bc 

116 

103 

14 

279 

391 

846 

846 

driverlog04cc 

208 

47 

7 

512 

161 

2082 

2082 

driverlog05ac 

8 

4 

<1 

63 

23 

230 

230 

driverlog05bc 

97 

106 

18 

472 

747 

806 

806 

driverlog05cc 

170 

52 

5 

703 

328 

2001 

2001 

driverlogs06 

5311 

1615 

528 

- 

6141 

1095 

1055 

driverlog08ac 

52 

16 

1 

671 

124 

735 

735 

driverlog08bc 

71 

49 

3 

851 

483 

1404 

1404 

driverlog08cc 

84 

47 

3 

892 

362 

2492 

2492 


Tableau 6.7 - Résultats obtenus par MAC* et LCPMAC* sur les instances extraites des 
problèmes de transport DriverLog. Le nombre des nœuds et nœuds multi explorés (en 
milliers), le temps CPU en secondes et le coût de la meilleure solution trouvée sont ra- 
portés. Le temps d’exécution a été limité à 4 heures par instance. Un tiret indique que le 
programme n’a pas mis fin au bout de 4 heures de temps d’exécution. 


instance 

nombre de nœuds 

temps 

(en sec.) 


coût 

MEDAC* 

DCPMEDAC* 

MEDAC* 

DCPMEDAC* 

MEDAC* 

DCPMEDAC* 

driverlog02ac 

12 

27 

54.99 

82.51 

2010 

2010 

driverlog02bc 

15 

18 

52.27 

46.78 

2085 

2085 

driverlog02cc 

6.38 

14.45 

22.63 

33.48 

2428 

2428 

driverlogs03 

2.6 

5.8 

1 

1 

1225 

1225 

driverlog04ac 

3.47 

10.46 

16.80 

20.14 

1790 

1790 

driverlog04bc 

3.65 

7.41 

12.47 

11.64 

1921 

1921 

driverlog04cc 

2.77 

2.98 

11.61 

6.13 

2932 

2932 

driverlog05ac 

3.9 

3.8 

32.60 

20.52 

1350 

1350 

driverlog05bc 

3.9 

3.5 

34.93 

14.56 

1344 

1344 

driverlog05cc 

3.9 

8.2 

33.82 

27.07 

2659 

2659 

driverlogs06 

388 

135 

821 

158 

1055 

1055 

driverlog08ac 

52.44 

29.85 

779 

327 

2220 

2220 

driverlog08bc 

53.06 

29.25 

645 

227 

2282 

2282 

driverlog08cc 

50.10 

28.52 

672 

233 

4052 

4052 


Tableau 6.8 - Résultats obtenus par medac* et dcpmedac* sur les instances extraites 
des problèmes de transport DriverLog. Le nombre des nœuds explorés (en milliers), le 
temps CPU en secondes et le coût de la meilleure solution trouvée sont raportés. Le temps 
d’exécution a été limité à une heure par instance. Un tiret indique que le programme n’a 
pas mis fin au bout d’une heure de temps d’exécution. 
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nombre de nœuds 

temps 

(en sec.) 


coût 

instance 

MEDAC* 

DCPMEDAC* 

MEDAC* 

DCPMEDAC* 

MEDAC* 

DCPMEDAC* 

10734 

228 

9.15 

185 

6.40 

150538 

150538 

11053 

54 

0.5 

232 

1.26 

140279 

140279 

12929 

- 

2.26 

- 

53.98 

1165000 

666797 

13106 

1.9 

0.11 

12.47 

1.16 

278934 

278934 

14481 

- 

59.53 

- 

178 

827112 

518137 

16864 

1.65 

0.14 

4.12 

0.21 

120356 

120356 

17289 

10.62 

0.21 

2.74 

0.04 

72638 

72638 

18827 

- 

0.47 

- 

25.87 

1611272 

920935 

19684 

0.07 

0.04 

0.13 

0.07 

96463 

96463 

20927 

8.74 

0.51 

4.20 

0.11 

82339 

82339 

21758 

0.28 

0.07 

8.26 

1.44 

243734 

243734 

22707 

- 

4.72 

- 

24.65 

287223 

287211 

25023 

0.06 

0.4 

0.10 

0.36 

132066 

132066 

26578 

- 

425 

- 

1547 

852467 

662633 

27116 

54 

2.0 

153 

3.07 

154277 

154277 

28188 

- 

0.45 

- 

3.45 

258747 

258747 


Tableau 6.9 - Résultats obtenus par MEDAC* et dcpmedac* sur les instances extraites des 
problèmes réels de TagSNP. Le nombre des nœuds explorés (en milliers), le temps CPU en 
secondes et le coût de la meilleure solution trouvée sont raportés. Le temps d’exécution a 
été limité à une heure par instance. Un tiret indique que le programme n’a pas mis fin au 
bout d’une heure de temps d’exécution. 


Le problème de sélection TagSNP est issu du domaine de la génétique. L’ob¬ 
jectif est de sélectionner un sous-ensemble restreint de polymorphisme nucléotide 
unique (SNP) qui capte l’information génétique critique des individus d’une même 
espèce. Le problème est étroitement lié à un problème, NP-difficile, celui de couver¬ 
ture par des sommets de graphes pondérés, Le problème peut être exprimé comme 
un VCSP binaire tel qu’il est présenté dans [1]. Le problème résultant est plutôt 
structuré connue le graphe sous-jasent à chacune de ses instances est construit sur 
la base d’une distance entre chaque paire de SNP. Nous expérimentons précisément 
sur plusieurs instances de tagsnp_r08 qui sont disponibles dans [2]. Le nombre de 
variables de ces instances varie de 33 à 93, les tailles de leurs domaines varient de 
15 à 74. 

Le Tableau 6.9 compare les résultats obtenus par les algorithmes dcpmedac* 
et medac*. Pour ces instances structurées, dcpmedac* domine medac* sur toutes 
les instances sauf une. L’instance pour laquelle MEDAC* domine DCPMEDAC* est 
(25023), elle est très facile à résoudre, comme elle est résolue en moins d’une seconde 
par les deux algorithmes. En plus, sur les deux premières instances, dcpmedac* est 
clairement plus rapide que MEDAC*. Finalement, DCPMEDAC* a résolu plus d’ins¬ 
tances en une heure de temps d’exécution. 
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6.6 Travaux connexes 

Dans [50], les auteurs ont proposé une théorie dont le but est d’étendre la classe des 
CSPs polynomiale max-fermés. Leur approche consiste à découvrir des permutations 
qui, appliquées à chaque domaine de valeur, produisent une instance de CSP max- 
fermé. Si de telles permutations existent, l’instance permutée est résolue en temps 
polynomial, puis une solution à l’instance d’origine est déduite par simple inversion 
des permutations. Le principal avantage de cette approche est que la recherche de 
permutations de domaine est, elle même, exprimée en tant que CSP dénommé l’ins¬ 
tance de contrainte liftée. Dans certains cas, le dernier CSP a l’advantage d’être 
polynomiale. C’est le cas pour les CSPs d’arités bornées avec un domaine booléen. 
Néanmoins, les auteurs ont prouvé que le problème de réduction dans le langage 
dont les contraintes sont max-fermées est NP-complet, même pour les CSPs binaires 
avec des domaines à trois valeurs. 

L’idée de permutation de domaine a également été utilisé par E. Chen et al en 
vue de transformer les CSPs binaires dans une classe plynomiale : la classe des CSPs 
impliquant des relations à rangée convexe connectée (CRC) [12]. Mais il est apparu 
que leur algorithme ne reconnaît qu’un sous-ensemble d’instances de CSP binaires 
qui peuvent être transformées dans la classe CRC. En fait, il a été prouvé que le 
problème d’identification de permutation qui rendent les relations de rangée convexe 
est NP-complet pour les domaines de taille quatre ou plus [50]. 

Certains algorithmes de consistance locale flexible, comme arc consistance flexible 
optimale (OSAC) [19] et arc consistance virtuel (VAC) [21], peuvent résoudre les 
VCSPs sous-modulaires obtenus par permutations de valeurs même sur des contraintes 
non binaires. Sur les VCSPs généraux, ces algorithmes agissent en tant que routines 
de filtrage de domaines où leur contribution est essentielle pour résoudre les cas de 
VCSPs durs. Outre cela, ils peuvent être utilisés conjointement avec un schéma de 
décomposition pour améliorer la résolution de problèmes, comme c’est le cas dans 
l’algorithme DCPMEDAC*. 

Dans [108], l’auteur a proposé un algorithme polynomial qui identifie les per¬ 
mutations des VCSPs binaires sous-modulaires. Contrairement aux algorithmes de 
consistance locale flexible, cet algorithme renvoie les permutations effectives qui 
rendent l’instance traitée sous-modulaire. L’auteur affirme que son algorithme ne 
fonctionne que sur des VCSPs binaires avec des valuations finis. La raison en est 
que l’une des propriétés assurant la justesse de l’algorithme d’identification, (voir 
Equation (7) de [108]), est perdue si des valuations infinies sont autorisées. Bien 
que, exprimé différemment, cette propriété indique ce qui suit : si deux couples 
de valeurs (u,v) et ( v,w ) vérifient l’inégalité sous-modulaire (l’inégalité 6.2), pour 
quelques fonctions de valuations fixes f> et une paire fixe u',v' alors ( u,w ) doit aussi 
la vérifier. Nous avons montré dans la Propriété 28 que cette propriété est vérifiée 
par des problèmes de valuations infinies à condition qu’elles soient arc-consistantes. 
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Ainsi, l’algorithme de Schlesinger peut être appliqué aux VCSPs binaires avec des 
valuations infinies, à condition qu’ils soient arc-consistants, une propriété qui peut 
être obtenue en appliquant un simple algorithme d’arc-consistance. 

6.7 Conclusion 

Dans ce chapitre, nous avons présenté un schéma de décomposition pour les 
VCSPs binaires qui tire profit de la classe des fonctions sous-modulaires, même 
lorsque le problème considéré n’est pas sous-modulaire. Le schéma proposé combine 
les techniques de partitionnement de domaine et de permutation de valeurs , deux 
techniques qui sont déjà utilisées dans le cadre de raisonnement par contraintes, 
mais de manière séparée. Le schéma proposé permet la décomposition du problème 
à résoudre en un ensemble de sous-problcmes sous-modulaires, et donc des sous- 
problèmes solubles en temps polynomial. 
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Conclusion générale 


Dans cette thèse, nous avons proposé un ensemble de techniques de décomposition 
de problème contribuant à l’amélioration de résolution de problème de satisfaction 
de contraintes valuées (VCSP). 

La première contribution (voir Chapitr 4) est un algorithme de résolution de 
VCSP qui s’appuie sur la notion de substituabilité directionnelle [92], Cette notion 
a permit d’identifier et de résoudre, en temps polynomial, des VCSPs vérifiant des 
propriétés de substituabilité. Nous avons également proposé une méthode de dé¬ 
composition de VCSP s’appuyant sur cette classe polynomiale. Cet algorithme de 
décomposition a prouvé son efficacité dans la résolution des VCSPs à Contraintes 
binaires Dures (voir Tableau 7.1). 

La deuxième contribution présentée dans cette thèse (voir Chapitre 5), consiste 
en un algorithme de résolution qui s’appuie sur la classe des fonctions dite mo¬ 
dulaires. Pour une classe de VCSPs binaires où toutes les fonctions de valuation 
sont exclusivement modulaires, nous avons proposé un algorithme d’identification 
de complexité 0(ed 2 ) et un algorithme de résolution de complexité 0(ed) [93]. Nous 
avons également proposé une méthode de décomposition de VCSP s’appuyant sur 
cette classe polynomiale. Notre méthode de décomposition a prouvé son efficacité 
dans la résolution des VCSPs binaires (voir Tableau 7.2). 

Enfin, comme troisième contribution (Chapitre 6), nous avons présenté une tech¬ 
nique de résolution de VCSP binaires qui s’appuie sur la classe des fonctions dites 
sous-modulaires [14]. Une super-classe de la classe polynomiale des VCSPs sous- 
modulaires est obtenue en calculant, en temps polynomial, des permutations des do¬ 
maines des variables d’un VCSP quelconque en vue de le transformer en VCSP sous- 


Algorithme 

MAC 

MACDS 

Domination 

10% 

90% 


Tableau 7.1 - Décomposition d’un VCSP-CbC en VCSPs-CbC Directionnellement Sub- 
tituables : sur 60 instances de Job-Shop. 
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Algorithme 

MEDAC* 

MMEDAC* 

Domination 

13.72% 

86.28% 


Tableau 7.2 - Décomposition d’un VCSP binaire en VCSPs binaires modulaires : sur 
12 instances de Max-CSPs aléatoires binaires, 9 instances de coloration de graphe, 14 
instances de Planification et 16 instanses de sélection TagSNP. 

Algorithme 

MEDAC* 

DCPMEDAC* 

Domination 

15.69% 

84.31% 


Tableau 7.3 - Décomposition d’un VCSP binaire en VCSPs binaires sous-modulaires : 
sur 12 instances de Max-CSPs aléatoires binaires, 9 instances de coloration de graphe, 14 
instances de Planification et 16 instanses de sélection TagSNP. 


modulaire. Pour les VCSPs de cette super-classe, appelée VCSPs sous-modulaires 
par permutation, nous avons proposé une solution algorithmique de complexité 
0(n 2 d A ). Une méthode de décomposition de VCSP s’appuyant sur cette classe po¬ 
lynomiale a été également été proposée. Notre méthode de décomposition a prouvé 
son efficacité pour la résolution des VCSPs binaires (voir Tableau 7.3). 


Les avantages et les inconvénients de chaque technique de décomposition présen¬ 
tée dans ce travail sont résumés dans la Figure 7.1. En se référant à ce diagramme, 
le lecteur pourra choisir la technique qui va mieux avec le problème qui est entrain 
de traiter. 


Décomposition d’un VCSP 


Dual(VCSP) = 
VCSP-CbC 

1 

Décomposition 
d’un VCSP-CbC en 
VCSPs(DS) 

Décompo 
d’un VCSF 
VCSPs(M 

sition 

J en 
ad) 

Décomposition 
d’un VCSP en 
VCSPs(SubMod) 

Limite 

Limite 


Limites 

- Le passage par le 

- Il existe une classe 

- Identification en 

dual peut aboutir à un 

plus générale. 

0(n z d 4 ) 

grand nombre de 



à chaque noeud. 

variables 

Avantaaes 

- Résolution en 0(n 3 d 3 ) 

Avantaae 

- Identification en 

(Dans le cas binaire) 

- Identification en 

Q(ed 2 ) 



0(nd 3 ) 

à chaque noeud 

Avantaae 

à chaque noeud. 

- Résolution en O(ed) 

- Une classe plus 


(Dans le cas binaire) 

générale 


Figure 7.1 - Avantages et inconvénients de chaque technique de décomposition. 
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Conclusion générale 


Une extension possible de cette thèse est de généraliser d’avantage les contribu¬ 
tions deux et trois en dérivant un schéma de décomposition similaire pour les VCSPs 
n- aires. 
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